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Berechnung von Unterschallstrémungen kompressibler Fliissigkeiten 
sine. um Profile. 
Von A. Betz und E. Krahn. 


1. Einleitung. Fir das Problem der Umstrémung eines Kérpers in einer kompressiblen 
Flissigkeit bei Unterschallgeschwindigkeiten kennt man keine exakten Lésungen. Man ist daher 
auf Naherungsverfahren angewiesen. Hier wird der Fall einer zweidimensionalen Potential- 
strémung behandelt. Ausgegangen wird von der Geschwindigkeitsverteilung in einer inkompres- 
siblen Strémung um denselben Korper. Als erste Naherung fiir die Geschwindigkeit der kom- 
pressiblen Strémung wird dann ein Ansatz gemacht, der nicht nur auf schlanke Kérper anwend- 
bar ist, und davon ausgehend eine zweite Naherung berechnet. Es interessiert dabei im wesent- 
lichen die Geschwindigkeitsverteilung an der Oberflache des umstrémten Kérpers. Nach dem 
gefundenen Verfahren ist eine Reihe vorgegebener Kérper numerisch durchgerechnet worden. 
Auf Grund der Ergebnisse lie® sich eine einfache Naherungsformel fiir die Geschwindigkeitsver- 
teilung bei iblichen Formen von Fligelprofilen aufstellen. 


_ 2, Ein bequemes Rechenverfahren. In einer nicht zusammendrickbaren F lissigkeit ist die 
Bedingung, da in ein Raumelement von gegebener GréBe ebensoviel hinein- wie hinausstromen 
mu, die sogenannte Kontinuitatsbedingung, in Vektorschreibweise durch die Gleichung 


Rae oe ad ne div 8=0 (1) 


gegeben, wobei G der Geschwindigkeitsvektor ist. Die Bedingung der Drehungsfreiheit wird 
ee oe 
iickt. Auf Grund dieser Bedingungen ist die ganze Stromung durch die Randbedingungen, 
die Anstromgeschwindigkeit im Unendlichen und die Form des umstrémten Kérpers ein- 


» der Flissigkeit ‘nicht konstant, so andert sich die Kontinuitatsgleichung in ~ 


Reet ae ESC) divgv=0; Pic (3) 


a <a } 


hung er Drehungsfreiheit bleibt unverandert bestehen: 
RAD eS. i : Fhe pe 2 3 t : 7 


Os oe age ee Ra ca) 
hnen wir die Lésung fir konstante Dichte, also fiir inkompressible Flissigkeiten 
(1) und (2)], mit » die fir kompressible Flissigkeiten [Gleichung (3) und (4)]- Um 
andbedingungen zu haben, miissen die Geschwindigkeiten im Unendlichen 8, und 
ichtu iben. Ihre Betrage kénnen verschieden sein. Wir wollen aber die Ge- 
Anstromgeschwindigkeit Eins im Unendlichen normieren und schreiben ._ 


BV =% und v/v, =w. ee od hy eae 

aN Gee Dichunastrcahoit; nicht aber die Kontinuitat erfillen: re fi 
oie he 0, din ges, So es 8 
div®4+@-gradg=W-gradg. = = 8) 


chstaben, z. B. b, Skalare mit lateinischen oder griechischen 
ktors wird entweder durch das Betragszeichen, z. B. |b| oder der — 
.v ausgedriickt. Die spater beniitzten komplexen GréBen werden _ 
elzeichen durch nterstreichen, kenntlich gemacht, z. B. v. Ihre 


sgedriickt, z.B. |v| oder vy. ly a 
s vektorielle Produkt: von %& und grad Q. ae PA 


da 
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Das Vektorfeld 938 enthalt Quellen von der Ergiebigkeit 28 - grad Q je Flacheneinheit. 
Wir kénnen die Lésung % fiir konstante Dichte auch naherungsweise fiir woe/0. setzen: 


wv BWoa/0- (9) 
Wir erfiillen dann die Kontinuitatsbedingung, aber nicht die der Drehungsfreiheit: 
Om div W=0, Om rot W/o +0. (10). 
Das Vektorfeld Wo./o enthalt ein Wirbelfeld 
Oo rot W/o =—Om W x grad 1/0. (11) - 


Der Vergleich mit bekannten Nahe- 


is rungswerten fiir die Geschwindigkeit 
cas - i sie an der Oberflache des Kreiszylinders 
Oa re und einiger tragfligelartiger Profile 
22 = A [ez zeigte, dah die Naherung w ~ W einen 
7) ree Fehler im entgegengesetzten Sinne 
30 + a \f- 76 See wie die Naherung tv &BWOn/oe ergibt 
. Ma.=0,3 J ts und daB beide Fehler etwa von der 
18 RES Sle Pia le a 4 gleichen GroSenordnung sind. In _ 
Me Abb. 1 sind die Ergebnisse fiireinen | 
16 I ina Ue : ek ene Kreiszylinder dargestellt, fiir eine | 
i . Machsche Zahl der Anstrémung 
ns : Va bs Ma,, =0,3, was einer maximalen 6rt- 
lichen Machschen Zahl Mamax =0,06 
[ | entspricht. Aufgetragen sind die 
‘ Werte von w tiber dem Kreisumfang, 
ausgedriickt durch den Bogenwinkel — 
4 ] gy vom Staupunkt an gerechnet. Als | 
richtigen Verlauf kann man etwa die _ 
ad eet Ey vierte Naherung der von Lamla! nach 
; ; einem Verfahren von Janzen und 
96 1] ; Rayleigh berechneten Werte an- ° 
i Sat sehen. Die unterste Kurve stellt die | 
O# : | Naherung w ~ 8, die oberste die | 
} Naherung w ® Wo,,/o dar. | 
G2 a 4 ! Auf Grund dieser Betrachtungen 
yp wurde von E&. Krahn? ein. Mittel- | 
- - = so = ae . weg vorgeschlagen, der eine bessere 
g RE tea ait id Cone sen Bouin wae eG Naherung ergibt. Er besteht darin, 
nas Geschwindigkeit am Rasecttec ass da8B man als Naherung . : 3) 
ce PBN Gold Dae ee 
“setat. Man erfiillt dann weder die Bedingung der Kontinuitat noch pe der Drehungsfreiheit. 
Es bleibt ein. Quellenfeld . eee | 
div % V elem = 28 - grad Veles 4 ee | 
a ein W irbelfeld eet Raran Ts yt <8,” a 
rot 3%. loaslo =—% x grad Veale oe ane (14) bi 


batchen Aber die Felder dieser Quellen und Wirbel heben sich an der Oberflache weitgehend | 
auf, stéren also das Ergebnis wesentlich weniger als die einzelnen Quellen-. und Wirbelfelder E | 
bei der Naherung » ~W baw. wx BWoa/o. ‘ : 


In Abb. 1 ist auch die = diese Weise fir den ‘Kreis ermittelte, Geschwindigheitaverteilung 
aufgetragen. Brag ee - ’ ve 

- Wenn auch die mit ‘ener einfachen Regel gewonnenen Pa xebrisss bei. den peln ccuties eB 
Verhaltnissen des Kreiszylinders noch merkliche Abweichungen von den wahren Werten vergeben, x 


so stimmen sie bei’ normalen Fligelprofilen so gut mit MeBergebnissen itberein, dafs man sich 
vielfach mit dieser Naherung begniigen kann. 


7. = 


ae henih Jahr). d. deutschen. Lauftfalitforschung S. 1165 (1939). 
ae E. Krahn, AVA-Bericht 42/A/01. 
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| Die Dichte 9 und damit auch is hangen vom Betrag der Geschwindigkeit v= v,,w ab, 
und zwar ist ae ed pe Lf pe \ DATTA) : ie Sie) ha 

{ o* Hi 2 2 ( c* ) : ( ) pir’ : ater 
wobei c* die kritische Schallgeschwindigkeit ist, das ist die Schallgeschwindigkeit, welche bei 


der Stromungsgeschwindigkeit v=c* herrscht; 9* ist die Dichte bei dieser Geschwindigkeit, _ ee ies 


und % ist das Verhaltnis der spezifischen Warmen bei. konstantem Druck und bei leadeeaivens 
Volumen. Wenn wir nun an Stelle von‘ w die erste Naherung 


w,=W oslo ; hn (16) 
setzen, so Lencen wir den Umrechnungsfaktor V ex/o /o als Funktion von VopW,/c* berechnen Sas 


und daont auch W= W, i 0/1. Fir den bequemen_ Gebrauch ist in Abb. 2 der Zusammenhang 


ae 02 EOE N OB GEN VADER AB) ERI OB ERE) a 


~ 


Abb. 2. cee a Tica ane wy der ‘Kouipreasibien’ Pearse ‘mit der’ Pee Ww 
RAGA Ile ; = : der mt Naariegs Strémung. : é 


t 


(17) 


t 
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Die Gleichungen (3) und (4) gehen daher tber in 


div (B+m’) eo =0, rot (W+w’)/|/o = 0 (20) || 


A ite eder 
Bes div Ww’ Vo= —div B Ve ; rot. w’/Vo= —rot W/o. | (21) 
Da nun allgemein fiir einen Skalar a und einen Vektor w gilt 


divaw=adivw+w- grada, pe (22) 


so ergibt sich unter Bericksichtigung, daB div %=0 und rot W=0 ist, aus (21): 


Ve 
rot ty’ = (%-tin') x Sd Ve) = — 5 (a8-+w') x grad Ing. (25) 
(i/o) 


Nun ist @ abhangig von |v|= v.o|w| und damit auch von Vin |W" | = | »| V olen: Wenn wir in 
Gleichung (15) v durch [23+ 1’ | ausdriicken, also — 


und 


hang zwischen (g/0*) und (m |%-+1'|)?.- 
Dabei ist 


Coo 


so ist m nur von der Machschen Zahl 
| Ma. =V._/¢. abhangig. Der Zusam- 
menhang zwischen m und Ma, ist in 
Abb. 3 dargestellt. 

Da sich 0, wie, sich aus- s (15) und 
(26) ergibt, als Funktion von mE ie 
darstellen. abe thes, 


ean: ait, hepyt : Ege 


_ Bedienhise Wir kénnen sie auch durch eine Potenzreihe nach ( m | X80! ie darstellen: ¥ 


Bevchitaken wir uns auf die ersten vier Glieder, so erhalten y wir eine oe ae bis 
etwa (m aaa \)?=1, 0 eee ie, Werte api iad ae 14, A, =0 As= =e (), 534. ier iiss ya 


| : ore pee = mise ae (26) 


setzen, so elaleen wir einen Auscn =e ay 


ae grad i O 7 
Abb, 3, Die Funktion m. | i ae ding d( x8 (28) 
SN oath : | Tee yer os apes 
- Wir setzen — tN St las edo eager ee ? 
Wo teme-a jetta ae (CAC Sy eae Aes ae, | 8 | 
und. finden fiir diese GroBe unter Verwendung von (15) und (26) — - 
: SRS at ape Pe ee ding _ : ms ie 
aha oak i , iN re ad (vje*)® Am Fw PTA Te} : i 
ae ape a fret ape. ao dino 1 ik et . 
ees el oh ugar As i ee deen ding a ne 
rics, ii 4 , ‘ . 14 f oe dele d(v/c* 8 ; co ( be 
rs prongs: ee a 
5 “Geb t) a lmi a a - a . 


rot aw—=arot w— wx grada, - (23) J 


erad O 1 ! © il 
div wm’ =—(W+w’) - grade =— 3 (B+ )- grad Ing ° (24) 


m= te. Se ]/te. (et) || 
Da 0a/9* und c.[e* gemaB (15) und | 
(17) selbst Funktionen von v,/e sind, — 


4 
. 


4 
af 
ca 
if 


are anaes eto 


ieee Funktion ist in Abb. 4  dartearene Sie bildet eine _wesentliche Grundlage fir ae! weiteren :f 


A= Ay + Ay (m|38-+10'])2 + dal rerariy vein pies ee (31) 
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Bis zu (m m |W+-’|)? =1,1 ist die Reihe mit vier Gliedern biauchbar, wenn wir Ay =0,264, 
A, =0,50, A, =—1,59, A, =1,89 setzen. Dabei miissen wir aber etwas groBere Abweichungen bei 
kleineren Werten (m [23 + ww’ |) “in Kauf nehmen. In Abb. 4 sind die sich aus diesen heiden Dar- 
stellungen ergebenden Naherungskurven mit eingezeichnet. Sie sind mit I und II bezeichnet. 

Unter Verwendung dieser Funktion 2 gehen (24) und (25) aber in 


div w= Am? [28-1] (W-+w’) - grad [3+ iw’| : (32) 
rot to’ = Am? |W+w'| (B+’) x grad |W+w'| . (33) 


Wenn wir tw’ gegen YW vernachlassigen, so ist uns in 
(32) und (33) von den GréBen der rechten Seite in 
der ganzen Strémungsebene alles bekannt, da ja & 
die bekannte Lésung der inkompressiblen Strémung 
und J eine gegebene Funktion von m |%| ist. Wir 
konnen daher fir jeden Punkt der Strémungsebene 
div w’ und rot tv’ berechnen. g 


Es ware nun eine sehr ighowicsine Arbeit, wenn’ 
wir die Rechnung so durchfithren wirden, daB wir die 
Quell- und Wirbelverteilungen punktweise berechnen 
und _daraus wieder punktweise die Strémungsfelder 
Penttela wollten. Wir wollen daher versuchen, diese - 
Rechnung” durch weitere Umformungen zu erleichtern. 
Dazu gehen wir von der Vektorschreibweise zur 
Darstellung der Vektoren durch komplexe GréSen 
(Blockschrift) uber. Einige Rechenregeln, die wir da- 
bei haufiger” gebrauchen, sind im Anhang zusammen- 
gestellt. ena ee 

Wir werden zwei Berechnungsverfahren zur Dar- 
tellung: bringen. Das eine beruht auf einer formel- 
maBigen Integration der vorkommenden Ausdricke, 
das andere auf einer numerischen oder graphischen 
ee Pes erste werden wir am Sapte der 


stout Roan 


ae schien n aber esha auch das erste Ver 


cs 


‘2 ae 
: ae 4. Die Funktion A. 


T i 


4 oa Gal ihrer Rechenregeln_ finbhwepedays: An-. 
r (32). ‘und (33), wenn wir tv’ gegen % vernachlassigen, 


» OW 2 aw ecpibetath GE sachs 
ee eG eis 


(mW) A 


chs 


mW) 8-4 (m W)e4 +++] (mW)? oo 


0 oe der Stromungsebene. Die iherrchenen GréBen | 
! venen. vies Anhang Gla B}fvin” : 


Far die Fae te Potentialstromung um einen * 
far das, Geschwindigkeitsverhaltnis W=V/V,, in einem 

le Achse (x-Achse, p=0) mit der Richtung der Anstrém- 
Nesee und den Koordinaten-Nullpunkt in den Kreismittel- 
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See ee SS eee 

— Yr, 2 2; 19% r 4 > : ah 
Ww? =w W=1—(=2) (c2ir 4 e-2in) 4 (2), (37) I 
| | dW _ gt _ 2 (%)\* 3:9 38) | 
oa Tee ae ee || 
dW. pdWy 2 (7938 356 fee 39) || 
eae (2) ¢ beatae BT 

—, dW 2 Tr N38 cig! To \oee tag NP aie (40 

2 Bai 1) hes Us Se Pitt Es ip "Os P\ . 
ie Zale)? 2 (i) p(t he ( ) 


Der letzte Ausdruck ist das mit m?/, zu multiplizierende Glied der rechten Seite von (34). Durch | 
Multiplikation dieses Ausdruckes mit dem fir W? nach (37) und mit m? erhalten wir den ent- 
Wins : sprechenden zu A, gehérigen Ausdruck, durch nochmalige Multiplikation mit m?* und der rechten 
Rat Seite von (37) den zu A, gehérigen usw. Die Zahl der Glieder wachst dabei stark an, weshalb- 
pees) wir hier die ausfihrliche Rechnung unterdriicken, nur ihren Gang schildern und das Endresultat 
ae (48) angeben. : an 
Rs ; Beim Einfithren von (37) und (40) in (34) treten als Faktoren Potenzen von m? auf. Diese 
aes hangen gema® (27) und Abb. 3 mit der Machschen Zahl der Anstrémung im Unendlichen der 
kompressiblen Strémung zusammen. Da 
ee’? = cos vp +isnvey (41) 5] 


ist, so erhalten wir durch Trennung des Real- und Imaginarteiles aus (34) 


div 1! = 2m? [(2)" cos 3p —2(2)' cos p + (2) cos g] + 2m* Af... Jp oe, (42) _| 
0 
“(rot to')y= 2 m? 2 [(22)" sin 3 p—2 (7)' sin g — (1) sin p]+2m* 2p. Jpen. 4B) 


Die erste dieser beiden Gleichungen stellt ein Quellenfeld, die zweite ein Wirbelfeld dar, 
bei dem auf jedem Kreise um den Nullpunkt die Verteilung durch ein trigonometrisches Polynom — 
ausgedrickt ist, das aus Gliedern von der Form cos yp bzw. sin vp besteht. Diese Darstellung 
hat das Angenehme, daB wir das zugehérige Geschwindigkeitsfeld ohne lange Rechnung angeben 
kénnen. Sind auf einem Kreise vom Halbmesser r Quellen mit der Ergiebigkeit ¢cos vp je 

_ Langeneinheit verteilt, wobei ¢ konstant sein soll, so ergeben diese in einem Punkte r, et PL 
ein Feld mit der Geschwindigkeit w,, fir das die einfache Beziehung. gilt 


| =1 7 De =1) ., 
. = = 5 ( + y° (eas YP, — 1 sin YY) = — s aoe Gen. . (44) 
ef. ‘ a 
- wenn r, <r ist, bzw. a 
1 We Tae Beale ein (45) 
eo Pi 2 Ty , : | 


fir r, > r. Dabei stellt. wegen der Division mit e'% der Realteil der rechten Seite, also das Glied a 
mit. cos, die Radialkomponente und der Imaginarteil, das sin-Glied, die Tangentialkomponente 
“von We: dar. © -«, i ORE iersiee Skee 
Sind auf dem Kreise vom Halbmesser r Wirbel mit der Zirkulation y sin vp je Langeneinheit | 
verteilt, so ist die Geschwindigkeit w, des zugehérigen Feldes im Punkte r, e'% gegeben durch | 


h ven y—1 ‘ es + 

BER ae 

er . ’ ' ” a de 

> . : i pei ei x ‘ . Ns ‘ 
? aif, ee (YY Plone eben eet i (47). 
| ean (4% 


S — 
‘| 


_ Dabei stellen wieder Real- und Imaginarteil die Radial- und Tangentialkomponente von w, dar. — 
_ Wie man durch Vergleich von (44) und (46) erkennt, heben sich fir r, <r das Quellenfeld | 
und das Wirbelfeld gerade auf, wenn ¢=y ist. Fiir ry >r heben sie sich auf, wenn ¢ =—y ist. | 
Da das gesamte Strémungsgebiet und damit das ganze Wirbelfeld auBerhalb des umstrémten 

_ Kreises liegt, also stets Ty <r ist, so heben sich bei der Berechnung der Geschwindigkeit auf Al 
der Zylinderoberflache, also auf dem Kreise mit dem Halbmesser ry alle Anteile der Quellen- | 
und Wirbelfelder auf, bei denen in (42) und (43) die einzelnen Glieder gleiche Vorzeichen haben, =| 
_ da ja ihre Betrage ohnehin gleich sind. ~~ ae js sass a | 
Die Verteilung der Quellen und Wirbel nach (42) und (43) setzt sich aus Gliedern von der | 
Form (ro/r)“ cos ym baw. (ro/r)“ sin yp zusammen, wobei (6 =3 ist. Das Feld einer solchen 
Quellen- oder Wirbelverteilung erhalten wir durch Integration wiber r aus den Ausdriicken (44) 4 
bis (47), wobei fiir ¢ baw. y jeweils (ry/r)* einzusetzen ist. Da hierbei der Integrand verschieden | 


1 
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ist, je nachdem r < oder > r, ist, so miissen w 


ir die Integrale von To (Zylinderoberflache) bis r, 
und von r, bis o getrennt berechnen. 


Fuhren wir die "eee in dieser Weise durch, so er- 
ies wir schlieBlich folgende Formel: 
ate af Hp OU 34 Y9 aoe To ‘le “al ro) ein (22) to] 
e g = tam?}| (2) (2) Je 72 |(72) ce) ee 3 (2) e le Pas 


+ Ay m4 ie: (3) Jerre [(2) (22)"] este 3((2) (22)"] etn 4 


Ty 


ole) edthealee |, ol elena ore 


ee eee ayers 
isl 
ee i) Bile 


Gyr (Bfer] +d yon ens 
(2) Retin — (2a) aie] 2 (20) o-tiy oe 
+ dame|— [(2)'-(2)"Jorr + 9[h (ys | 
AG ey esns [ey _myenn 
a. ne 15 | (24)°— (2)'Jeer— a0 
we - ae AG pee ley) eto 
pees ds a7 Ta, 2 wt : 
ee ee 
2 


5 [ (To 

7 1 
Glew ele] 
a 


Blo —— 
emer A 
iS 
an 
|e 
J be 

o 
as 
a 
Pt 
BH lo 
So 
(eee 
ts) 
_ 
~ 
Ss 


10| (22) as ( 4 ‘ e | 
Bena o aes 
; rRaeeRe = , Bint 
3S (20) aig Gye ‘| oa eee 
4 hry 7 9 Ty. 4 Ty 3 


ny ne eon ie foe T,\4 3 y. 

{io > 1 To \8 Bulg PHL ed 4 3i¢ to eR gur (=) e ae 

cao (eyeio] : 

GEA ocd WE USS US AeA Gaki e T 
3 ope eae (to 4 ee Se = on8¥e. 2 | Te * Sip as 
(2) er ag e / rid Ss 

fl sieceieeee NLO : 
OC a DS ia aera eee mer “)- : 
a =e pes ime ae 


i 


9 2 Sangam) ae 
He ps ee 
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Als zweite Naherung erhalten wir 


wi tw. = (W+ w’) | e0/e- (50 


Das Ergebnis ist in Abb. 5 fiir den | 
Kreiszylinder dargestellt. Es stimmt | 
ziemlich gut mit den von Lamla mit — 
vier Naherungen nach dem Verfahren 


von Janzen- Rayleigh berechneten Wer- 
ten uberein. 

Um auch ein Bild zu cuhalien: wie |} 
sich die Felder im Inneren der Flis- 
sigkeit verhalten, sind fur die Winkel 


aI 


g=0', p=} und EN le die Fehler 
w |/00/o fiir verschiedene Radien r, 
berechnet und nach GroBe und Rich- 
tung in Abb. 6 dargestellt. 


Man kénnte das geschilderte Be- 


aeeennes inkompr. 

—— L Naherung 

——— If. Naherung 
°- nach Lamla 


°) Naherungswerte nach Formel &¢ 


rechnungsverfahren grundsatzlich auch 
bei beliebigen Profilen anwenden, in- 
dem man diese jeweils konform auf — 
einen Kreis abbildet und die trigono- 
. metrischen Polynome in dieser Kreis-— 


0 Qt 02 O83 OF 
us 1 ie Abb. 5. Naherungswerte beim Kreiszylinder. % 


wierig ‘und auch. unibersichtlich. Man 
_wendet daher. hierfiir besser das im fol- 
pcpden -dargestellte Mecegs an. 


5. Strémung um ein Fliigelprofil ohne 
Zirkulation. Tn (34) 
Ow’ 1 


aig. dW 
2 
W) iz 
z abhangige Funktion. Wenn wir tiber z 


-integrieren, geht sie ‘unverandert vor ‘das 
eee 


tentialstromung, die ebenfalls nur von z 
und nicht von,z abhangt und daher die 
Rolle einer Integrationskonstanten spielt, 
indem dU/dz=0 ist. U wird am SchluB 
so bestimmt, daB die Randbedingungen 
erfallt sind. ; 

Dabei ist 2 eine F unktion von |W? 
und da |Wi?= -W W ist, hangt es sowohl 
x von Z wie von z ab. Da wir tiber z zu inte- 
; grieren haben, miussen wir hi in Funktionen 


D 


‘ist dW/dz eine nur von Z, nicht aber von 


ain ist. hel tite  Geschonindigkei « einer Po- 


Abb. 6. Grobe get Richtung ms Rowektie wv’ | Coo im Pelde eines 


ebene berechnet. Da hierbei aber zu-~ 
den bisher aufgetretenen Funktionen 
auch noch die Abbildungsfunktion hin- 
zutritt, wird die Zahl der Glieder der — 
trigonometrischen Polynome sehr groB 
und damit die Rechnung sehr lang- | 


waa 


002 004 


Kreiszylinders bei Mag =0, 36. ri eet 
Die Kurven eehon den Betrag der Korrektur im » Deigegsbonen: MaBstab an, 
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von z und Z aes da wir nur bei einer nicht von z abhangigen Funktion von z die In- 
tegration ausfihren kénnen und nur nicht von z abhingige Funktionen von Z bei der Integration 
als Konstante zu behandeln sind. Wenn wir fir A die Potenzreihe in W2 wie in (34) einsetzen 
und beriicksichtigen, dah W?"—W”" W” ist, so haben wir in jedem Glied die gewtinschte Tren- 
nung in eine Funktion W”, die nur von z, und eine Funktion W”, die nur von z abhangt. Dem- 
nach wird nach (51), wenn wir nur vier Glieder Te cpHery 
v4 


te 

wes hi, me [ Wedz + Aamtw [ Weds ay mo we f Weaz 4 
Z9 s . ; ee : 29 (52) 

eee | he 4+ U(2). 
Z9 

- Anstatt die Potentialstrémung U fir den gesamten Ausdruck zu ermitteln, kénnen wir sie 
auch auf die einzelnen mit An multiplizierten Glieder aufteilen und far jedes so bestimmen, daB 
die Teilstroémung tangential am Profil verlauft. Wir miissen dann diese Potentialstrémung zwar 
so oft bestimmen als wir Glieder angesetzt haben, also viermal. Wir haben aber den Vorteil, da® 
wir dann aus diesen Gliedern die Lésung fiir alle Machschen Zahlen, die sich im Werte von m 
ausdriicken (Abb. 3), leicht angeben anne: Sonst miBten wir fir die Gesamtfunktion die Po- 
tentialstromung U fiir jede Machsche Zahl neu et) was im allgemeinen mehr Arbeit macht. 


i 
“3 


Bp mcemae gesegen wir w’ in Teilbetrage 


: h re a ‘ ; Si ‘ ; _ . Z9 , 
ba Is Aasgangepinky Zo fair die Adtegeale wahlen wir ir zweckmabig die F higélhintorkante und 
als Integrationsweg zunachst die Profiloberflache, die wir entgegen dem Uhrzeigersinn umfahren. 


(55) 


y 


se on aes an ‘das Profil mit eae x-Achse (Abb. 7). Das aha Zeichen 


sich auf die Oberseite, ae auf 


; Abb. ce Fogo 4 der Beceichmungen. 


ETS, 


Ryle Aine? wy + A, mt ws + dam? ws + Asm’ Wis (53) 
bei eR La tee eet : 
f Cy : sein Key é tay wi, be La aM we f Wer) dz + oe ) | (54) 


Die ee ete von. der Paper onte bis zu dem Beeee betrachteten Punkt P bezeichnen 
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Wir kénnen die zu integrierenden Werte anstatt langs der Abwicklung des Profilumfanges 
auch langs seiner Abbildung auf den Kreis der c Ebene durch z =f(¢) auftragen. Der Halbmesser 


dieses Kreises sei a, sein Rand durch ¢ =a ei? gegeben. Wir setzen also 


ds=alf'(0)| dq, 5 eS EE Ge) 
beachten, da® fiir die Geschwindigkeit 
ae ‘, 2 sin p (59) 
ci betes) 
gilt und erhalten anstatt (57) _ 
74 v2 : d 
Wes) az af Wt) dz =4 ai Z sinn?d (—- W)" singdp = 2i Jn. (60) 


29 Zo 


Dirck die Integration erhalten wir die Funktion J, abhangig von der Bogenlange s baw. 
dem ihr bei der konformen Abbildung des Profiles auf einen Kreis entsprechenden Winkel Q | 
des Kreises fir jeden Punkt der Profiloberflache. Wir miissen nun diese Werte noch mit 


dW w tered) OW Wy 08 (n—1) pi(n+1)0 (= 1)" 61 | 
Fal or riW SO) (F}) | Oy) 

multiplizieren. Dabei_ ist 
Oo (62) |} 
ds R 

der reziproke Kriimmungsradius der Profilkurve an der betreffenden Stelle, und 
ad OW (63) a 

ds on 


ist der Bt ah suk aie yon W in der es dn 1 ds (Abb. 7). Real- und Imaginarteil ’|f 
~ des Ausdruckes 


(w,—U,) ec! = fn We pause Saat Wee $2) i IF (64), 


stellen die Tanbentials und Normalkomponente Tr, und Nn a Geschwindigkeit | iW, Us auf 
der Profiloberflache dar. Diese ergeben sich zu 


aw. W | ee 
Th =- we- J, (2 ys 5m no ae Re 69s nd) (+ 1)", (65) 
Na= We In (2™ cos nt — sin nd) (Fy. (66) 


Um- fir die Gecchausieler w, die Normalkomponente zum Verschwinden zu bringen, 
miussen wir die Potentialstromung Un so bestimmen, daB ihre Normalkomponente gleich —N, 
wird. Durch diese Forderung und die Bedingung, daB die Strémung frei von Singularitaten sein 
und im Unendlichen verschwinden muf, ist U, eindeutig bestimmt. Die Berechnung dieser 
‘Potentialstrémung i ist nach bekannten Verfahren leicht méglich. Sie wird noch dadurch erleichtert, | 

daB von der Berechnung der Potentialstromung Ww, 
von der ja das ganze Rechenverfahren. ausgeht, in der 

_ Regel ohnehin schon Hilfsmittel vorliegen, wie z. B. 
‘die konforme Abbildung des Profiles auf den Kreis. 
Zur Erhéhung der payee Se empfiehlt es sich, — 

: - von dem Integral J, (60) noch eine Konstante pire | 
Abb. 8. (aresbneeata etoe aa  Stromung, ziehen. ZweckmaBig wahlt man dafiir den Wert von Jn 

: im Punkte gy =z. 3 

Im Aligeatainer ist das peaktische Bedirfnis mit der Berechnung der Ge anotdice ee baw. . 
Druckverteilung lings der Profiloberflache befriedigt. Wenn aber in manchen Fallen auch die 
 Geschwindigkeit im Innern der Flussigkeit. interessiert, so ist die Berechnung ebenfalls ohne — 
- besondere: Schwierigkeit méglich. Wir wahlen dazu den Integrationsweg lings der Oberflache, 
also langs einer Stromlinie, bis zu einem Punkte P, der gleiches Potential hat wie der Punkt P, 
_dessen Geschwindigkeit berechnet werden soll, und setzen dann den. Integrationsweg langs der be- 
treffenden Potentiallinie fort (Abb. 8). Fir den ersten Teil des Integrationsweges gilt das bereits ~ 
geschilderte Verfahren. Bei der Integration langs der Potentiallinie ist wie bisher W = = We-id- 
und dz=—ie'? ds’, wobei ve die positive bzw, it ha Hiobcune der Stromung ney und 
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ds’ das Wegelbnent: langs des Integrationsweges, also —_ langs der Potentiallumie, ist. Es wird 
also entsprechend AS2), Neal : 


Jf Weds f (wy in mop ees nd?) ds’. (67) ‘ ee 


An Stelle des cnesprockiadon Faktors in (61) tritt wegen der Differentiation langs der Potential- ‘< 
linie : pepe Sse 


r) 


tea dW ot a AW. 

. BU ar ra 1 rae 

ie a io (i Qs! R) . oe a ee eae 

--wobei R’ bata Krimmungsradius der Potentiallinie und ae eran 
oo 0) 


¥ d. h. der negative Diffeontialquotient von W senkrecht zur r Richtung des Integrationsweges ist. 


6. Strémung um ein 1 Fliigelprofil mit. Zirkulation. Wird ¢ ein Kreis um den Nullpunkt. mit den pee 
Halbmesser a unter einem Winkel « angestrémt und gefordert, dak der hintere Staupunkt un- 
be verandert im Schnittpunkt des Kreises mit. der x-Achse bleiben soll, ‘so mu der Strémung one Se ie 
4 nach (36) eine Zirkulationsstromung mit der Zirkulation 


BOSE ANR BEAD actA ot ite door atthe : Fe he (10). = Te 
eae werden, die auch bei konformer ap tane auf das ‘Biapelprofl unter Beibehaltung eal 
des” Unendlichen erhalten bleibt. © ce 


z Die Geschwindigkeit im Aufengebiet des Kreises bei Heit Sons unter dem Winkel « und < ae 
der Zirkulation Ae aoe mEOy cS aaa : aes ee 


‘ P. 
’ = 


1 


me gt eae — i(a/é) 2sin a. . eS ; (a) 
i i poe - + , 


Wie] = thang —a) + sin a eae (72). 
Die ie Geschwinigheit am -F Fligelprof erecatarr sich hieraus mites der Abbildungsfunktion 


: cht” mehr adenees sondern ean beliebig viele, ‘sich 
von. We P unterscheidende Werte annehmen. In jenteprechender, eee sind dann 
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7 


Bei der Integration ist nur das Glied mit der Potenz 1/¢ mehrdeutig. Es ergibt 
| ode ee eal a3 c ! 

2(n+1) tae—'"* sin a f a = 2 (n+1) tae*"% sin la : Been i i 

: Um die stérende Mehrdeutigkeit dieses . 
Gliedes zum Verschwinden zu bringen, — 
addieren wir das bezitiglich ¢ dazu kon- 


jugiert komplexe 


ao, 
2(n+1l)iae—***sin « n—, 


das die gleiche Mehrdeutigkeit mit 
Maco dem umgekehrten Vorzeichen hat. Es 


Rags inkompr. ist eine Funktion von z baw. C, kann 
—— I. Naherung ; 


r » -. daher bei der Integration uber ¢€ als 
¢ = Naherungswerte nach Forme! 8¢ 


=07 


$A BR AF ORIN cl et BIL FE 


Integrationskonstante hinzugefigt und 

durch entsprechende Anderung der noch 

offenen Funktion U ausgeglichen wer- | 
_den. In analoger Weise ergibt die nach 

(57) hinzugefigte Funktion f/We+) dz } 


den mehrdeutigen Anteil 3 
0 6 here 98, OF 95 Pee ay 


4 — 


peso FY santana it 


Abb, e Ney mx am m elliptischen Zylinder vom Dickenverhaltnis 0,2. ae (n+ 1) j 


Dieser ist selbst schon eine Funktion von a und kann daher durch Subtraktion beseitigt werden. | | 
Wenn wir noch E 
“€=re? bzw. in 


= a 


= In—— —1ip . gy ts Sots 


setzen, erhalten wir an Stelle von (57) bow. (60) die er ORO Gleichung 


“ol 


fs ee f Wer) dz — f wendz eae tae a(e—?n* — e'"2) sin « In 
Zo Z0 f 


_ Bei der Integration langs des. Kreis- 
randes ee vad wird — 


af 


IE 


2a flew "sin nd x 


x [sin (p— a) a dees 
Res 2(n-41) ay sin na sina.) 


i UD weitere Rechnung- ‘ist ‘Mie weithe 
wie beim Flagel ohne Zirkulation. 


aa dm allgemeinen wird die so berech- | 
nete. -Korrekturstrémung. eine Umstr6- : 
_mung_ der Hinterkante ergeben, und 
nae es diese bei der i caeuiay 


=< (Pkompr 
EES G Naherung — 


— ” 


Ue bein Or, ie Kasei SENS Foi Et tila ihe Naherungswerte 
er Zirkulation anbringen. -Ergibt die _ ‘ast4— ree: tas we 
-Potentialstrémung_ in der ¢-Ebene in ‘ ait) se ere can | 
dem der Hinterkante entsprechenden 
Punkte | eine Geschwindigkeit Uy so ist 
eine Zirkulation . 

cats Ta-2naUy @ 


om _uberlager, er tease 


Abb. st. _ Goshen am oelliptscen iat vom Diekenverhitais 06. 


Poy 
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7. Eine praktische Rechenformel. In Abb. 9 bis 14 sind die nach ten geschilderten Verfahren 
berechneten Geschwindigkeitsverteilungen fiir einige Ellipsen und einige Fligelprofile dar- 


ij gestellt, und zwar jeweils fir die inkompressible Strémung (Ma <1) und die erste und zweite 
Naherung fiir die kompressible 


‘§ Stré6mung bei einer Machschen 
8 Zahl, die eime maximale 6rtliche 
‘i Geschwindigkeit ergibt, die nur 
wenig unterhalb der Schall- 
geschwindigkeit liegt. Man er- 
‘sieht aus den Kurven, da® die 
Differenz der zweiten und ersten 
"| Naherung im allgemeinen ent- 
a gegengesetzt gerichtet und klei- 
‘inj mer ist als die Differenz der 
\@ersten und nullten Naherung 


ge /nkompr. Fi 
| (wenn man als solche die Ge- eta. I Naherung 
i schwindigkeit derinkompressib- 4 ya fe 


len Strémung bezeichnet). , es pa ga 


mM 8 Man ersieht aber auch, daB 
Tr | die durch die zweite Naherung 
@ gegeniber der ersten gegebene 
_Korrektur einen ziemlich ‘ein- 0 DIR AOR LOTR TAOP ST OS DOG OTe OS ROOT TO 
fachen Verlauf iiber die Fligel- 
)tiefe hat und nirgends besonders 
‘groBe Werte erreicht. Dies legt den Gedanken nahe, eine Formel zu suchen, welche die be-. 
rechneten Korrekturwerte einigermaBen gut wiedergibt, aber nur von einigen wenigen leicht 
om _ feststellbaren. GroBen abhangt. 


Abb. 11, Geschwindigkeit am Profil NACA 230—12. 


Zunachst lieB sich an Hand 
— der berechneten Strémungen um 
Ellipsen fir die Maximalwerte 
der Geschwindigkeit eine Be- 
ziehung finden, welche nur die 

— .Machsche Zahl der Anstrémung 
und die értliche inkompressible. 
Geschwindigkeit enthalt. Fir 


AS Ue aie pigeon die Punkte auBerhalb der Sym- 
ff -metrieachse ist noch ein EinfluB 
- des Geschwindigkeitsgradienten 

. _dW/dx zu bericksichtigen. Fur 

: _ Profile, bei denen die Geschwin- 
| (seeseses Anhompn ne digkeits verteilung unsymme- 
ST Naherung aN trisch zur Profilmitte ist, tritt 

4 nae Niheriagsivarhe hierzu noch ein Faktor, der 

| 1. nach Formel 84 von der Lage des Geschwindig- 

Sasa paces keitsmaximums abhangt. Die 

f _ Formel fir die endgiltige Ge- 


schwindigkeit w, die auf die- 
sen Parametern aufgebaut und 
: nai Ke durch Prohieren den genauen — 
“Abb. 125, Coschwindigkeit ary Profil NACA 0009--1,130. Werten méglichst gut angepabt 
Be Agadir i y Bead ist, lautet. aay SH 


TC a SR AE YR 


- OCeye cae ee gy) He gy i 8 as AW : | 
4 bs 7 a) Caer pe et) ee ae) A ” ——__—, j 
[w2—1]? Mats — 0.4 at )( 2 [) Mate Wey OM): 
Berechnung-der Strémung ual Ellipsen und der Folgerungen daraus 
E Krahn erloleene ent ah in ae : 

Bn Sy heim Reker ey Fe toe picts *y aa KN 2, on : ae 
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400 


al Sins WAG! Gt oe inkompr. 
_ —— LNaherung 
: — 7 ” 
950 ¢  Naherungswerte 4 Sy 


nach Fortne/ 8¥ 


G3 


02 


MQ oo = 9,55 
Petes 

pose at seca /nkompr 

: -—— £ Naherung 
ey ”- . 
¢ WNaherungswerte 
nach Formel 84 


06 


Abb. 14, Geschwindigkeit im Profil NACA 0030—1,130. 


Abb. 15. Lage des Geschwindigkeitsmaximums. ya 


H > a 4 = 
/ : fim 


_ trigonometrische Polynome aus- 
driickt, deren Felder leicht an- 


das. Flachenintegral auf ein Rand-— 


‘ziemlich - umstandliche Berech- a 


-ersparen, wurde auf Grund der | 


_ gegeben [Gleichung (84)], welche 
_ die berechneten Werte gut 
_ wiedergibt. Se eee 


eee : 
Ingenieur-Archiv ~ 


Hierin bedeutet 1 die Profiltiefe, 
— die Entfernung des Geschwin- 
digkeitsmaximums von der Vor- 
derkante und x die laufende Ab- 
szisse (Abb, 15). 


In Abb. 5 und 9 bis 14 sind 
einige nach dieser Formel .be- 
rechnete Werte eingetragen. Sie 
zeigen eine recht gute Uberein- 
stimmung mit der zweiten Nahe- 
rung. 


8. Zusammenfassung. Zur Be- 
rechnung. der .Strémung einer 
kompressiblen Flissigkeit mit der 
Geschwindigkeit » kann man von 
der entsprechenden Geschwin- 
digkeit & der inkompressiblen 
Stroémung ausgehen und in erster 


Naherung »/v. = (¥/V.o) Veale | 


setzen. Eine zweite Naherung Jj 


ergibt. sich aus dem Felde von 
Quellen und Wirbeln, die bei der © 
ersten Naherung vernachlassigt 
sind, Die zur Berechnung dieses — 
Feldes nétige Integration uber 
die ganze Strémungsebene kann 
wesentlich vereinfacht werden, 
indem man entweder in einer 
Ebene, in welcher der umstrémte 
Rand ein Kreis ist, die Belegung 
der Quellen und Wirbel durch 


zugeben sind, oder indem man 


integral zuriickfihrt. Nach der 
ersten Methode wurde die Stré- 
mung um-einen Kreiszylinder, 
nach der zweiten die um Ellip-— 
sen und F lagelprofile berechnet. 
Um fiir die Praxis die immerhin 


nung der zweiten Naherung zu 
berechneten  Korrekturen eine 


praktische Rechenformel zur Ab- 
schaitzung dieser Korrektur an- 


aN Re ; 
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Wy Be ier eek Anhang. 


- Kinige Vektoroperationen i in Koniploser Schreibweise’. Es sei a ein Vektor mit den Kompo- 
-henten a, und a, im rechtwinkligen (x, y)-Koordinatensystem. Entsprechendes gilt fir 6. Die 
BE pedeviang der Zeichen - , x, div, rot, grad ist durch folgendeVorschriften festgelegt, wobei (') 
ii \/eeediezur (x, y)-Ebene normale Komponente bedeutet: 


* 


N 


lb Usa . ime a-b = a,b, + a,b, (ax b)y = a,b,—a5b,, : (1) 
< 4 diva Sp Se,  - (rota)y— 2 — oe, () 
i Der Ausdruck grad |a| hat die Komponenten ue und a 

Wir fihren noch folgende Bezeichnungen. ein: 
ma etek Niele a, ta, =a a, —ta,='4, (3) 
fees wo i die.i imaginare Einheit ist. coe ist die Bedeutung von b und b. Fur den Betrag 


“des Vektors a gilt dann ‘ % 
tl Yon ’ a] = =|/aa. a a * ; (4) 
iin |e Wie man leicht nachrechnen kann, gelten ae ee . 
itn ; a: b+i(axb)y=ab und a-b —i(axb)y=ab. (5) 


An Stelle von x, y wollen wir durch folgende Transformation die a A Groen z und 
z etniobreh « 


i f RiySz,  wetye Ze. ete (6) 
at as ee Zee . 
EE EEO | i ERS RO Ace Oar 3 (7) 
Sos jean ye meer a Se kaa 
paz le ek | Foe aera Bee oe! 


Fur diva und rot a gelten folgende soe ye 


ss diy afi (rot a) = 255 - : div. a—t (rot On= = (9) 
Die dem eat eateprechende komplexe GroBe TaBe : sich so darstellen 
ene I 
. Oye “Afal 8 ial, a fal i h@a > | (aa), ap 
wes EGS ER Ree pes et ee at by ey ee ays Bene Cw Gin con { 
| 1 (5) und (10) Jassen sich nun ae Formeln ableiten: . sj 
ind Dees We aba tO) 
Sige O1Dio. | Cae cat bbyse: “pericytes 
Spe sad he i (a gad is = ] ae Sree Ee SEALS Bes 


der BiSdiee coindighers des pe upretilen Cutes mit den Kom- 


b | u und wv. Jin na und bi sRichtung - ist, und beachten, daB die komplexe | GroBe i 
- (13) 
VEN i a 
(15) 


| Ibert Bets, (20) Gottingen, Hotchesses LanantieBe 39a. oe x 
srahn, 2 Zs Z, 224 Uieidest Road, Henpten Hill eee eee 
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Allgemeine Theorie der Stabilitat erzwungener Schwingungen 
elastischer Korper. 


Herrn R. Grammel zum 60. Geburtstag am 3. Marz 1949 gewidmet. 
Von E. Mettler. 


1. Einleitung. In einer vor kurzem in dieser Zeitschrift veréffentlichten Arbeit’ (im folgenden 

bezeichnet als ,,erste Mitteilung“) habe ich mit der Grundlegung einer allgemeinen Theorie 

: der Stabilitat der elastischen Bewegung begonnen. Diese vorausgegangene erste Mitteilung hatte 
den Sinn einer Problemstellung. Sie warf die Frage auf, ob die durch zeitlich veranderliche Krafte 


Mitteilung zu unserem Thema beschaftigt sich mit der praktischen Beantwortung der gestellten 
Frage, und zwar fast in derselben Allgemeinheit, in der das Problem damals mathematisch 
formuliert wurde. Ich mache im folgenden’ keinerlei weitere Einschrankungen hinsichtlich der 
Art des Kraftangriffs oder hinsichtlich der Form des elastischen Korpers, sondern setze lediglich 
von jetzt ab voraus, daf der zeitliche Verlauf der auBeren Belastung der einer harmonischen 
Schwingung ist. Da man fir den praktischen Ausbau einer Integrationstheorie der Stabilitats- 
gleichungen bestimmte Annahmen uber die Zeitabhangigkeit der Lasten machen mufB, so er- 


nommene Voraussetzung einzufiihren, mit anderen Worten die Stabilitat der erzwungenen 
harmonischen Schwingungen elastischer Korper zu untersuchen. 


vorgefiihrte Theorie arbeitet namlich zum Teil mit Hilfsmitteln und Ableitungen, deren mathe- 
matische Zuverlassigkeit nicht von vornherein feststeht. Es ist notwendig, dies hier deutlich 
-auszusprechen und den Sinn, der den nachstehenden Ausfihrungen beigelegt werden soll, klar 
anzugeben. ; 


Der grundlegende erste Schritt der folgenden Theorie, namlich der. Ubergang von der das 
Stabilitatsproblem formulierenden Variationsgleichung zu dem Differentialgleichungssystem, an 
dem die praktische Rechnung ansetzt (I. Abschnitt), benutzt eine ,,direkte Methode‘t der Va- 
riationsrechnung. Es ist eine Naherungsmethode, deren Konvergenz und Brauchbarkeit vorerst 


nicht streng bewiesen wird. Ferner werden im Verlauf der zahlenmaBigen Festlegung der In- 
stabilitatsbereiche (III. Abschnitt) unendliche Reihen aufgestellt, uber deren Konvergenz nichts — 
bekannt ist, ja bei denen man sogar mit Semikonvergenz rechnen muf. Zwar kennt man auch 
hierzu verwandte Entwicklungen — namlich aus der Stérungstheorie der Himmelsmechanik —, 


liegenden Anwendungsfall missen solche Erfahrungen doch erst bestatigt werden. 


genauen zahlenmaBigen Lésung des Stabilitatsproblems gelangen kann. Im gleichen Sinne ist . 
_ etwa eine allgemeine Formulierung des Ritzschen Verfahrens zur Lésung von statischen Stabili- 
tatsproblemen nicht als Existenzheweis fur die Verzweigungspunkte des elastischen Gleich- 
_ gewichts aufzufassen, sondern als Aufstellung einer allgemeinen Rechenanweisung. Im Einzelfall 
_ erst kann aber die numerische Brauchbarkeit des Verfahrens und die erzielbare Genauigkeit der 


Cee eens see 


x 


-durchzufihren, ist eine notwendige Erganzung, durch welche die vorliegende Theorie ihre eigent-. 


_weshalb ich mich hier auf die allgemeine Theorie und ein einziges, eng begrenztes Zahlenbeispiel 
beschranke. . ; may Oe RS ON Belt mages: Ue EERE OM 


f 
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erzwungene Bewegung eines elastischen Kérpers stabil ist oder nicht. Die vorliegende zweite | 


scheint es besonders zweckmafig, zuerst die genannte, der traditionellen Elastokinetik ent-~ i 


Mehr noch als die erste Mitteilung mu die folgende als ein Versuch Baselines ‘werden, Die 


lediglich aus ihrer Verwandtschaft mit dem altbewahrten Ritzschen Verfahren geschlossen, aber “a 


deren numerische Brauchbarkeit erfahrungsgemaB eine ausgezeichnete ist. Aber in dem vor-~ | 


uss Der strenge Mathematiker hat also Grund zu kritischer Einstellung. Thm sei gesagt, daB al 
_ die folgenden Ausfiihrungen keineswegs im Sinne eines mathematischen Existenzbeweises ver- i 
- standen werden wollen, sondern den Zweck haben, dem praktischen Rechner in allgemeiner a 
Form einen Naherungsweg zu zeigen, auf dem er in jedem Einzelfall zu einer mehr oder weniger 4 


_Naherung nachgeprift werden, und dies in méglichst vielen verschiedenen F allen tatsachlich — 


\ liche Abrundung erfahren, muB. Allerdings wurde jede einigermaBen grindliche numerische . 
Durchforschung der aufgestellten Lésungen den Rahmen dieser Arbeit ungebihrlich erweitern, — 


| 


; 
1 


me 


ek ein da oe 


waits 
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I. Absehnitt. Der Lésungsansatz. 


2. Ausgangsgleichungen und Kinfiihrung eines Naherungsansatzes. 
Ergebnisse aus der ersten Mitteilung tbernehmen; 
dort. Das in Rede stehende Stabilitatsproblem wi 


So eEe wir im folgenden 
verwenden wir dieselben Bezeichnungen wie 
ird formuliert durch die Variationsgleichung 


of Lar=0. (1) 


Darin ist nach Gleichung (12), (13), (19), (20) und (22) der ersten Mitteilung 
L = T@) Aigass Avonis Me 


G : du; , du)? 4 du; du 
Au = = ff{f > ee Se a 0 x; Far | dm dm dx, 
3 
oe ou, du (2) 
=) ( Tr r 
Ams. = 5 ff > k;) hae dx, dx, dx, , : 


TQ = Ff ed CRY ax ax das. 


Die ki, die Spannungen des Grundbewegungszustandes, sind als gegeben anzuschen. Gesucht 
sind ‘aie Zusatzverschiebungen u;, welche die Nachbarbewegungen zur Grundbewegung uy”? 
darstellen und tiber die Stabilitat des Grundzustandes Auskunft geben sollen. Sie unterliegen 
gewissen geometrischen Randbedingungen, die durch die-den Gesamtverschiebungen | 


w= Ut uz (i=l, 2, 3) 
vorgeschriebenen Blddbediveanéen unmittelbar gegeben sind. : 


ig Uber die AuBeren Lasten X; und pi machen wir in dieser Arbeit im n Gegensatz zur ersten 
Bes Mitteilung eine bestimmte Annahme. Die Lasten sollen als Summen eines zeitlich konstanten 
und eines zeitlich harmonisch veranderlichen Anteils gegeben sein: 


Xi=AX} Pex; cos wt 


=Ap; + ep;' cos wt 


ee x, x! Pi ‘und pi seien ae Ortsfunktionen. Berk Multiplikation mit den freibleiben- 
den Parametern / bzw. ¢ erhalt man daraus die in (3) stehenden statischen baw. veranderlichen 
Lastanteile in jeder gewiinschten GréBe. ele 


Gag 5: ve “Oe (3). 


© Ausden Lasten (3) mu8. der zugehérige Grundbewegungszustand berechnet werden, und 
zwar muB dies, da wir uns wie in der ersten Mitteilung auf ,,linearisierte Grundzustande“ be-. 
@ schranken: wollen, auf dem Boden der klassischen Elastizitatstheorie geschehen. Da in dieser 
Theorie das Superpositionsgesetz gilt und mehrdeutige Lésungen nicht vorkommen, so erhalt 
man die. Spereotecs des_ eee eindeutig in der Gestalt . 


~ 


— eS bee nee KD = 3 oy + 60% cos wt (j=); 2, 3) a oie i (4) 
“Die seanephos Spannungen hoy sind: darin ausschlieBlich durch i statischen Lastanteile be- 
‘stimmt,. die pulsierenden Spannungen £0% cos wt sind die Spannungen der von den pulsierenden 
Lastanteilen erzwungenen Schwingung, die von allen Rigenschwingungen befreit sein soll. 

| ' Die ojjund Oy sind somit Ae aousbone Ortsfunktionen zu betrachten. Die Oi j hangen” auBerdem 
von der erregenden Frequenz @ ab. Fir gewisse Werte von-w werden die oj; unendlich grof. 
Das tritt dann ein, wenn bei Zusammentreffen von @ mit einer Eigenfrequenz des Kérpers 
Resonanz fiir die Grundbewegung entsteht. Solche Werte von w schlieBen sich von selbst dadurch 
aus, -daB fiir sie die folgenden Rechnungen, in denen bestimmte endliche Betrage der oj; voraus- _ 
_geset nt werden, gar nicht. durchfihrbar sind. Wichtig ist, daB durch diesen Umstand im all- 


n eee samtliche i ae ear ae fir w verboten werden. Beispiele sind leicht 


nid 


eit 


ORG wn ae Uig Wig nines ae DN Win (i 1,8 3) “ wae - 6 
Neve peed : . 28 
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Darin:sollen die wij =w,(t) unbekannte Funktionen der Zeit und die ug =wi (x, %, x3) bestimmte 
fest gewihlte Funktionen des Ortes, sogenannte ,.Koordinatenfunktionen*’ sein, welche die den 
u; vorgeschriebenen geometrischen Randbedingungen erfiillen. Wie man sieht, stellt (5) eine 
einfache Prdclierine des von dem bekannten Ritzschen Verfahren her gewohnten Ansatzes} 
dar. Nach Ritz setzt man die Lésung eines Variationsproblems als Summe von bestimmten 
Funktionen an, deren jede mit einem freibleibenden konstanten Koeffizienten multipliziert ist. 
Wenn wir hier statt konstanter Koeffizienten unbekannte Zeitfunktionen stehen lassen, fiir die 
wir dann spater ein System von — den Ritzschen linearen algebraischen Gleichungen entsprechen- 
den — gewohnlichen Differentialgleichungen ableiten, so schlagen wir damit einen Mittelweg 
ein zwischen dem Ritzschen Verfahren als konsequenter direkter Methode einerseits und der 
Ableitung der Eulerschen partiellen Differentialgleichungen als klassischer Methode der Va- 
riationsrechnung andererseits?. Von gréBter Wichtigkeit fir die Gite der mit dem Ansatz (5) 
erzielbaren Ergebnisse ist zweifellos auch hier wie bei allen direkten Verfahren der Variations- 
rechnung die geschickte, dem jeweils vorliegenden Problem angepaBte Wahl der Koordinaten- 
funktionen. ; 


Um die folgenden Formeln tbersichtlicher zu gestalten, wollen wir den im Anschlu8 an das | 
Ritzsche Verfahren gebildeten: Ansatz (5) in eine auBerlich etwas andere Form umschreiben. | 
' Wir setzen an Stelle von (5) mit anderer Bezifferung der ui 


Uy = Uy, Wy, + Uyg Wg ee! + Uy, N—2 WN-2 
us = Uo, We + Ups Ws + +++ + Us, N-1 Wy_y (6) 
ug = Uz3 Ws + Ugg We + +++ + UsyWy 


-..worin nun die w;(t) die unbekannten Zeitfunktionen darstellen. Es ist klar, daB durch die Er- 
_ setzung von (5) durch (6) die Allgemeinheit des Ansatzes nicht eingeschrankt wird. 


_ Geht man mit (6) in (2) hinein, so kann man die raumlichen Integrationen ausfihven, und © 
: erhalt, wenn man eine zeitliche Ableitung durch einen Punkt bezeichnet, 


dtr tet eae Sa 


ne 


“_ 


aot te 5 Afm Wf Wm — ¥ Dim We Wm —E COs wt y Efm tay Wm 5 : 
“ fiesta reat | f,m=1 ' fym=1 
-worln Dfm die Form 


— 
Pe see! 


Dias = Bym +1 Chm 


hat. Die Koeffizienten Asm bis Eym lassen sich auf Grund des niaeraes (6) in der Form schreiben: 


ee olf 0 > Uif Uim Ax, dx, des, 


ba jup Huy Jum , du; 4y dug duj 
Br veer. ol, (( J. Jm if Im 
f > Sl Pal age cera ars ie 0x; )+ l—2y Ox; dx da, doty day, 


> 


3 3 a 
0u,f 0 Urm 
Che Coe 5) oil ae ae aay 3 dx, x5 dxs, . 
ts a 1 3 u : 3 sae 9 Urn a 
ees I} 28 > ie anew 2 eee ate é 


In Riesen ee dradken ist fir diejenigen Funktionen ux, die in (6) nicht vorkommen, Null ein- : 
; ausetzen. Es gilt offenbar_ a rah ee aris. 


ay 
aoe 


Ve Afra = Ans Bim = = Brp Chm = = Cop, Dim = Dry, Ejm = En 


‘Die ik oaffizionten Ae Be und Cfm sind auf. eee unserer. Festlegung der Koordinaten- 
funktionen konstante Leman die Eos F unktionen von Olas 


t 
. h 


ai Vel. z Z. oB. CB Biting und R, ‘Grananil: Technische bya "Bistin 1939, Se 136. h A 


. 2 Unser Verfahren ist grundsatzlich beka Es ‘AL. Poritsk: A; | 
“bridge (Maas) pooner Me mn ee ry, Ver a. (on Ant. ae i . Techn. Mech. Cam- i 
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3. Transformation des Problems in die Normalform. Die Form (7) von L eignet sich nicht 
gut zur Ermittlung der unbekannten Funktionen w;. Wir miissen L erst auf eine einfachere 
Normalform bringen. Wir suchen dazu eine lineare reelle Transformation 


W,' = Aj, Vy + Ayn Vo + +++ + 4:n UN (t=1, 2,...N) 
auf, welche die beiden quadratischen Formen 


(8) 


N N 
A= Da Ajm Wf Wm und D= De Dim Wf Wm 


f,m=1 f,m=1 
gleichzeitig auf die Gestalt 


N N 
A= > und D= > apy 
f=l f= 


bringt!. Eine solche Transformation (8) existiert, weil die quadratische Form A vermoge ihrer 


mechanischen Bedeutung als kinetische Energie positiv definit ist. Die w; genigen der Deter- 
minantengleichung 


| Dm —co; Ajm| = 0. (9) 


Wie man die Koeffizienten vo 


(8) praktisch auffindet, soll hier als bekannt vorausgesetzt 
werden. 


Durch die Transformation (8) geht die Funktion (7) in die Normalform 


N N N f 
' i > vy — >) wF F— e cos wt » Fm Uf Um _ (10) 
fal fal fm=l 


mit 


oe : 
Pym a yi) Ey aif Ajm, Pim = F'ng 


i,j=1 


iber. Darin hangen die w7 von A und die Fym von A und @ ab. Aus (10) folgen die Euler-Lagrange- 


; é d (aL aL he 
schen Gleichungen oF ( =) rei 0 des Variationsproblems zu 
N 
dp + wf vf + e cos wt >) Fim Um =0 (f= 1525...) 5 (11) 


_ deren Integrationstheorie uns in den nachsten Abschnitten beschaftigen wird. 


Zuvor ist es niitzlich, sich die mechanische Bedeutung der Transformation (8) klarzumachen. 


 Setzen wir in (10) ¢=0, so steht statt (11) einfach 


vet opup=0 (f=1,2,...N). 


e Die Lésungen dieser Differentialgleichungen sind (fur o> 0) harmonische Schwingungen 


‘(sin aft 
ee sin f (PHT 25. S5N):. 
cos wt 


_ Es handelt sich offenbar um die Eigenschwingungen des zeitlich konstant vorbelasteten Kérpers, 
_ mit anderen Worten um die Nachbarbewegungen zur statischen Grundgleichgewichtslage. Die 

GréBen ws sind die Eigenfrequenzen. Die geometrischen Eigenschwingungsformen findet man 
sofort durch Einsetzen von (8) in (6): Es entsteht 


yo: N : peel, N i 
ui vy > Uj, Ary oir V2 > Mir Gyo + ae a UN a Uj, A, (i=l, 2, 3). 
r=1, aah oer r= 


Die k-te Eigenschwingung hat also die Schwingungsform 


Pec TENE ) 
uk = >i ura (i =I, 2, 3). 
rea i 


ee. : 1 Vgl. zu dieser ,,Transformation auf Haupt- oder Normalkoordinaten“ E. T. Whittaker, Analytische 


Dynamik, Berlin 1924, § 77—78, oder J. Horn, Gewohnliche Differentialgleichungen, 3. Aufl., Berlin und 
Leipzig 1937, § 28. Save : mle eye ny 


28* 
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cK 


Wir haben somit (jetzt wieder fir das allgemeine Problem e+0) vermoge (8) den Lésungs- 
ansatz (6) in 


tip wh oy th ay he ay te eB) (12) 


iibergeftihrt, d.h. wir setzen jetzt fir die Zusatzverschiebungen u; eine Summe aus Produkten 
der Eigenschwingungsformen uj, des statisch vorbelasteten Kérpers mit unbekannten Zeit- 
funktionen vy. Durch diesen besonderen Ansatz wird die Funktion L auf die Normalform (10) 
gebracht. 

Die Eigenfrequenzen @; stehen in einem wichtigen Zusammenhang mit den statischen Knick- 
lasten des elastischen’ Kérpers. Setzen wir in der Gleichung (9), die ja identisch ist mit 


zunichst A=0, so sind die N Wurzeln w; als Frequenzquadrate des unbelasteten Kérpers alle 
positiv, die @; also alle reell und +0. Daran andert sich aus Stetigkeitsgriinden auch nichts, 
wenn / kleine Betrage annimmt. Die statische Gleichgewichtslage ist also fir geniigend kleine 
Belastungen stabil. Wachst |A| weiter an, bis A schlieBlich die erste (positive oder negative) 
Wurzel A, der Gleichung 

| Bow A Gpo| = 0 (14) 


erreicht, so wird mindestens eine Wurzel w; von (13) Null, wie man durch Vergleich von (13) 
und (14) sofort sieht. Die fur A=A, nicht verschwindenden Wurzeln sind positiv geblieben. 
Offenbar ist nun auf den Kérper gerade die erste statische Knicklast aufgebracht. Um dies 
einzusehen, braucht man nur mit dem Ritzansatz (6), jetzt mit konstanten Koeffizienten wx, 
in die statische Gleichung 6(Ajiass + Aus) = 9 (mit ¢=0) hineinzugehen; man kommt: dann 
sofort auf die Gleichung (14) fir die Knicklasten. 


Was weiter geschieht, wenn A die Wurzel A, iberschreitet, soll hier nicht verfolgt werden. 
Vielmehr lassen wir fiir A weiterhin nur solche Werte zu, fiir die der Kérper im stabilen Gleich- 
gewicht ist; denn es erscheint uberflissig, sich mit der Stabilitat dex erzwangenen Schwingungen 
eines elastischen Kérpers zu befassen, der schon statisch iiber die Stabilitatsgrenze hinaus be- 
lastet ist. Wir halten Jediglich als Ergebnis der Uberlegung fest, daB mindestens eine der samtlich 
reellen Eigenfrequenzen @; auf Null absinkt, wenn sich die statische Last von Null her einer 
Knicklast nahert. 


_ Die Gré®e der statischen Last wird tbrigens unter Umstanden nicht nur durch das Erreichen 
der statischen Stabilitatsgrenze beschrankt, sondern auch durch die Bedingung des ,,lineari- 
sierten Grundzustandes‘‘, namlich der Forderung, daB sich die Beitrage des statischen und des 
schwingenden Lastanteils zum Grundspannungszustand ohne gegenseitige Beeinflussung ein- 
fach tberlagern missen. Was das heiBt, macht man sich am besten an einem Beispiel klar: 


An dem in Abschnitt 5 der ersten, Mitteilung behandelten Stab stéren sich offenbar der 


_statische und der schwingende Grundbiegungszustand in y-Richtung gegenseitig so lange nicht, | 


als es sich um eine reine Querbelastung des Stabes in y-Richtung handelt. Die statische Last 
_darf bis zur Kippgrenze ansteigen, ohne dafs sie die erzwungene Schwingung beeinflu8t; ins- 
besondere verschieben sich dabei die Biegeeigenfrequenzen des Stabes in y-Richtung, die allein 
als Resonanzfrequenzen der erzwungenen Schwingung in Frage kommen, nicht. Davon un- 
_abhangig andern sich die Biegeeigenfrequenzen in z-Richtung und die Torsionseigenfrequenzen 


mit A in der dargelegten Weise. Anders ist es, wenn der Stab auch noch durch eine statische | 


Langskraft (in x-Richtung) beansprucht wird. Dann sinken namlich auch die Biegeeigenfre- 
quenzen in y-Richtung in bekannter Wcise ab, wenn sich die Langskraft der Knicklast nahert, 


und damit werden die Resonanzfrequenzen der erzwungenen Schwingung verschoben. Die 
klassische Theorie gibt diesen Effekt allerdings nicht wieder. Damit sie aber fir die Berechnung 


der erzwungenen Schwingung genau genug bleibt, muB die statische Langskraft auf Werte weit 
unterhalb der Knicklast begrenzt werden. — 


. 


4, Kin Beispiel. Zur Veranschaulichung der bisherigen allgemeinen Austahrun een sei hier 
ein konkretes Beispiel behandelt, das schon in Ziff.5 der ersten Mitteilung betrachtet wurde: 
der flache Rechteckstreifen in Gabellagerung mit Belastung durch Endmomente 


My =AM, +My, cos wt ya uy 
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(Abb. 1, Bezeichnungen wie in der ersten 
Mitteilung!). Wie frither abgeleitet, kann 


in vertikaler Richtung als unwesentlich 
auBer Betracht lassen (vgl. dazu die Be- 
merkung am Schluf dieser Ziff.) und hat 
somit nur die Querverschiebung w der 
7 Stabachse in waagerechter Richtung so- 
Maley ete Bear vehdndchie’ steven. wie dan\Vecdsohvantel @ des Querschnitts 


in den Stabilitatsgleichungen zu _beriick- 
eetineh: Wir haben also nach den Gleichungen (30), (31) und (32) der ersten Mitteilung 


I 


l 
done = Boe (Beans 22 (Dan as 
0 0 
R 1 5 
Cae eae eae (Mijdx, (16) 
ry 5 B Alles \ ey Jax. | | (17) 


_ 


Die Grundbew egung des Stabes ist die von den auBeren Momenten M,, erzwungene Biegeschwin- 


gung in lotrechter (v-) Richtung. Sie hat nach den allgemeimen Aussagen der Stabkinetik? das 
Biegemoment 


ES ; 3 ka x 

oY Eo etl : cO y . 
ee MAM, (Oe cos is ae Ee (18) 
. anes 

oe = ! PY 

- _wobei der Strich am Summenzeichen besagt, da® nur ungerade ‘A auftreten, und 

oa 5; Sie i Ra \Vae re J, 

a ae ; ate oe ( 1 ) V7 be ety 
“die ee eiccsien des Stabes in y-Richtung bedeuten. Gleichung (18). entsp-icht der 
_allgemeinen Gleichung (4). 


Wenn wir nun den Naherungsansatz (6) aufstellen, so ist. zu beachten, daB die Cetiecehilen 
- Ausdriicke (15), (16), (17) nicht mehr wie (2) aus den drei allgemeinen, von drei Koordinaten 
cabbaagier Zusatzverschiebungen gebildet, sondern soweit vereinfacht sind, daB sie nur noch 

die zwei von x allein abhangenden Funktioen w und @ enthalten. Das Problem ist also im spe- 
~ ziellen Fall unseres Beispiels zweidimensional. Wir setzen deshalb an die Stelle des allgemeinen 
- dreifachen Ansatzes (6) einen entsprechenden zweifachen fiir w und #, in dem die Koordinaten- 


- funktionen Funktionen von x sind, welche die ee fir x—0 und x =I (Gabel- 
lagerung) erfiillen: 


w= me ae sin Aes = tus sin - 


a een Seb + wy yasin ae | can 


aa i se ec = msn * 4 wgsin = Bri > ++ Wey Sin 


Oe uns Be ay DF kA es kat dx ae ae (7) ye | wth foo =e ae 


/ 


Ae AEM Dy wit St el (F me ee cere te (20) 


(18) vorbehalten bleiben mu8. Die Faktoren 2 und é sind in den Ausdruck fiir das ine cia Mr eingefiigt, 
ih Einklang mit.(3) zu_ bringen. ve 
ee Val etwe I henenser ae We: Prager okapais der Stabwerke, Kap. IIB. Berlin 1933. 


man die Verschiebung v der Stabachse' 


a Wir schreiben jeu. silerdities Mr statt des dort verwendeten Puplapsbens M. der fiir das Biegemoment 


a a 
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Ebenso kommt aus (17) i i 
T(2) = pl ye ok Hla | 105 L (21) 


Die Arbeit (16) berechnet man am bequemsten, wenn man sie in der Form 


Mies ee 
0 


schreibt. Man erhalt dann aoory 


N 
Awa = AM, (2) 5 Se | Wee + € cos wt Y" Ey i,24—1 W2i W2k-1> (22) 
4 i,k=1 ‘ 
worin entsprechend der Schreibweise (7) abkiirzend 
4 2 | kax & som os 
3 Sane hice aoe , 
Eo ign = My — (+) k? if sin —— sin —>— SS; : dx (23) 


\ . i 7 Ll = . (23) 2 
0 ll Cjy 


gesetzt ist (wahrend wir die leicht ersichtlichen ubrigen Konstanten von (7), namlich Ajm; Brym 
und Cyn nicht besonders aufgeschrieben haben). ; ; 


Beachten wir nun bei der Auswertung des pee in (23) die elementaren Formeln 


sin a sin B = 5 cos (a —f) — + cos ( (a +f), ee 4 


sin « cos B =~ sin Cee obs > sin (a — B); 
pate * Neh F 2 
ee hace a : : ‘far gerade k, : | ed 
l eo far ungerade k, . 

kn 


0 


\ 


so folgt sofort ; 

(sy eae eae ctr gies a falls i t und k beide gerade oder ungerade, : 
. Eoiena | 

: =0, fallei Oe k gemischt. 

In Worten ausgedriickt: die Gesamtheit dex Efm zerfallt in unserem: Beispiel in zwei verschiedene 
Gruppen. Die eine Gruppe ist ausschlieBlich gebildet aus den ,,geraden‘‘ Koordinatenfunktionen 
sin kax/l mit geradem k, die andere Gruppe aus den ,,ungeraden‘* Koordinatenfunktionen mit 
- ungeradem k. Eine »»Mischung*‘ der verschiedenen Koordinatenfunktionen tritt nicht ein. Dieselbe 
Feststellung kann man auch fir die ibrigen Arbeitsausdriicke (20), (21) und den mit 4 multi- — | 
plizierten Teil von (22) (die Koeffizienten Afm, Bfm und Cfm) treffen. Daraus folgt weiter, daB 
= (a aL ee OL. . 

0 WK 0 Wr 
von Differentialgleichungen kommen wird, von denen man jeweils das eine oder andere allein 
_erhalt, wenn man den Ansatz (19) ausschlieBlich mit ,,geraden” oder ,»ungeraden* Koordinaten- 
_funktionen einfuhrt. 


Wir wollen in Fortfihrung des Beispiels_ nur noch die Gleichungen betrachten, die aus den : i 


-_tngeraden Koordinatenfunktionen entstehen, . da die ibrigen dann nichts Neues mehr bringen. 
Wir schreiben daher, um die wesentlichen Formeln kurz zu wiederholen, mit anderer Boa ane 


man auch bei der Bildung der Gleichungen —— =0 auf zwei vollig getrennte Systeme Ae 


ales M1, statt (19) Be > d So 


+ ates + Wn] sin in =D ne ; Ae 


3 x 
w = sy sin “7 + mn sin 


ee a sae O <u sin 9% 4 1 sin oe ty in CDE “i 


4 


‘ Har and. folgt, wenn wir den mit. A multiplizierten Teil \ von Aas mit Acne abkarzen, nach 20, 
= (21) und (22) oe = 


Jy ak ali A Naat eas ea 0s ee mR Cae a 

An + Aba = Fp 5) ES (GE) a+ 2AM, + Cena]. 5) 
ie MWR erred Git at Me aces ar. MADR a ay A Dak hae 

i 9). Bh pe BU Te eSATA te ; § 4 

Te) eA Wi, + on Whar Bek el cc) 

k=1- Bure k=1 x « i 


mn = 
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Auf den mit ‘e cos wt multiplizierten Teil von Aj; kommen wir sogleich zuriick und erledigen 
zunachst die Aufgabe der Einfihrung von Normalkoordinaten nach Ziff. 3. ae 

Der Ausdruck (25) besteht aus einer Summe von n/2 quadratischen Formen mit je zwei 
Variabeln w,, w;,,;- Die Hauptachsentransformation dieser einzelnen Summanden laBt sich leicht 
ausfiihren. Man erhalt, wie aus der analytischen Geometrie bekannt oder durch Rickwiarts- 
einsetzen leicht zu bestatigen ist, die den Gleichungen (8) entsprechenden Transformations- 


formeln 
we = note (v, COS %, — Vz, SIN O,), z 
(ey say ee, cos!) 
ip Yul (eal, Sean 1): (27) 
tp ay 


wahrend (9) in die Gleichungen 


| AN no} =r) 
| EJ,(=) ee AM, (=) 0 eh 35. on I) (28) 


AMT) 6s(“P) —wipo} 


tibergeht. Die Wurzeln von (28) sind @,; = @,%, Wy... ° Mit (27) und (28) geht (26) in 
T) = of +--+ +02 und (25) in Ayes; + AAms = 02? 0? + w2 v2 + --- +2 v2? ber. 

Die Uberlegungen, die an (14) anknipfen, gestalten sich folgendermaBen: Fur k=1, wo; =0 
liefert (28) sofort 
by = aT VEC I Sa 
als kleinsten A- Wert, fiir den statische Instabilitat (Kippen) eintritt. Nun ist AM, das tatsachlich 
aufgebrachte statische Moment, also Bed oe 

May =I, 2 VEC, 
das (bekannte!) kleinste Kippmoment. . 
SchlieBlich haben wir noch den mit € cos wt multiplizierten Teil von A,,, zu betrachten und — 


bezeichnen ihn zur Abkirzung mit A,,,. Es ist jetzt 


= higgneren : im x Ep £2 ae 
wie inky A ee eptecy 4 2 : . red eed ses . eel eR v 90 
Aims = & M, cos wt— (7) ee k? w;,. W, i: sin —— sin — 2; | ae dx. (29) 
inn | 0 jJ=1 4] Les (--) ‘ 
: Djy. 
Die Integrale kénnen ausgewertet werden, worauf wir noch nicht hier, sondern erst in dem 
Zahlenbeispiel in Ziff. 12 eingehen. Fihrt man dann (27) in (29) ein, so erscheint Aims in trans- 


- formierter Gestalt mit den Variabeln v, und kann mit Ajjas +} AAms und T() zu der Normal- 
form (10) von L zusammengesetzt werden, aus der die Differentialgleichungen (11) folgen. 
Ein wichtiger Sonderfall tritt ein, wenn w so klein ist, da8 man 1 —(w/ajy)? =1 fir alle jens 
- setzen kann. Dann ist namlich auf Grund der bekannten Fourierentwicklung 
| OR RON Se Sa 5 
== 8 ae far 0< x <l (30) 


j=! Z 
die Summe 2%” in (18) ebenfalls gleich 7/4, und aus (29) wird einfach 

a n—1 

5 — Efe 

3 Ars aio i eee cos wt 2 k? wy Wy 15 
‘worin noch die Transformationsformeln (27) einzusetzen sind. Man erkennt jetzt, daB die 
_ Differentialgleichungen (11) in die getrennten Gleichungspaare 

Dy, + wpry + cos wt (Fav, + Fi, aires) = 9, 
Meir + Whi Mex + COS wt (Fria, 4 Ue + Frypans eri) = 9 


1 Auerbach-Hort, Handb. d. phys. u. techn. Mech., Bd, IV, 1, $.122. Leipzig 1931. 
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zerfallen. Das ist der Fall, den Woinowsky-Krieger! in etwas anderem Zusammenhang und mit 
anderen Hilfsmitteln direkt behandelt hat. Wir wollen hier auf diese Dinge nicht naher eingehen, 
sondern zum Schlu8 aus den vorstehenden Uberlegungen lediglich die Tatsache festhalten, 

daB das in Ziff. 3 allgemein abgeleitete Differentialgleichungssystem (11) in konkreten Einzel- 
fallen unter Umstanden in mehrere getrennte Systeme zerfallt, also hier beispielsweise in die 
beiden zu den geraden und den ungeraden Koordinatenfunktionen gehorigen Systeme und dann 
unter bestimmten Voraussetzungen sogar noch weiter in die genannten Gleichungspaare. Es _ 
gibt sogar Falle (der axial pulsierend belastete Stab gehért dazu), in denen sich (11) in lauter _ 
unabhangige Einzelgleichungen, also in Mathieu-Gleichungen, zerlegt, was aber hier nur erwahnt 

sei. Wenn wir nun in den nachsten Abschnitten die Integrationstheorie der Gleichungen (11) 
behandeln, so haben wir dabei, falls der erwahnte Zerfall eintritt, stets eines der selbstandigen” 
Teilsysteme vor Augen. Die Zah] N, welche den Ansatz (6) charakterisiert und die Zahl der Glei- 4 
chungen und Variablen von (11) angibt, kann also gréBer sein als die Zahl m der Gleichungen 
des Systems, fiir das wir in den nachsten Abschnitten die mathematische Theorie entwickeln. 
Deshalb haben wir schon von (24) ab n statt N geschrieben und rechnen von jetzt ab allgemein 
mit der Zahl n = N weiter. 


Noch eine Bemerkung sei der Vollstandigkeit halber hier angefigt: Um ganz korrekt zu sein, 
miiBte man den Ansatz (19) baw. (24) durch einen entsprechenden Ansatz mit sinus- Koordinaten- 
funktionen fiir die lotrechte Zusatzverschiebung v der Stabachse vervollstandigen. Dieser Ansatz 
fiihrt, in die in (15) und (17) weggelassenen Glieder 


Apis = a ex Jax und T(?) = “(8% ae 
0 0 


(vgl. (30) und (31) der ersten Mitteilung) eingesetzt, zu weiteren linearen Differentialgleichungen, 
welche von den anderen unabhangig sind und die freien Eigenschwingungen des Stabes in v- 
Richtung beschreiben. Den periodischen Faktor cos wt besitzen diese Differentialgleichungen 
nicht, weil der Ausdruck. A,,, die Verschiebung v nicht enthalt, so daB also das pulsierende 
Biegemoment M gar nicht eingeht, oder weil anders ausgedriickt die Arbeit der Zusatzver- 
schiebung v an den Grundspannungen keinen quadratischen Anteil A, liefert. Deshalb kann 
die Zusatzbewegung in v-Richtung keine Instabilitat zeigen und ist eben aus diesem Grunde 
uninteressant. 


II. Abschnitt. Allgemeine Betrachtungen. 


5. Theorie der Differentialgleichungen mit periodischen Koeffizienten. Wenn man in die 
Differentialgleichungen (11) die neue Variable 


s=ot (31) 
einfihrt und gleichzeitig n statt N schreibt, so kommt ; 


79 A iy n ~ 
wort Te + ofr, + cos s 2, Fim tm = 0 (fal 2,...m). (32). 

Das ist ein System von n, homogenen linearen Differentialgleichungen zweiter Ordnung mit 
- periodischen Koeffizienten. Seine Lésungen geben den zeitlichen Verlauf der Nachbarbewegungen 4 
zum Grundbewegungszustand wieder, lassen also unmittelbar das Stabilitatsverhalten der Grund- 

_ bewegung erkennen. Die allgemeine Theorie solcher Gleichungen ist an sich bekannt. Doch gelten 
fiir das vorgelegte System einige Besonderheiten, und wir miissen als Wegweiser fiir seine im’ 

Ill. Abschnitt folgende praktische Auflésung die Hauptpunkte der allgemeinen Theorie hier 

entwickeln?2, fe aa Sa i ae nea 


Wir betrachten ein System. von Partikularlésungen von (32): 
vf = of; (8), tyo(s) 2+ +. M% 2n(8) (FSL...m). 2 eee 
1 §, Woinowsky-Krieger, Ing.-Arch. 13 (1942), S. 197. 


? Vel. au den Ausfithrungen dieses Abschnitts J. Horn, Gewohnliche Ditieeeatial leichunge ane a 

ssi ; : 5 p i ungen, Sammlung ~ 

oe § 55, Leipzig 1905, oder H. Poincaré, Les méthodes nouvelles de la aetanigud ets oa § 29. ey) | 
‘ aris 92. ie: : be aaah ' : ae . a ea . | 


\ 
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ee ee 
_ Die vse und ihre Ableitungen oe sollen den Anfangsbedingungen 
A fork =f 
vx (0) = 
Oe ON S| 0 Mie bef 
SRLS aie bie des Ay Byes : 
dof, (0) fen if fir k=ftn (f fy 3 aH) (34) 
ds HO fit k+ftn 


 gentigen. Das System (33) bildet dann ein Haupt- oder Fundamentalsystem von Lésungen, die 
sich, mit 2n willkirlichen Konstanten multipliziert, zur allgemeinen Lésung zusammenfigen 
lassen. Der Begriff des Hauptsystems von Lésungen wird in der Differentialgleichungstheorie 
_ gewohnlich nur fiir Systeme erster Ordnung aufgestellt. Die Ubertragung auf das vorliegende 

System zweiter Ordnung bedarf aber wohl keiner naheren Erklarung. — és 


; Da sich (32) nicht andert, wenn man s durch s+2 7 ersetzt, so mussen die Losungen (33) 
~ Integrale von (32) bleiben, wenn man in (33) s+-2 a an die Stelle von s setzt. Diese neuen In- 
“tegrale miissen sich mit bestimmten Koeffizienten «;, linear aus dem Fundamentalsystem (33) 
_ gusammensetzen lassen: 


Uk (S+2 at) = a4 % Vpy (S) ee Upo (8), ++ ite n, Opa (8). (F225 miko t, 229 n). (35) 


- Setzt’ man darin s=0, so folgt aus (34) fir die Konstanten x, sofort " 2 ioe 


ae ape = vpk (200) 3 Ze Cee 
A (fale oy cat hea la 6h : ee 
Oftnk = a (2m) : : ee 


sf Wir suchen nun besondere partikulare Loésungen V;,(s) von (32), die an Stelle von (35) den 


einfachen Beziehungen . 


US Vaden) = 8. Ves) fo? my oe peer 


mit konstantem S, gentigen. Der Index r wird zunachst ohne nahere Bestimmung mitgefihrt. 
‘Auch die Vf miissen sich als lineare Kombinationen der Vek Mit gewissen Konstanten Bir dar-— 
stellen lassen: ; rt : 


Vie = Bay Uy Por Ufo a aoe + Ban,r Vp2n+ oes (38) | 
Schreiben wir mit diesem Ansatz die aus (37) folgenden Gleichungen Vp-(2 2). = S, Vjr(0) und | 


£Th (2 2) =S ate (0) hin, so erhalten wir unter Beachtung yon (34) und (36) das Gleichungs- 
ystem fiir die Koeffizienten Birs +s Bane ie 
Ne Pat Ai bie tia On ae 
ee ee Pas Oo ats Bar (%22 —'S,) sae t Ban,r 2,2 Teal FN a (SO 
oe : Siege : : Bir zs gee : + Baran, 2 +. . Soa Ban, r(%9n,2n— Sr) = 0 
dk die Bestimmungsgleichung fiir Sees Seep ean Me Siegal 
: Ge as sie Osha [oar tas Ve Olen : | : 


a 


tzen wir weiter 


4 


(42) 


ach (37) eine periodische Funktion von s mit der Periode 2.7 ist. Die . 


-charakteristische Exponenten’ genannt. — 


: Ver = & pp (s) “ (f=1,2,...m). 


wy 


) 
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Weitere Ergebnisse gewinnt man, wenn man in den Differentialgleichungen (32) und ihren 
bisher diskutierten Losungen die Veranderliche s durch —s ersetzt. Das System (32) andert sich 
durch diese Substitution nicht. Folglich missen alle Lésungen von (32) Integrale bleiben, wenn 
man in ihnen —s an die Stelle von s setzt. Was zundchst das Lésungssystem (33) betrifft, so 
erkennt man aus der eben angegebenen Tatsache sofort, daB jedes Teilsystem v);,..2, Un~ mit 
k <n lauter gerade Funktionen-von s enthalt, fiir die also vsx(s) =vsx(—s) gilt, weil sich weder |} 
die Differentialgleichungen noch die Anfangsbedingungen (34) andern, wenn man s durch —s 
ersetzt. Ebenso ergibt sich, daB fir k>n die Beziehung vs.(s) =—vsr(—s) gilt, daB diese vx 
also ungerade Funktionen von s sind, weil bei Ersetzung von s durch —s die Anfangswerte 
der Ableitungen das Vorzeichen andern. 

Betrachten wir ferner die Losungen Vj, welche die Beziehungen (37) erfiillen. Vr(s) hat” 
die Eigenschaft, daB és sich mit einer Konstanten S, multipliziert, wenn man mit s um 27 ,,vor-_ 
warts“ schreitet. Folglich multipliziert sich Vy-(s) mit 1/S,, wenn man mit s den Schritt 27 
rickwarts, d. h. in Richtung negativer s geht. Also erfillen die Funktionen 


Vi.(s)=Vel—s)  (f=1,2...-n), 


die nach dem oben Gesagten gleichfalls Lésungen der Differentialgleichungen (32) sind, die 
Beziehungen 


: 
Vit (s + 20) = 5 Vp (s) (f= hiZeo kn) (43) |] 
Es gibt also zu jeder Lésung (42) eine zugeordnete Lésung 


| Aan t7 =e rg (so (fal, 254. 0), (44) 
worin qj; (s) = gr(—s) periodisch mit der Periode 27 ist. ‘ 


Setzen wir auch Vj, in der Form (38) an, so folgt auf Grund der eben ausgesprochenen ices! | 
schaften der Funktionen vyx a 


‘ bef + 


oe Vir Barty + + Bar Ofe — Busiyr Unter — 9 ~ Bamr Mone PLAN HI 
Aus diesem Ausdruck ergibt sich auf demselben Wege, der von (38) zu (39) fihrte, ein Gleichungs- — 
system, von dem wir nur die ersten n Gleichungen brauchen: 5s 

Par (“4 as EBay Qin — Bazin Gini > “t+ — Pony Cen = 0, 1° 


; : 1 j : . 
: Bir Omi + $52 OB wr. (ann =<) — Bia, 7.0 aaa mee eg — B2nr On, 20 = 0. 


Addiert man die f-te Gleichung (45) zur f-ten Gleichung (39) fiir pack, 2,...n, so kommt et 
sofort auf die einfachere Gleichung fiir S, i . 


(45) | 


rae 
| sakeye ; | ia 
| Oy iP Oa Sevan ++ Ain 
‘| or’ Dogacierrrte « Son, l=0 . he ‘4 
pee % ‘ ; ea ’ rf 6) fs 
at aay : e ny CEE ha heen elie I SD 
worn ~ : ae 
vas ir i z=—| —— Smee. gn 
ee he al ay | eet 


KSPR DERN EL LA BLAND COLAO OLE SE 


os Besetat ist. ee ee corse Gleichung. Sie ist in der Ausdrucksweise depo Alwetice 
- eine reziproke Gleichung fir S,, die mit jeder Wurzel S, a die: 'S, -besitzt: Dasselbe 
Ue Groaihimiek eto arent e briegin rks bie 
_ In z, ist (46) vom n-ten Grade. Der Index r soll, wie wir jetzt nachtraglich festsetzen, die _ 
Nummer der Wurzel z, von (46) angeben, muB also von r =1 bis r=n laufen. Falls (46) n verseined 
dene Wurzeln 2, +1 hat, sind also 2 n verschiedene Werte S, vorhanden, fir die man aus (39) 
- 2n linear unabhangige Lésungssysteme der fir findet. Damit fuhrt (38) zu n Paaren von parti- 
e kularen Lésungen Vj, Vi (f, r=1, 2,...n) der Form (42), (44), die zusammen ein neues F unda-— 
_ mentalsystem bilden. Ist dagegen eine Wurzel z,=+1 oder —1 oder eine mehrfache Wurzel z, 
vorhanden, so liegt mindestens eine mehrfache Wurzel S, von (40) mit mindestens zweifacher 4 
a Vielfachheit vor. Dann liefert (39) nach den Satzen der linearen Algebra nur dann 2n linear. 
_ unabhangige Lésungssysteme der Bir und damit Integrale (42), (44), wenn die Determinante (40) 
MA, fiir jedes solche S, héchstens den Rang ,.n minus Vielfachheit der Wurzel S,“* hat. Ist der ae 


Sls Sia cated nny’? 


ae » 


7 ‘, 
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groBer als diese Zahl, so 1laBt sich das Fundamentalsystem nicht mehr vollstandig aus Funktionen 
der Form (42), (44) aufbauen. In diesem Fall leitet man, wie wir ohne Beweis erwahnen, in der 
allgemeinen Theorie der linearen Differentialgleichungen mit periodischen Koeffizienten andere 
Integrale ab, die das Fundamentalsystem vervollstandigen und die sikulare Glieder, namlich 
_ gewisse Potenzen von s auBerhalb der periodischen und der e-Funktionen, enthalten. Wir brauchen 
, auf die genaue Form dieser entarteten Lésungen nicht im einzelnen einzugehen, sondern merken 
_ nur an, daf} vermége der sakularen Glieder stets gewisse Teile der Lésungen mit s dem Betrage 
nach unbegrenzt anwachsen, was ja fiir die Beantwortung der Stabilitatsfrage allein wesentlich ist. 


a & Formulierans eines allgemeinen Stabilitatskriteriums. Die Integrale Vj; und V#, welche 
die in Ziff.5 entwickelte traditionelle Theorie der linearen Differentialgleichungen mit periodi- 
schen Koeffizienten in der Form (42), (44) geliefert hat, sind vielfach nicht reell. Das hangt damit ie, 
zusammen, dah (40) und (46) keine Sakulargleichungen sind, sondern komplexe Wurzeln be- me 
_ sitzen koénnen, was auf zweierlei Arten denkbar ist: Entweder liegt ein konjugiert komplexes 
sf _Wurzelpaar z,, %, von (46) vor. Dann ist auch S, und S; nach (47) konjugiert komplex. Oder 
_ %, ist zwar reell, aber die zugehérigen Zahlen S, und 1/S, sind konjugiert komplex. Auf Grund 
a2 algebraischen Beziehungen (38), (39) werden im erstgenannten Fall die Lésungen Vj, und 
Vy konjugiert komplex ausfallen, im zweiten Fall die Lacunees Vj und V7. Aber selbst wenn 
Vj, (und damit auch V;;) als Ganzes reell wird, bleiben die Faktoren e®,* und Py, (s) ima allgemeinen 
doch komplex, denn “ie Gleichung (41), die wegen der Periodizitat der e-Funktion nicht ein- 
deutig nach 0, aufzuldésen ist, hat stets komplexe Wurzeln 0, mit verschiedenen Imaginarteilen. 
Im III. Abschnitt wird es: ‘sich zeigen, daB man bei der praktischen Berechnung der 0, zweck- 
-mafig solche komplexen Werte zugrunde legt. Ohne jedoch auf die komplexe Form der Lésungen 
-und die Frage ihrer Verwandlung in eine reelle Form im einzelnen einzugehen, kénnen wir die 
-wesentlichsten Aussagen tiber unser Hauptproblem, namlich die Stabilitatsfrage, schon auf 
Grund der - jetzt vorliegenden Ergebnisse machen. Wir kénnen diese Aussagen wahlweise auf 
awei KenngréBen stiitzen, entweder auf die Realtcile der charakteristischen Exponenten @Q, 
' oder auf die Wurzeln z, der charakteristischen Gleichung (46). Die Exponenten 0,, die tiblichen 
_ Hilfsmittel der einschlagigen Stabilitatsuntersuchungen, werden wir im nachsten Abschnitt aus- 
-fihrlich diskutieren. Zunachst ist es aber einfacher, das abschlieBende Ergebnis des vorliegenden 
3 Abschnitts: cunmittelbar an die Wurzeln z, baw. S, anzuknipfen. 


‘Die Zahlen S, geben vermége (37) und (43) Auskunft tiber das Verhalten dér Integrale (42), 
. (44) far wachsendes s (d. h. nach (31) fiir wachsende Zeit), also tber die Stabilitat der Grund- 
bewegung. Sobald. eine Zahl S,, sei sie reell oder komplex, einen Betrag +1 hat, wachst nach 

a 43) entweder V,,| oder | Viz | unbegrenzt mit s an, so daB mechanisch also Instabilitat 
Zwi ischen z und S,° bestehen nach ss ed ig ag 


; i ‘ = ae i Sat ea Ss. } t : , i 3 4 
CMGI Gata ate mag erase peg FAS ine western! (48) 


ennt aus haar ices dab fiir ie, me Boa A (<1 pve Beeide! |S, | =Hys afl haleed: tea 
r reelle %, mit |z[>1 und fii komplexe z, ein Wert |S,|+1 entsteht. Daraus und aus 
n uber dig: ee Waurzeln eee ergibt sich flcendee ma sa Stabilitats-— Ek 


che Gleichung (46) iy veluchiedeile reelle Warsi z. mit (ee sO 
abil. Ist eine Wurzel z, von (46) entweder komplex_ oder reell mit 
veg n: Hat Nees eine Wurzel A mole oder ab ne eine Pei oe 
nderfall, terr 
Rang ae iner od 
treffende. PEntcheding wie sich zeigen wird, pean von. 


lenmaBig ge, mee ans ‘$0 iegt 

erschiedene Raa der Og ieutoseed (11) bzw. 
eines der bekannten. numerischen Integrationsverfahren 
toh: immer emai: die Popmmuon (33) und 
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damit die Koeffizienten (36) der charakteristischen Gleichung (46) [sowie, falls notwendig, der 
erginzenden Gleichung (40)] aufzusuchen, so daB auf Grund des eben zusammengefafiten Kri- | 
teriums die Stabilitatsfrage eindeutig entschieden werden kann. Ferner liefern die Wurzeln |f} 
von (46) tiber (47) und (41) die charakteristischen Exponenten o,, deren Realteile als ein Mab | 
fiir die ,,Starke‘’ des Anwachsens der Lésungen (42), (44) betrachtetswerden kénnen. Die: @, jf 
geben gewissermafen den Grad der Instabilitat einer elastischen Bewegung an und werden | 
insbesondere dann von Wichtigkeit sein, wenn es sich um die Beurteilung der praktischen Be- | 
deutung der Instabilitat im konkreten Finzelfall handelt. 


III. Abschnitt. Praktische Lésung des Stabilitatsproblems durch Reihenentwicklung. 


7. Ansitze fiir die Reihenentwicklungen. Dic in dem vorhergehenden Abschnitt gegebene 
Lésungsmethode des Stahilitatsproblems ist zwai allgemein anwendbar, hat aber die unver- | 
kennbare Schwache, keinen direkten Einblick in die Abhangigkeit der Stabilitatsverhaltnisse | 
von den physikalischen Daten zu liefern, Die im Schrifttum vorhandenen Untersuchungen! | 
sind in dieser Beziehung weiter vorgedrungen. Man findet bei Durchsicht der vorliegenden 
Einzeluntersuchungen immer wieder als wesentlichstes Ergebnis ein Diagramm, in dem die 
beiden hauptsachlich wichtigen Parameter, namlich die Lastamplitude und die erregende Fre- | 
quenz oder gewisse aus ihnen gebildete andere GréBen, auf den beiden Achsen eines ebenen |} 
Koordinatensystems aufgetragen sind, wahrend die Koordinatenebene in bestimmte Stabilitats- | 
und Instabilitatsbereiche aufgeteilt erscheint. Ein solches Stabilitatsdiagramm ist mehr als eine |} 
mathematisch korrekte, nimlich eine auch physikalisch befriedigende Lésung. ¥ 

Wir setzen uns nun fiir unsern allgemeinen Fall ein ahnliches Ziel: Das Stabilitatsverhalten 
der erzwungenen Schwingung soll in seiner Abhangigkeit von den Parametern w und ¢ unmittel- 
bar in Evidenz gesetzt werden durch Abgrenzung von Stabilitats- und Labilitatsgebieten in 
der (@, €)-Ehene. Zwar lassen sich diese Gebiete zahlenmaBig natiirlich nur fir konkrete Einzel- J 
falle festlegen. Es wird aber méglich sein, ihre Form und Verteilung in gewissem Umfang all-. 
gemein zu kennzeichnen und dazu ein Verfahren anzugeben, das die Grenzen der Instabilitats- 
_bereiche sowie die ma gebenden charakteristischen Exponenten geschlossen aus den Ausgangs- 
daten zu berechnen gustattet. Diese Rechnung beruht auf einer doppelten Reihenentwicklung. 
namlich as . 

1) der Entwicklung der Integrale von (32) in trigonometrische Reihen, 
2) der Entwicklung aller vorkommender Koeffizienten nach Potenzen von e. 


Punkt 2) 1aBt erkennen, daB die Betrachtung von nun an kleine Betrage von ¢ voraussetzt. 
Die Rechnung, die sich notwendigerweise auf endlich viele Reihenglieder beschranken muf,~ 
liefert Naherungslésungen. Die allgemeine Theorie des vorhergehenden Abschnitts kann daneben 
eine Bedeutung als Kontrollmittel behalten. ; 5 
_ Wir beginnen mit einigen Vorarbeiten zu der eigentlichen Integration der Gleichungen (32). | 
Zuerst sollen die Lésungen (42), (44) in der angekiindigten trigonometrischen Reihenform ge- 
schrieben werden. Wir setzen dazu die periodische Funktion (y als komplexe Fourierreihe? . 
& 3 ee . of ae 
Or > Cfrn e7® 
? : F ' Ui ita ae) . 
an und fihren auBerdem statt or eine GréBe h, gemaf zi aed se, 
. | a 0 =ih, ay (49) | 
ein. Dann kommt sofort aus (42), (44) ra 
— 7 +40 : : 2 
Vp Say eter 


N= = © 


ore (60) | 
TE = Sage ten. ad 
n=— 0 ie 


Die Koeffizienten der Reihen (50) sollen durch Potenzentwicklung nach é eines werden. 


eee: die Zusammenstellung am Anfang der ersten Mitteilung. Hinzugekommen ist inevehen lediglich * : 
eine Arbeit tiber den achsial pulsierend belasteten Stab: K. Kejwal, Diss. Hannover 1947, ig 

a Der Buchstabe i bedeutet die imaginare Einheit. Eine Verwechslung mit dem friiher benutzten Index i 4 
diirfte nicht zu befirchten sein. . : : die A aah Ns ‘ ; as 4 


) 
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Zur Vorbereitung dieser Entwicklungen betrachten wir die Waacaiccs fiir ¢ =0. Die ersten 
n Lésungen (33) werden fir ¢=0, wie man sofort nachrechnet!, 
cos 


r? = Zs fir f=r | 
0 Kenic 


- und somit die. Konstanten afr nach (36) 


22 w y 
%; cos — fir f= | 


Kf, — 


= (Wy 


Cf ipeale 25 tet) 


CRIA, Qeoccn) 


» Damit hat man nach (46) und (47) die Wurzeln 


22H 
ft Oe 


- und | 2 | Gr decay. . (51) 


0 Ds 2 
SO = (a) 


_ Nach (41) und (49) setzen wir deshalb 


0 =i i= aa HOr 2 Bixee (Redes TE) (52) 


@ 


Den Fall wr =x fir r+k, d.h. das Zusammenfallen zweier Eigenfrequenzen mit verschic- 
; dener Nummer, schlieSen wir von nun an aus. Dann erlauben die Wurzeln (51) der charakteristi- 
__ schen Gleichung (46) fir ¢ =0 einen ersten Uberblick ber das gesuchte Diagramm der Stabilitats- 
_ und Instabilitatsbereiche in der (w, ¢)-Ebene. Zunachst ist zu beachten, daB die Koeffizienten Ofr 
nach einem bekannten Satz der Fidercnvalaleihungsdhendis stetige Funktionen von @ und ¢ 
_sind, und da dasselbe nach einem Satz der Funktionentheorie2 deshalb auch fir die Wurzeln fe 
~ von (46) gilt. Nach dem allgemeinen Kriterium von Ziff. 6 sind die Instabilitatsbereiche in dér 
' (@, €)-Ebene dadurch gekennzeichnet, daB in ihnen entweder mindestens ein Paar komplexer 
Wurzeln oder mindestens eine reelle Wurzel, deren ess die Eins tbertrifft, vorhanden ist. 


- genommen <1, auBer in den diskreten Punkten 


O.= = “ (p> 0, ganz), a (53) 


der w- ~Achse treten auf in den diskreten Punkten 


stent St are ‘ Pane |\@, + wp| p>, ganz woe 
ees Le eng ( r+hk ). . (0) 


fir welche die r-te Wurzel 2 mit der k-ten Wurzel z,” zusammenfallt. Das Absolutzeichen muB 
an a» (RED: stehen, weil die Frequenz. mw ihrer Natur nach positiv ist. 


Die etwa vorhandenen Anstabilitatsbereiche kénnen also wegen der stetigen Abhangigkeit 
der Wurzeln z, von @ und ¢ nur in den durch (53) und (54) gegebenen diskreten Punkten auf die 
»-Achse paftreticn’. In allen anderen Punkten der w-Achse hat (46) n verschiedene Wurzeln 
mit Betragen <1, so daf hier und auch bis zu einer endlichen Héhe e>0* dariber Stabilitat 
herrschen mu. Uber den Punkten (53) dagegen wird die Wurzel z, vielleicht in einem gewissen 
B reich tiber Eins ansteigen. Wir nennen éinen solchen Bereich, falls er existiert, ein Instabili- 
sgebiet erster Art. Ebenso ist mit der Méglichkeit zu rechnen, daB iiber den Punkten (54) 
_konjugiert komplexes Wurzelpaar z,, x, auftritt. Derartige Instabilitatsbereiche nennen wir 
n solche zweiter Art, oe die” besonderen Falle, die eintreten, wenn gewisse Werte oe und 


ie uablidehoreie erster: und ae Art sind tbrigens gemeinsam dadurcty gekenn- 
net, daB der zugehérige Koeffizient h, in ihnen komplex wird, wahrend im Stabilitatsgebiet 
hr see sind. Sat ae ae das. sofort aus ater aus (41), (47) und: ae) jelpreden Gleichung 


' Der Tndee ‘nied (0) ‘ik stets ag Fall ¢=0 bezeichnen. 
ae Knopp, ‘Eunktionentheorie, Bd. II, Kap. V, 2. Aufl. Berlin u. Leipsig 1926, 
paths: eenust selbstverstindlich, ese € WU. Betirohten, a 


Auf der w-Achse sind die Wurzeln von (46) nach (51) alle reell, und sie sind tiberdies. absolut 


i denen ae r-te Wu re) re (+ fir gerades, — fir ungerades p) wird. Doppelwurzeln auf 


\ 
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% =cos 27h, oder auch unmittelbar aus der Form der Lésungen (42), (44), die nur dann eine 
Exponentialfunktion mit reellem Exponenten als Faktor abspalten, wenn 6, einen nicht ver- 
schwindenden Realteil, mithin h, einen nicht verschwindénden Imaginarteil besitzt. Diese Fest- | 
stellung wird fiir unsere Reihenentwicklung von Bedeutung sein. Auf den Grenzlinien der In- 
stabilitatsbereiche erster Art liegt eine Wurzel z,=--1 oder —1 vor, auf den Grenzen der Bereiche 
zweiter Art eine Doppelwurzel z=; hier sind also im allgemeinen ausgeartete Lésungen mit — 
Sakulargliedern zu erwarten. Man sieht jetzt, warum diese Losungen kaum von Interesse sind: | 
sie zeigen zwar Instabilitat der Grundbewegung an, aber diese Feststellung ist gerade auf der | 
Grenze belanglos. 

Damit kénnen wir die Vorarbeiten zur Reihenintegration von (32) beenden. Die Integration 
soll sich auf Parameterwerte w, ¢ beschrinken, die in der (0), e)-Ebene in der Nahe der diskreten 
Punkte w, der w-Achse liegen. Also wird insbesondere w in der Nahe der Grundwerte @) bleiben, 
und man kann Fym, das eine Funktion von q@ ist, nach Potenzen von w—q@, entwickeln: 


Fem = FP, + (@ — Wy) ER -+ (@— @)? Fan +: :- ; : : (55) 


Die Fe sind darin bekannte Funktionen von w,. Ferner dividieren wir (32) mit @> durch und 
fihren zur Abkirzung die Bezeichnungen 


; Opi OS AG go Sse q 
Wo = Gr> Oo q> po h, q: . (56) i 
FO FH ) 3 
bc imeey al() ee me ad 2) 2) 72 ¢ ite) 
op, “BOR, BR=GR, FRo.=6R.... 6M 


wt 


ein. SchlieBlich ist es klar, daB die Bestimmung von Vj allein ausreicht, da ja Vy, aus Vj, eiaireh 
durch Umkehr des Vorzeichens von s entsteht. Es geniigt also, die Differentialgleichungen fir | 
vg = Vr zu betrachten. Sie werden nach (32), (55), (56) und (57) (a | 
BV, ak 7 
gH + GV + 


“$2008 8 >’ Vn (G® +.(g—1) GR + (q—1)26BR +++) =0 (f=1,2...-n). 
f m=1 $ ; ' } 


(58) 


\ 
ie! 
‘ 


Die GréBen g, q und cf,,, werden nun allgemein als Potenzreihen von ¢ angesetzt. Wir ver- 
langen, daB w fir e=0 in wy tbergeht, um den im Punkt @, auf die w-Achse auftreffenden | 
Instabilitatsbereich erfassen zu kénnen. Folglich mu8 die Reihe fir q nach (56) mit 1 beginnen. — 
Ferner wird h, fir ¢=0 nach (52) und (56) gleich gq, so daB die Reihe fiir q mit g, anfangen mu8. © 
ae eres der Reihe fur cy, bleibt vorerst noch offen. Unsere Potenzreihenansatze | 
auten also S ct : | 


Siem = Sen “he Cnn +e Chen eae is 4 
: i q . ! 6 gO a? g@) oe sy at Te \ i (59) 
NR eGe ig = @ be gt) tatg os, | 


ie Sie sind noch durch eine Aussage uber das Anfangsglied g, der Reihe q zu erganzen. Wenn @, 
einen Wert (53) oder einen Wert (54) mit dem oberen Vorzeichen (+) bekommt, so ist nach (56 
Ve ecatla cen a Sas Te aee qe + Oe = Pp (r=k oder r+k; p>0, ganz). (60) a 


* Nimmt | dagegen @y einen Wert (54) mit dem unteren Vorzeichen (—) an, Siiaaea nach (56) 


ee Cee apes k= P (r+k; p>0, ganz) ‘ (ewer fed (61) 


Wir haben in (61) allerdings das Absolutzeichen von (54) nicht ibernommen und damit im — 
__Falle des Minuszeichens die Voraussetzung gq, > qx oder @r> Wk gemacht, die offenbar mit i 
_ einer entsprechenden Numerierung der beiden durch (6l verbunderien’ Kipentrequenzen ohne 

_ Beschrankung der Allgemeinheit zu erfiillen ist. } Tiger ek he ae ees 


Va Zum Verstandnis der weiteren Ableitungen ist. es ferner wichtig, au beachtens dap: daveh.: 
die Wahl eines bestimmten @y nach (53) oder (54), d.h. eines bestimmten zu untersuchenden 
_Instabilitatsbereichs, drei bestimmte Ziffern r, k und p herausgestellt werden und gegeniibe: ‘ 


er ae | 


5 


den veranderlichen Indizes f, m, 7 vorlaufig festzuhalten sind. Der Index r ist gleichzeitig die 


Nummer des Partikularlésungssystems Vr, das ‘im folgenden berechnet werden soll. 
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8. Das Rekursionsverfahren zur Berechnung der Reihenglieder. Fiihrt man zunachst die 
Reihen (50) in (58) ein, faBt die einzelnen Glieder e'(,+”)* zusammen und setzt ihre Faktoren 


gleich Null, so folgt 


q+ 2g Cfrn — | 
L(g ”4) G | ern (fa bic2ee heres 


. 62 
BD Cmny—1 + Omang) (GQ + (g—1) GR + (@=1)? GR + ---) <0 ee Ce ae (02) 


In diese Gleichung gehen wir mit den Ansitzen (59) hinein und erhalten durch Vergleich der 
Potenzen von ¢ nach kurzer Rechnung die Rekursionsformeln zur sukzessiven Berechnung der 


Reihenglieder 
[ia +1)? —g] o?, = 0. | (63,0) 


L(g. is m)* —q| CF n =—2 (qr a, n) (q® at nq”) CP n re 


Pa (63, 1) 
| ? og 1 Ds Gr, (eo. n-1 i Gre Lt) ; 
; [(ae +)? — af] ¢f?, =~ 2 (ae +) (G® + 4%) ef, — 
| : — [2 (@ +) (G+ 4) +(G@ + q®)?] , + 
yates 
ex SOR 2s Groth afta Oneal e (082) 
ge s Go (c + ) 
) col fm \“mr,n—1 mr,n+1/? 
lw + nt-gl ef, - KY, 
- mit ; é ; 
KP, Be prue Se oe SS (q® + nq) (q@-9 + nq —*)) —- 
oe t=1 &=0 | 
‘s ede ja ne 63,3 
: op Sraleen na: + Gaal ee 
m=1 
zi 1 j-l1 on ; 5 4 p : pays 
epee 2 FAT (Cl sf as Cima 
t=1 m= 


D (63,j) ist q =1, gq =q, zu setzen. Die GroBe ae) stellt den Koeffizienten von €” in der am 
chluB von (58) stehenden Reihe (q—1) G2 +(q—1)? Gin +--+ dar. Es ist also nach (59) 


TR = op, 
. TR = gp +a ew, . 
TS = q® GD +2 g® ¢9G2 + go cw, 


\Ilgemein ist J") ein Polynom Siler q®) mit St. Die Folge der Gleichungen (63) setzt sich mit 
apt); 63,j+2),... unbegrenzt fort. ide r ees wae | 
Damit das durch (63) vorgeschriebene Rekursionsverfahren eindeutig wird, miissen wir noch 
einige F estsetzungen treffen. Die Anfangskoeffizienten c¢;), sind nach (63,0) nur dann von Null 
ee : fst cape Ol et | (64) 
Mind Biches damn taberlaunt. Die Beziehung (64) ist nach (60) und (61) erfiillt ay ae 
SF : cers | dnd ede | a (65) 
i ata ee oa Cay a 
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Wir verbieten ferner die Erfillung der Gleichung (64) fir alle andern als die in (65) zusammen- ; | 
gefaBten Ziffernpaare f, 7. Dann wird 


Ore 
fry Zh 


fir alle f, 7 auBer den in (65) angegebenen, wahrend c¢??, und cy, _, unbestimmt bleiben. Uber |] 
einen der letztgenannten Koeffizienten darf beliebig verfiigt werden. Wir setzen deshalb 


Crowd | 
und, denken uns dainit die Lésung (50), die ja nur bis auf einen konstanten Faktor bestimmt ist, | 
mit c,,, durchdividiert. Dagegen darf man, um die Allgemeinheit zu wahren, Sas ie nicht auch = 1 | 
annehmen; wir schreiben daher vereinfacht 


(9) =n. (0) 
Cher; =p =-¢€ 


und miussen diesen Koeffizienten als vorlaufig unbestimmte Konstante in der Rechnung mit: | 
filhren. 
Als. Zwischenbemerkung sei hier kurz ndtivtt, was. das eben ausgesprochene Verbot der +) 


Beziehung (64) bedeutet. Nach (53), (54) und (56) 1aBt sich (64) in 
“(@r+ ox) + p(@r+wr) =0 (ky f,r=1,2,...n; p,7n ganz) (66) if 

umformen. Damit sind also gewisse lineare Gleichungen mit ganzzahligen Koeffizienten zwischen | 

je drei Eigenfrequenzen wx, wr, we ausgeschlossen. Dieses Verbot miissen wir in Erweiterung der | 


schon oben eingefiihrten, Forderung w,+qax fir r+k aussprechen; seine Verletzung wurde in | 
den folgenden Entwicklungen die Bildung der Reihenglieder unméglich machen. 


Nachdem iiber die Koeffizienten c/?, in der angegebenen Weise verfigt ist, lassen sich: die 
Koeffizienten de, mit j>0 durch sukzessive Berechnung der rechten Seiten von (63, 1), (63,2), . 
unmittelbar biscimen, soweit ihre Indizes f, 7 von (65) abweichen. Dagegen steht in (63) fiir 
die Ziffern (65) links Null, weshalb fir eben diese Ziffern f,.y die Ky verschwinden miissen: | 
und die entsprechenden ee. wiederum unbestimmt bleiben, und ene far alle j- Da aber die 
Tesung: Mis Wut C3 durchdividiert sein soll, muf 


3.0) furs 4,>0 
gesetzt werden, wahrend die 
Se (7 = 0) 


als zunachst freibleibende Konstante in der Rechnung mitzufihren sind. Das Verschwinden. ] 


_ der KO fir die Ziffern (65) 1aBt sich durch grelgncte Bestimmung der q@ und q) erzwingen. 
Die Gletchungen 


KW), = 0 ‘(erste Hadneglcithune); 


KW) 


ies (67): 
kr,-p —9 (zweite Hauptgleichung) ( ) 


sind somit als Bestimmungsgleichungen fir gq) und q\) zu benutzen. 

‘Damit ist das Rekursionsverfahren eindeutig vorgeschrieben. Es liefert durch sukzessive | 
- Auflésung von (63) und (67) die Koeffizienten eo » q2) und qU) und damit die Lésungen ee in | 
z Abhangigkeit von den -vorlaufig noch frcibleibenden Parametern e(), olen 


9. Entwicklung de Hiauptgleichnnsen,. Es aitaiat im a Polgevden Kaipmichinr darauf. an, 
die durch (67) eingefihrten Hauptgleichungen naher zu studieren und daraus die Grenzen der | 
gesuchten Instabilitatsbereiche sowie die GréBe der charakteristischen Exponenten abzuleiten. 

-Dazu miissen wir zundchst den Mechanismus der Gliederbildung in unserem Rekursionss 
verfahren naher betrachten und bedienen uns dazu des in “Abb. 2 wiedergegebenen Schemas. | 
Wir betrachten ein Quadratnetz der skizzierten Art und beziffern darin die Zeilen von der 
oberen Begrenzung an nach unten mit j und die Spalten von einer beliebigen mittleren Spalte 
_ aus mit 7 (nach rechts positive, nach links negative 7). In das quadratische Feld, in dem sich | 

die y-te Spalte und die j-te Zeile kreuzen, denken wir uns die n Koeffizienten cl) (f=1, 2, .. .n)| 
eingeschrieben. AuBerdem soll in der Spalte 4 =0 der Ausdruck KY) 28 ) 


rro. (die linke Seite der ersten | 
Hauptgleichung (67)) und in der Spalte 7 = —p der Ausdruck KV) - Ly (die linke Seite der zweiten | 


Hauptgleichung (67)) gusatzlich in die. j-te Zeile eingetragen werden. SchlieBlich setzen wir 
noch fest, daB bei der Ausfillung der Tabelle durch sukzessive Auflésung von (63,1), (63, 2) 


x mach cf,» ef), 


fry? Chr n? Au We) erganzende Roe der: See mis ~ nach a und 4 Bs (i). vor- 


NY 


bt So 
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—<+— 7 ——_> + erst noch zurickgestellt werden 


Be een cen ante soll. Vielmehr sollen simtliche 
1 Az Gr6Ben explizite in ihrer durch 
LY (63) gegebenen’ Abhangigkeit 
\ von den gq) und gq) einge- 
SE + tragen werden. Dann _ lassen 
sich aus dem Schema sofort 
BRRCUAHRERORGCROnEM die folgenden  Feststellungen 

I TL ablesen: 
1) Die Ausfillung der Ta- 
Soo So oe Ua uaeneeBReeuE belle erfolgt Zeile fir Zeile von 
oben nach unten. Die héchste 
[GY HH Zeile 7 =0 enthalt in dem Qua- 
Whee THEE drat 7=0 allein die GriéBe 1 
EN und in dem Quadrat 71 =—p 
J LG) ELL nur die GréBe c.. Alle andern 
GY 1 ' Quadrate der Zeile 7 =0 sind 
LULL LL) leer. Ein beliebiges Glied an 


irgendeiner Stelle der Tabvlle 


Ke), . wird, weil in einer Gleichung 

(63,j) neben dem Index y nur. 

) : : noch die Indizes y—1 und y+1 
A 9 ESO EAEaCe _auftreten, nur aus solchen Glie- 
g Al Platze Wes dern erzeugt, die lotrecht tber 
MN Keofizienten 36 ihm und schrag uber ihm auf. 


den 45°-Linien stehen. Die 
Tabelle enthalt daher Eintra- 
gungen nur in den_ schraf- 
fierten Gebieten. Die weil gelassenen Quadrate bleiben leer. Die linken Seiten der Haupt- 
- gleichungen (67) haben ihren Platz in den beiden -ausgezeichneten Spalten 4=0 und 7=—p. 


Die vereinfacht. mit cli) bezeichneten Koeffizienten cy) p stehen in der Spalte 7 =—p. 


Abb. 2, Tabellenschema zur Berechnung der Reihenglieder.. % 


2) in dem: lotrecht schraffierten Gebiet gehen alle Clieder von dem Anfangsglied 1 aus. 
Der Anfangskoeffizient c) kann nach dem unter 1) Gesagten keinen Beitrag liefexa, Gipasowents 
einer” der tiefer-stehenden Koeffizienten c(). Alle Glieder-in dem lotrecht schraffierten Gebiet | 
sind daher frei von allen cl) ( = 0). ‘Umgekehrt tritt in dem waagerecht schraffierten Gebiet 
kein Glied ohne ein c® auf. In dem schrag schraffierten Gebiet stehen Benispit Glieder, die 


at ad 


aus Teilen mit CO) und aus Teilen ohne (8). gusammengesetzt sind. . ee 


ee Die unter 2) genannte Abhangigkeit der Tabellenglieder von et (6) ist ioe Das ist 


ne mr,n— 


, linear sind. Eine Aussage iiber die Faktoren der cS) in den Tabellengliedern aoe 


ae ef). eae ef) riwie r(c= OT 2, -.J)> multipliziert mit bestimmten Faktoren. 
van Mapsbel, die Cleans (63, i+é) é > 0, aaa so oe man darin dieselben_ 


fSt 9) 


Cnr, n— 1? 
; Das: laBt sich durch-eusfuhrliches ‘Anschreiben von (63, j) ad (63, jay’ Neicht veri- — 

und folgt im tibrigen a aus dem Umstand, da® (63) durch Ausmultiplizieren der Potenz- ae 
) in den Gleichungen. (62) entstanden ist, Wenn also der héchststehende Koeffizient — 

Verla r Rekursion i in dem unter ihm liegenden’ Gebiet der Tabelle, begrenzt durch die 
beid ihn ausgehenden 45 °-Linien, bestimmte Anteile in den Tabellengliedern mit dem — 
Pak ei zeug ul der um & Zeilen tiefer stehende Koeffizient c an entsprechender 2 


+e 7 


-gegentib dem obengenannten_ um & Zeilen tiefér gertickten Gebiets dieselben Glied- 


urch Berechnung aus den Gleichungen (63, J), (63,2) 2) .. .(63,j) eis 
nde. weitgenannte ‘mit dem Faktor (8) wegen der angegebenen 

In in’ | genau dem gleichen Rechnungsgang aus den Gleichungen 
1 anderen Worten: die Glieder i in einem beliebigen Feld ‘ 


1c , selbstverstindliche Folge des Umstandes, | daB die Gleichungen (63) in den oe FG Ree bates? 


ae Note bemerkenswerten Eigenschaft der Gleichungen (63) her. In (63,j) stehen : 


en, aber etzt mit dem F aktor cS) statt c, Ist namlich einer der erstgenann- ent 
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des waagerecht und des schrag schraffierten Gebietes der Tabelle sind lineare Funktionen der ~ 


Parameter c, c®, c2), ...cW), wobei die Faktoren von c(S) gleich den entsprechenden Faktoren 
von c in dem um é Zeilen hoher liegenden Feld sind. ; 
Auf Grund der Gleichungen (63) und der eben getroffenen Feststellungen 1) bis 3) kann 


nun der Bau der Hauptgleichungen (67) genauer angegeben werden. Wir stellen zunachst die 
Ergebnisse zusammen und bringen anschlieBend die erforderlichen Begriindungen. 


Die erste Hauptgleichung lautet 


Bw) fir O<j<p, (68) 
2 gr qu) = BY) + x cl far ‘fj =p, (69) 
BO+x2cU-P) fir j>p, (70) 
und die zweite Hauptgleichung wird 
DU) fir <p <p (YG 
CE) 2 ge (GF — pal) e = m oO) + x far j =p, (72) 
DG) + Dl) cli-P) 4) 2 q, (q”— pq) ci-P) fir j > p. (73) 


Darin sind die BY) und DY) in ihrem analytischen Aufbau 


rationale Funktionen der Gleichungskonstanten Gene Ge» een Fa 


Polynome von q, q(, Gey Geto qe), 

fur j =p frei von allen cS) (E20), 

fiir j > p linear abhangig von c®), c, ... climP—)), 
° far j =1 gleich Null. 


Die Konstante x wird unten angegeben. In (71) bis (73) treten doppelte Vorzeichen auf. Die 


oberen, nicht eingeklammerten Vorzeichen stehen, wenn zwischen q, und q, die Gleichung (60) 


gilt, wahrend die unteren, eingeklammerten Vorzeichen aus (61) folgen. Stets, wenn von nun 


ab das Doppelvorzeichen (+) oder (+) auftritt, bezieht es sich in derselben Weise auf die Glei- 
chungen (60) und (61). 


Zur Begriindung dieser Behauptungen genigt der Hinweis auf einige wenige Punkte. DaB 
sich in (68) bis (70) 2 gq) und. in (71) bis (73) (+) 2. q (q4) — pq) c auf die linke Seite des. 


Gleichheitszeichens bringen laBt, erkennt man sofort, wenn man die Ausdriicke KY. bzw. 


rro0 


‘ Ki) _p nach (63,j) unter Beachtung von (60) und (61) ausfihrlich hinschreibt. Fur j<p folgen 


‘dann die linken Seiten von (68) und (71) unmittelbar aus der Tatsache, da® K), in der Tabelle 


Abb. 2 in dem lotrecht schraffierten und KY) _, in dem waagerecht schraffierten Gebiet steht. 
Gleichung (71) lautet nach der vorhergehenden Feststellung 3) tiber die Eigenschaften der 


Tabellenglieder ausfihrlich 


(#) [2 "Tk (q pets pq) c(9) + 29% (qu-D ees pqi-) ce) aie cee 2 Tk (qv) — pq) ci) 
F : ‘ = DU) ¢ + DGD) ce. + D@) eG-), 


Da sich hierin aber die Glieder mit c(i—), cU-2),... auf Grund der vorhergehenden, d.h. fur 


eae, 


j=l, j=2,... aufgestellten Gleichungen (71) herausheben, bleibt die Gleichung allgemein in 4 | 


der Form (71) stehen. 


Fur j =p tritt in (69) ein Glied xc und in (72) ein Glied % ohne Faktor c(®) auf, Die Platze § 
dieser Glieder in der Tabelle Abb. 2 sind die doppelt umrandeten Quadrate, und man findet die 


GréBe des Koeffizienten x offenbar dadurch, da8 man von den Anfangsgliedern c) bzw. 1 aus- | 


fihrt, hier angegeben. 


gehend nacheinander die mit dem Faktor c( bzw. 1 behafteten Anteile derjenigen Tabellenglieder | 
berechnet, die in den Feldern auf den beiden gestrichelten 45°-Linien liegen. Da es auf den — 
_ Koeffizienten x im folgenden besonders ankommt, sei die Teilrekursion, die zu seiner Ermittlung 


Wir behalten in dieser Rechnung der Rinfachhett halber die Bezeichnung © bei, oe : 


frm 


co) ,=1 und schreiten auf der in Abb. 2 nach links unten laufenden gestrichelten Diagonale — 


fort, haben also nacheinander cep, a cr, _9)... Zu bestimmen. Aus (63) folgt sofort 


die angeschriebenen GréSen nur Koeffizientenanteile. darstellen. Zuerst gehen wir aus von ‘| 
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(0 
t ® 46 
fils eke 2 2)? 
(qr —1) ak Fe 


(G) met 
ey (qr —j)? —¢ 
Nach p—1 Schritten ergibt sich ie Anteil von K(?) |, der aus den berechneten GréBen zusam- 


mengesetzt ist, zu 


nh 
72.8 GO) (e-)) 
km Cnr, —(p-—1) 
oder nach Rickwartseinsetzen der vorstehenden Ausdriicke zu 


n GO-EOr ere) ral) 


rry ~TyFro "p—2"p—1 Ta 


=1 [(a—1)?? — 4} + [(-@—9)? — 47-4] 


k 


(74) 


_-Nehmen wir zweitens die AusgangsgréBe ef) Se =¢() vor und gehen von ihr aus auf der 


anderen in Abb. 2 gestrichelt eingetragenen Diagonale weiter, so kommt entsprechend 
16. <0 


aces 
Pe OO) (gq (p-= Ya? 


: = (0) (j—1 
in : z a, Co mr, —(p—j+1) 
eo: Oy a ee sae 5 ’ 
$e) Oe 


und nach p-l Schritten der rogehige Anteil von n K@) 
: . | 0) ¢f 
25: 2 Grn en —1? 


: _woraus “sich durch Rickwartseinsetzen und Herausheben von c() 


1 li Bnet 0) G0), 0) AO) 
Rees de ae Ge, Ce ls MY a pec nor 48) 
Peace pee es We iP) ah): fpr teet BAY 
She. Die Indizes rs Biv ee pian der Summe (74) durchlaufen ebenso wie die Indizes f,, fy, . - ay Ney 


n der Summe (75) unabhangig voneinander die Zahlen | bis n. Zu jedem Glied von (74) mit der 
Indexfolge r4,+--fp—1gibt es daher in (75) ein entsprechendes Glied mit der Indexfolge f, =rp_1, 
te STp—2-- aes Diese beiden Glieder haben gleiche Zahler und gleiche Nenner. Die beiden 
Summen (74) und (75) sind also gleich. Erganzend sei noch vermerkt, da® sich x fir p=1 nicht 
nmittelbar aus (75) baw. ie entnehmen laBt. Aus der Ableitung dieser Gleichungen folgt 
aber dogleigh 


oe Snake fiir p=l. eee ee meee (76) 


SchlieBlich ist aa den F all J = pw betrachten. (70) ub. dineae cli-P) enthalten mit dem- 
I Faktor w, den c® in (69) besitzt. Ferner zeigt die Stellung von K\, in der Tabelle 
2 an, da® BU) auBerdem cc), ... cG-P-) linear enthalt, aber kein “0 mit &2j—p. 
8 in (73) cU-P) mit dem pleiohen Faktor multipliziert sein wie c in (72). Daneben 
Di als. lineare F unktion von c, ec, ...cU-P-), cG-P+), .. .c). Darin ist der Faktor von 
‘aktor von. ¢ i in der um & ‘Selivitte orhernenenden Gleichung. Fir &>j—p 
DG-8(4)2 gq. (4-9 —pq—®)) und verschwindet auf Grund von (71). 
are Funktion yon ¢(, c(), ...c-P-)) stehen. Endlich folgt BO = DO = 0 
ist die oben angegebene - orm. der Hauptgleichungen wohl aus- 


Wises 


ry 


29* 
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Um gq) und q?) aus den Hauptgleichungen getrennt zu erhalten, setzen wir weiter qu) 


aus (68) bis (70) in (71) bis (73) ein. Dann kommt statt (71) bis (73) 


E® far j= 'p ; (77) 
; qk fete Bice 
Bt, puede |S ol (cey <a| fir =p. (78) 
Bis apes Gate 
EG) + 2 ci?) ee (F) ao] firi>p- ee | 
mit 
Do) fir j=p,; 


7 o 
i Steebae +4) pura far j > p. 

Nach dem oben tiber die BY) und D& Gesagten sind die Ms J) und gq) in der durch (68) bis 
(70) und (77) bis (79) gegebenen Form Polynome der g und q©) mit § <j. Diese Abhangigkeit 
lat sich offenbar durch sukzessives Einsetzen von q@), q® in q@), ¢@), von g®, g@ und qg@), g@) i 
in g®), q@) usf. eliminieren. Die Gleichungen (68) be (70) ssid (77) bis (79) behalten a 
diesem Verfahren ihre auBere Gestalt, aber ihre rechten Seiten bekommen einen anderen ana- 
lytischen Bau: 

Die BU) und EY) werden fiir j Sp Konstante, die nur von den Gleichungskonstanten Gc, 
.. und gq abhangen, und fiir j > p Polynome von 1/c, c, cl, . ..c(i-P-}), Deshalb sind 
die q() und q\) fiir j <p Konstante, wahrend sie fir j =p Polynome von 1/e(, c, c®, .. ci?) - 
werden. Die Abhangigkeit der q'), ¢) von cl ist durch (69), (78) explizit angegeben, ebenso | 
die Abhangigkeit der g), q) (j > p) von cU—P) durch (70), (79). Auf diesen Punkt kommt es. | 
im folgenden besonders an. Bemerkenswert ist noch, daB wegen BO) = DY)—0 auch E)=0- 
ist, woraus gq = q@=0 fir p> 1 folgt. Al 

An diese umgeformte Gestalt (68) bis (70) und (77) bis (79) der Hauptgleichungen knipft 4 
die weitere Diskussion an. ; 


RN 


_ 10. Bestimmung der Grenzen der Instabilititsbereiche und der charakteristischen Exponenten. |f 
Die Funktion q=q(e)=l1+eq%-+--- stellt geometrisch gedeutet eine Kurve in einem (gq, €)- — 
Koordinatensystem dar. Durch Multiplikation von q() mit w, erhalt man nach (56) eine Kurve — 


Wir sprechen im folgenden kurz von der Kurve q. 3 
Die Seah elena Ziff. 9 brachte zum Schlu8 das Ergebnis, daB die q9) und q{) far j <p a 
Konstante und fir j = p Funktionen der Parameter ce‘) (¢=j—p) sind. Sind alle cl reell, jf 
so werden alle q/) und q) reell, d.h. man hat nach (59) und (56) ein reelles h,, so daB die Parti- jf 
_kularlésungen (50) nicht unbegrenzt wachsen, sondern sich stabil verhalten. Im Instabilitats- ! 
fall dagegen muB es komplexe cl) geben, derart, daB fir j =p die q'/) komplex und die q\) — 
gleichzeitig reell werden. Dann ist namlich nach (59) ¢ und nach (56) hr komplex, so daB also 
die Lésung (50) einen anwachsenden Faktor abspaltet, mithin die zugehérige reelle Kurve qin — |} 
einem Instabilitatsbereich liegt. Wir stellen die Frage, ob dieser Fall méglich ist. 


Zur Pourrorane knipfen s wir an (78) an und schreiben darin ¢( allgemein als komplexe | Zahl 3 


. ce) = R(cos p+isin g). 
Sie liefert in (78) eingesetzt den Inidginarted von 2 qupq”) . a 


‘ LAS ; TRapd) =x|2 Ry zIsing, 
Fee ifereelaindlel mu, wenn q() reell bleiben soll. Er ist Null eet 


ue | entweder fiir R? = cy = y heliebig 
oder fir mp =0, RB beliebig. 


ie Zweite Méglichkeit fibrt allerdings vom Wege ab, denn sie bedeutet reelles O, 


Wir betrachten nun zuerst den Fall, daB (61) erfillt ist und somit in den vorstehenden | 
Se arainh das eingeklammerte, untere Vorzeichen gilt. Da die Gleichung R? =— qr/q, durch kein — 
-reelles R befriedigt wird, kann c() nicht komplex sein, sondern es liegt jetzt der Fall eines reellen 
fe : ~€ vor. Dann miissen aber nach (79) auch cl), c®), ... reell sein, damit kein oe ae 


@=w(e) =p q(é ) in der (q@, e)-Ebene, dem Ausgangs- und Crake unserer Betrachtung. — | 
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_ auftritt. Dadurch ist die oben gestellte Frage verneint. Man kann den q) fir j= p beliebige fae 

Werte vorschreiben und die Gleichungen (78), (79) sukzessive durch gecignete reelle Zahlen eee 

- c®, cM, ... stets erfiillen. Die in dieser Weise unter der Voraussetzung (61) aus den Haupt- ie 

_ gleichungen folgende Schar von Kurven gq (q\/) fiir j < p festgelegt, fiir jJ=p beliebig) enthalt Cee 

also keine einzige Kurve; auf der Instabilitat herrscht. Vergleicht man diesen Sachverhalt mit paren 

dem im folgenden fir den Fall der Gleichung (60) diskutierten, so wird man ihn ohne Bedenken ; 
- dahin deuten, daB von den-Punkten Wy =(@,—@,)/p, €=0 der (w, €)-Ebene keine Instabilitats- 
_ bereiche ausgehen. Ein mathematisch strenger Beweis fiir die Nichtexistenz dieser Bereiche 
liegt allerdings nicht vor. Er ist mit unsern Hilfsmitteln, namlich der Aufstellung rein formaler 
_ Reihenentwicklungen ohne Konvergenzuntersuchung, auch gar nicht méglich, ebensowenig wie 
die folgende Berechnung von Instabilitatsbereichen auf Grund von (60) als strenger Existenz- 
_ beweis gelten kann. Wir lassen aber den Fall (61) im vorliegenden Abschnitt von jetzt ab beiseite. 


ee Anders liegen die Dinge, wenn wir die Gleichung (60) zugrunde legen und somit in allen 
_ Gleichungen das obere, nicht eingeklammerte Vorzeichen bericksichtigen. Jetzt ist R= V aed 


reell und’ aac : ive 
; oye Be Nae = /« (cos p +i sin g) (80) 
-komplex mit beliebigem reellem Q- Aus (78) folgt nun i ‘ 


a Beep ge) = EO 42% 1/E cosy ae any Loe 


< 2G ge) = Be) + x y= (cos +4 sin @). 


ae (82) 


komplexes ria) bei reellem qe) aufgezeigt, was auf das Vorhandensein eines In- 
ichs hinweist: Wenn es gelingt, auch die weiteren ql) (j > p) als reelle Zahlen mit 
xen q) anzugeben, so ist damit eine reelle Kurve q gefunden, die eindeutig 


damit gestellte Aufgabe, indem wir sogleich noch_ einen Schritt: weiter — 

len Kurven im Instabilitatsbereich auch diejenigen Kurven aufsuchen, die den — | 

I s stabile Gebiet abgrenzen. Wir behaupten — und das soll nachher bewiesen — 
man die | Koeffizienten gi) d er gesuchten Grenzkurven nach der folgenden Rechen- 


ndex G an einem Formelzeichen soll dabei besagen, | daB die betreffende 


¥ 


joke det eal hfoleediion, cba tes fice Kectiienten q0? ae : 
eS pe Oe ee ae 

n iP) sind sukzessiv als Wurzeln der Gleichungen | 

De ee lan 

st aus (78) die Gleichung dg(?)/8c = 0 au bilden und nach 

von ¢') ist, folgt sofort “ oe a S51 3 Sas were a 


ung von (84) firj>p 


Rare 


ee Eline 
ew, 
ronzwe oan (Ke) 
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Die Gleichungen (84) bis (86) sind als Rekursionsformeln zur Berechnung der qi zu betrachten. 
Wenn man in (84) und (86) einmal von dem oberen und eimmal von dem unteren Vorzeichen 
ausgeht, gelangt man eindeutig zu zwei bestimmten Wertefolgen q!??, gie*, ..., die zwei be- 


stimmte Kurven 


q=qa=lteq)+erg? +: 
as) 


ergeben. Da die c{}) nach (84) und (85) reel) sind, werden die zugehérigen qd) = qi?’ gleichfalls 


reell. Die beiden Kurven qg schlieBen ein Gebiet ein, von dem wir zu zeigen haben, daB in ihm — 


Instabilitat, auBen-aber Stabilitat herrscht. 


Um diesen Beweis zu fuihren, betrachten wir neben einer Kurve qg eine weitere Kurve q, 


in deren Potenzentwicklung alle q‘J) fir j = p mit den festen Koeffizienten q?) elner der oben- 
genannten Kurven iibereinstimmen, wahrend q'?) verainderlich und durch (78) bzw. (81) 


gegeben sein soll. Wir geben nun c() zuerst einen reellen Wert + cf) und den c, c@), ... die’ 


damit festliegenden, aus (79) folgenden Werte 


| Lok pas Eee 
pelos Sree aes (J=ptl, p+2,...), : (87) 


qk 1 
“(So | 


worin die Klammer im Nenner wegen c() + cf?) nicht verschwindet. Die c&) sind damit alle — 


reell, d.h. die Kurve q liegt jetzt im Stabilitatsgebiet. Nun ist aber fur reelle c 
1 


) 


qk (0) 1-.| 
Ch . 
ai | dr e r cy | 


Te ¢(0) 
= € 
Fe Fle 


Also liegt gq) nach (78) auBerhalb des von den beiden Grenzwerten g\?) eingeschlossenen q(?)- 
Intervalles und damit die Kurve q auBerhalb des Bereichs zwischen den beiden Kurven qg 


(Abb. 3). Das Stabilitatsgebiet liegt also auBen. Zweitens geben — 
wir c() einen komplexen Wert + cl) und den cf), c@), ... wie- » 


1 


Grenzkurven % 


derum die Werte (87). Da diese jetzt gleichfalls komplex 


urve g fir Nach Gleichung (81), die fir komplexe c®) gilt, liegt q(?) inner- 
reelles oc? halb des laut (86) von den beiden Grenzwerten q!?) eingeschlos- 
senen q'?)-Intervalles und somit die Kurve q in dem -Gebiet 
g zwischen den beiden Grenzkurven qg (Abb. 3), das sich somit 


g=1 _ als Instabilitatsgebiet herausgestellt hat. Lassen wir ce) - 
streben, wobei c(°) reell oder komplex sein kann, so geht (87) — 
zunachst in die Form 0/0 tber. Durch Differentiieren des Zahlers — 


Abb. 3. Zur Herleitung der Instabilitats- - 


sa: und Nenners nach c(°) folgt aber aus (85) sofort 


lim. P= cP) (f= p+1,p+2,...). 
(0) =p 


is (0 . 7 sae ' a kee : : Tea fee 
Mit c(—¥c schwenkt also die Kurve q von auBen oder innen her stetig in eine der beiden oben 
bestimmten Kurven gg ein. Diese Kurven stellen also in der Tat die Grenzlinien eines Instabili- - 


tatsbereichs dar. 


Innerhalb der Instabilitatsbereiche, die damit in der (w, €)-Ebene iiber ded Pankten 


- Wy= (@,+,)/p nachgewiesen sind, ist, wie schon mehrfach erwahnt, q und deshalb h, komplex. | 
Fur die Starke des Anwachsens der Lésung Vy, mit s ist der Betrag des Realteils von Or oder | 


nach (49) und (56) des Imaginarteils von h, =q/q maBgebend;: dieser soll daher etwas naher be- 


trachtet werden. Offenbar ist der Imaginarteil von g/q auf verschiedenen Kurven q im Instabili- — 
_ tatsbereich verschieden groB und geht mit stetig in eine Grenzkurve einschwenkender Kurve q A 
> stetig in Null uber. Es mu8 eine bestimmte Kurve q geben, auf welcher |J(q/q)| ein Maximum | 

hat. Man kénnte diese Kurve und das zugehérige Maximum durch eine Extremalbetrachtung | 
nach den Regeln der Differentialrechnung bestimmen, was allerdings ziemlich umstandlich. : 


ware. Einfacher ist es, auf diese genaue Rechnung zu verzichten und, was sicherlich nicht viel 


fehlgehen wird, die Annahme zu machen, da das Maximum von J(h,) etwa auf der Mittellinie 


‘Kurve gy tir werden, was dieselbe Eigenschaft. fiir die q) (j 2 p) nach sich ~ 
Koamplexes c@ zieht, fallt die Kurve q nunmehr in das Instabilitatsgebiet._— 


\ ater 
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des Instabilitatsbereichs, genauer gesagt in der Mitte der waagerechten, parallel zur q- baw. 
w-Achse gefihrten Schnitte, zu liegen kommt. Diese Mittellinie 


g=lteqh +2 qi +---, 
deren Koeffizienten qj) fir j < p mit den Konstanten q')=q\) des Bereichs tbereinstimmen 
und fir j= p gleich der halben Summe der beiden Grenzwerte qi) sind, ist als bekannt an- 
zusehen, und die Aufgabe besteht jetzt darin, den Imaginarteil von q/q auf dieser Mittellinie 
_anzugeben. Wir wollen dabei, wie es an den Koeffizienten q\) schon geschehen ist, alle zur 
Mittellinie gehorigen GréBen mit dem Index M versehen. 
In der Entwicklung 
G= q+ eG) +g +... 
sind die a fir j < p durch (68) gegeben und somit reell, weil unabhangig von den cl). Nach 
(81) wird g'?) = qi fiir cos ~=9, p=2/2. Derselbe Wert y liefert EAS (82) 
=—(p) __ Be) . % 


Iu 9 t ea 
qr 2 V ana 
und nach (80) ee 
(0) gly fe: 
OMe Apa. 
Die weiteren Zahlen c{), cM, ... folgen aus (79) zu 
x : 2 qnp qi) — EY : 
(j—p) _ IME v/§ Gea a 
Cc - = 
He 2 %¢ Qk /qr Pee sl: 


Sie werden komplex, da E™) von j=p-+1 ab von c® abhangt, und ergeben in (70) eingesetzt 
| Pee | 

au 
Ti 2a a 


% 
: 2 ar 
~ Damit ist der gesuchte Imaginarteil von h, auf der Mittellinie 


cli “eb 


u 
: 2 V qk ar ; 
J (hr) 3 t+ eq) + tq?) +... 


-bekannt. Man kann noch einen Schritt weiter gehen und den Quotienten (88) durch Ausfihrun 
der Division in eine Potenzreihe verwandeln und erhalt die Entwicklung 


= Fleece as, 
cir) Maka each th ne 2) gear , 
von der aber die folgenden Glieder nicht angegeben werden sollen. ; 
act, Selbstverstandlich labt sich auch auf jeder anderen vorgegebenen Kurve q im Instabilitats-— 
_ bereich eine Potenzentwicklung von q und J(hr) aufstellen, deren Ableitung ganz entsprechend 
_ der fiir die Mittellinie gegebenen verlauft, Es wird dann allgemein 


ep + aybhe 


Pt TCG ED) ert 


(88) 


eP + eae ‘ : : \ (89) 


Ee Ey 
= ee ae 2) aKa 


ye 


. mit beliebigem @: Daraus geht hervor, daB das erste Glied der Entwicklung von J(hr) seinen 


be Maximalbetrag |x/2 V4, | eP genau auf der Mittellinie annimmt. 

= Die Entwicklungen der: vorliegenden. Ziff. 10 setzten «+0 voraus. Es ist’ mehrmals mit % 

_ dividiert worden. Sollte in Sonderfallen x =0 werden, so treten, wie hier nur beilaufig erwahnt 
sei, in der Ableitung der ausgezeichneten Koeffizientenfolgen q\J), q¥) Anderungen ein. Wir 

~ wollen hier, “wo es sich um die Darstellung des allgemeinsten Falles handeln soll, derartige Einzel- 

falle, ‘die immer nur in Sonderproblemen auftreten kénnen, nicht ausfihrlich verfolgen. Es sei 

Iediglich mit wenigen Strichen skizziert, in welcher Richtung die Anderungen gehen. Fir x =0 

sind nach (69) und (78) und dem am Schluf von Ziff. 9 tiber die BY und El) Gesagten q®) und 

q®) konstant, d.h. unabhangig von allen cl. Dies bedeutet, daB die beiden Grenzkurven qg(é) 

les Instabilitatsbereichs nicht nur wie bisher gleiche qi) mit j <p haben, sondern auch noch. 
sleiche q(?), ferner daB q) noch keinen von Null verschiedenen Imaginarteil hat, Zwei nicht 


eng 
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gleiche Grenzwerte q)?) und einen Imagindrteil von q”) kann man erst fir dasjenige p> Pp 
erwarten, welches in (70) zum erstenmal ein lineares Glied mit ¢) und in (73) ein von cl) freies 
Glied anzeigt. Der Koeffizient c) tritt dann in (70) und (79) nicht schon fir j=p+l auf, sondern 
erst fir j=p’+1, ebenso der Koeffizient c@) eine Gleichung spater usw., aber immer linear. 
Der Aufbau der Hauptgleichungen bleibt auch in derartigen Sonderfallen grundsatzlich er- 
halten, nur ist der ausgezeichnete Index j nicht mehr gleich p, sondern gleich p > p» Die Auf- 
gabe besteht dann darin, die Hauptgleichungen dem gegebenen Sonderfall entsprechend im 
einzelnen zu prazisieren, vor allen Dingen festzustellen, welche GroBen nun in (70) und (79) 
De. fir j=p’ an die Stelle des in (69) und (78) ausgefallenen x treten. Dann labt sich an die Haupt- 
gleichungen des Sonderfalls cine ebensolche Theorie zur Bestimmung der Grenzlinien der In- | 
stabilitatsbereiche und der Extremwerte der charakteristischen Exponenten anschlieBen, wie  |j 
das oben im allgemeinen Fall geschehen ist. Die Einzeldiskussion ist allerdings nicht ohne vor* 
herige Kenntnis der spezicllen Bedingungen durchzufihren. 


11. Zusammenstellung ‘und Diskussion der Ergebnisse. Die mechanische Bedeutung der Ab- 
leitungen des Abschnitts III ist-in der bisherigen Darstellung nicht ganz leicht zu ubersehen. 
Deshalb ‘sollen die wichtigsten Ergebnisse hier besonders herausgestellt und mit einigen Er- 
ganzungen versehen werden. 


Wir knipfen an die Tatsache an, daB die Zusatzverschiebungen u; des elastischen Korpers 
nach (12) als Summen von Produkten der Eigenschwingungsformen ujy mit den zunachst un- 
bekannten Zeitfunktionen vy angesetzt sind. Fiir die Bestimmung der vf steht ein System von N _| 
linearen Differentialgleichungen (11) mit periodischen Koeffizienten zur Verfigung, dessen L6- 
sungen den zeitlichen Verlauf der Nachbarbewegungen und somit die Stabilitatseigenschaften 
der Grundbewegung des elastischen Kérpers wiedergeben. Das Differentialgleichungssystem (11), — 

das in der allgemeinen Theorie als ein geschlossenes System erscheint, kann in Sonderfallen 
in mehrere Teilsysteme (32) mit n Gleichungen zerfallen. Darunter kénnen auch Gleichungen 
ohne periodische Koeffizienten vorkommen, namlich dann, wenn die zugehérigen Verschie- | 
bungsanteile an den Grundspannungen eine Arbeit leisten, die keinen quadratischen Teil A,u, 
besitzt (Ziff. 4). Derartige Gleichungen sollen allerdings als belanglos fur die Stabilitatsfrage 
auBerhalb der Betrachtung bleiben. — . . 


Die Zeitfunktionen x;(f=1, 2,...n) sind nun als allgemeine Lésung von (32) aufzufassen — 
und setzen sich als solche linear aus dem Fundamentalsystem der 2n partikularen Integrale — 
Vj, Vj- zusammen. Diese Integrale haben, wie in Abschnitt I] erértert wurde, im Normalfall 
~.die Form (42), (44), in Sonderfallen, die als Grenzfalle unwichtig sind, eine ausgeartete Form — 
mit Sakulargliedern. Der Sinn der Reihenentwicklungen des Abschnitts III ist der, in der Ko-— 

ordinatenebene der mechanisch mafigebenden Gleichungsparameter w, ¢ diejenigen Gebiete auf- 
zusuchen, in denen ein Partikularlésungssystem Vir, V7-(f=1, 2,...n) mit bestimmter Nummer r ; 
(aus den Zahlen 1, 2,...n) imstabil wird, weil der charakteristische Exponent 9, einen nicht — 
-verschwindenden Realteil bekommt, sowie innerhalb dieser Gebiete den genannten Realteil wirk- | 
_. lich zu berechnen. Wir haben also keineswegs die Berechnung des vollstandigen Fundamental- — 
systems fir bestimmte Parameterwerte @, ¢ durchgefiihrt oder auch nur angestrebt, sondern — 
_ haben lediglich die cinzelnen Partikularlésungen soweit studiert, als es zur Klarung ihres stabilen — 
a bzw. instabilen Verhaltens in der (w, €)-Ebene erforderlich war. Ret 


eee 


- 


shine Uber die Form der zu der Lésung Vrs Vj. gehérigen Instabilitatsbereiche und tiber den | 
Sis charakteristischen Exponenten Q, lassen sich einige allgemeine Aussagen machen (vgl. dazu die | 


schematische Skizze Abb. 4): 


near 1) Die, FuBpunkte der Instabilitatsbereiche auf der w-Achse haben die Abszissen = = 


wobei k neben dem jetzt herausgehobenen r eine beliebige Zahl von 1 bis n (einschlieBlich r) ist. 


Pes), Die Grenzlinien der Instabilitatsbereiche lassen sich durch Gleichungen der F orm a 
ie (91) | 


Og= Wy + eno + eta?) + see : 
onnenen Entwicklung dc= I+e ge) + e%q2)+ . a 4 


= 
4 


Py aie darstellen. Sie folgen nach (56) aus der in Ziff. 10 gew 
_ durch Multiplikation aller Glieder mit w). Jeder. Bereich hat zwei Grenzlinien (91). 


| 


oe 


= = ie ; Bias f ; ‘ 
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3) Far p=1 sind! die zu den beiden Grenzlinien gehérigen Koeffizienten oW, dh. die 
_ Steigungen dieser Linien fir ¢=0, nach (76) und (86) 


oy 
Z V Wr OL 


“Die Grenzkurven miinden also fir p=1 von links und rechts oben her unter zwei endlichen 
gleichen Winkeln gegen die Lotrechte in den FuBpunkt w, = Or + Oy ein. 


4) Fir p> IJ ist wl )— 0; beide Grenzkurven treffen senkrecht auf die @-Achse auf. Ferner 
sind die J) beider Grenzlinien bis j =p —|einander gleich, d.h. es liegt eine Berthrung (p —1)-ter 
Ordnung der zwei Kurven vor. Die beiden w?) dagegen unterscheiden sich nach (86) um die 
GréBe 2 wrx|p Vo, 0% » wobei nach (74) oder (75) 


Wane 0) 0)... yn(0) (0) 
se 4 “hy pe FE T9 supe Tp—1 Fp 
is \ Mo * = oP Sige tae =F i 
rae aief aes hep 1 ZI, [(a,—m @)? — ro 


é ist.’ as . ae ; : 
A) Der Charalsacieustba rea ‘0, 1aBt oh gleichfalls in eine Reihe 
Le =e + co + 620% + eee 


€ ntwickeln. -Er wird auBerhalb der Br clitepeies rein imaginar. Innerhalb ace Instabilitats- 
bereiche bekommt 0, einen Realteil, dessen Betrag ein MaB fiir die Starke des Anwachsens der 
Lésung Vy, baw. Vir ist, und der gegen Null geht, wenn sich der Punkt mit den Koordinaten’ 
w, € dem Rand eines Bereichs nahert. Der Exponent oe ist nach (52) stets imaginar. Fe ur p=1 


hat's schon er” einen ‘Realteil, und zwar hat dieser nach An?) (64) das Maximum 


ae Br fete A oy 1 ae ae : Sei ue nes = 40 | On or 


tbh 


IZ 


erster Art ee a Fir p sac} sind die od) bis j=p—l 


p=7 rein imaginar und erst ‘in oe) tritt: 


ein Realteil auf mit dem Maximum 


[vel oy) 


schied zwischen den ee 
bereichen erster_ und aweiter Art — 


Pi Zade Cin, den Raia 4 bis a 
gar nicht. hervortritt. WA heh 


ake ee Die Tnstabiliteuebiete fir die 
; = pares verschiedenen _ _ Partikularlésungs- — 
pi Sey at systeme. Vie Ve (r= 1, are 4; tt), die 
“aueiler Art wir uns zundchst_ far jedes r auf — 


ein besonderes” Blatt. gezeichnet ge-. 
at aa. me - alle iibereinandergelegt und el sich- 
te zeitig: betrachtet werden, da ja in 
der allgemeinen Lésung von (82) 


; i 


eihnngen 2) bow. raya aus. 


Aor Hicere aware ist, cae ate nese. ; 


_@ — dacht haben (Abb. 4), miissen nun. x 
alle ee alain eos 


a. es Stabilitat, ‘te : 
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Unter den mannigfachen auf diese Weise in dem Diagramm verzeichneten Instabilitats- |} 


bereichen fallen selbstverstandlich diejenigen mit den FuBpunkten Wy) =(@,+@,)/p, die za dem 
Partikularsystem Vj, V#. gehéren, zusammen mit denjenigen, welche dieselben FuBpunkte 
W) =(@,-+o,)/p haben und die Instabilitat des Partikularsystems Ves V fx anzeigen. Es handelt 
sich ja, wie aus Abschnitt II hervorgeht, um diejenigen Gebiete der (w, e)-Ebene, in denen die 
zu Vs, gehorende Wurzel z, und die zu Vy, gehérende Wurzel z, der charakteristischen Glei- 
chung (46) konjugiert komplex ‘sind. Es kann aber dartber hinaus der Fall eintreten, das die 
Fu8punkte von weiteren Bereichen auf der w-Achse zusammentreffen. Das hangt eng mit dem 
Verbot von (64) bzw. (66) zusammen. Die Gleichung (66) ist namlich, wie man sofort nach- 
prift, erfillt, wenn auf einen der bisher betrachteten FuBpunkte (90) der FuS8punkt eines zweiten 
Instabilitatsbereichs mit Wy = (0, +-@4)/(— n) (f£+k, 4 < 0) fallt, ferner speziell fir den Fall r=k 
auch dann, wenn der FuBpunkt w,)=(2/p)@, mit einem anderen FuBpunkt ®» =[2/(2y+p)] wf 
(f+r,24+p > 0) zusammentrifft.: Derartige Ubereinstimmungen sind also verboten. AuBerdem 
dirfen die FuBpunkte (90) auch nicht auf die Punkte w,»=(w,—qy)/(—n) (7 290) fallen, von 


denen zwar keiner der angegebenen Instabilitatsbereiche ausgeht, die aber aus rechnerischen Griin- 


den zu meiden sind. : 

Diese Einschrankungen muften getroffen werden, um die Durchfihrbarkeit unserer Rech- 
nungen in der gegebenen allgemeinen Form sicherzustellen. In den Einzelproblemen der Elasto- 
kinetik wird man sich aus praktischen Griinden auf Ansatze mit verhaltnismaBig kleiner Glieder- 


zahl, also auf wenige Gleichungen (32) mit den niedrigsten Eigenfrequenzen wy und auf kleine p — 


beschranken miissen. Damit kommt man auf nur wenige Instabilitatsbereiche niedriger ,,Ord- 


nung p, fiir welche die ausgesprochenen Verbote keine allzu starke Beschrankung bedeuten. |} 
Geht man aber zu immer gréBeren n, also zu héheren Eigenfrequenzen, und gleichzeitig zu | 


groBeren p tiber, so findet man in der Nahe der erwahnten Bereiche niedriger Ordnung eine 
groBe Anzahl von FuBpunkten (90) und von zugehérigen Instabilitatsbereichen hoher Ord- 


nung p, die wegen der (p—1)-fachen Beriihrung ihrer Grenzkurven sehr schmal sind und deren ' { 


charakteristische Exponenten einen mit ¢? multiplizierten, also sehr kleinen Realteil haben. 
~ Mechanisch ist mit Sicherheit anzunehmen, daB diese ganz schwachen Instabilitaten durch die 
hier nicht beriicksichtigte, aber in Wirklichkeit stets vorhandene Dampfung verwischt werden, 
denn anschaulich gesprochen wird sich die theoretisch geforderte langsame Aufschaukelung 


gegen das von der Dampfung bewirkte Abklingen der Nachbarbewegungen aus Stetigkeits- : 


griinden nicht durchsetzen konnen. Mathematisch aber tritt die Frage auf, was geschieht, wenn 
-n baw. N ber alle Grenzen wachst, also die Frage nach der Konvergenz des Verfahrens. 


Das Konvergenzproblem soll allerdings im Rahmen dieser. Arbeit nicht behandelt werden; — 


schon in der Einleitung wurde ja angekiindigt, da es sich hier um die Angabe eines Naherungs- 
weges zur zahlenmafigen Lésung des Stabilitatsproblems handeln soll. Dieser Naherungsweg 


grindet sich auf formale Reihenentwicklungen, namlich auf den Ansatz (6) bzw. (12), der fir. 


No in ‘Wahrheit ein Reihenansatz ist, und auf die Potenzentwicklung nach ¢ der Koeffi- 
zienten von (12). Der in der Mechanik Erfahrene wei} aber, daB fiir praktische Rechnungen haufig 


_genug Reihen verwendet werden, deren Konvergenzverhalten unbekannt ist, und da® fir die ot 
Abschatzung der erzielten Genauigkeit ein strenger Konvergenzbeweis viel weniger wichtig | 


‘ist als etwa die Beurteilung der Gré®enordnung der Reihenglieder, der Vergleich mit streng 


lésbaren Sonderfallen und ahnliche Mittel. Im ubrigen wird das Konvergenzverhalten unsrer | 


-Reihen sowohl in theoretischer als auch in praktisch-numerischer Hinsicht sehr von der GroBe. 


der nach (63) in den Gliednennern yorkommenden Faktoren (q,-+7)? —qj abhangen. Das exakte 


-- Verschwinden eines Nenners ist durch die Erfillung der Hauptgleichungen (67) und das Verbot 


__ von (64) fiir alle andern als die in (65) angegebenen. Ziffern zwar ausgeschlossen. Aber fiir die 

Konvergenz und die numerische Brauchbarkeit der Reihen ist es wichtig, da (q,+-n)? — q | 
_ auch nicht sehr klein wird, da man sonst groBe Betrage der entsprechenden Reihenglieder zu — 
_ erwarten hatte, was besonders unangenehm ware, wenn es sich um hohe Glieder (mit groBem 7 — 
- oder grofem f) handelte. Man wiinscht ja die Rechnung auf wenige Glieder zu beschranken. | 


Eine nahere Diskussion dieser F ragen soll hier, wie gesagt, unterbleiben. Erwahnt mag noch 
werden, da® unser Verfahren der Entwicklung der Lésungen von (32) nach Potenzen von ¢ 


letzten Endes mit der hauptsachlich in der Himmelsmechanik verwendeten Stérungstheorie. 4 
identisch ist, und da dort ahnliche Teiler in den Gliednennern auftreten wie hier. Diese Ver- _ 


wandtschaft JaBt es sogar als méglich erscheinen, da® unsere Reihen semikonvergent sind!. 


} Vel. R.Grammel in Handbuch d, Physik, herausgeg. v. Geiger u. K, Scheel, Bd. V, 8.358 ff. Berlin 1927. j 


t ; . 
fi . 


= 


f é i \, 2" s 5 * a od 
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Unsere bisherigen Aussagen tiber das Stabilitatsverhalten des elastischen Kérpers sollen 
noch durch einen Zusatz mechanischer Art erganzt werden. Zunachst sei nochmal wiederholt, 
dafs der Grundbewegungszustand, d. h. die erzwungene Schwingung, fiir solche Erregerfrequenzen 
instabil wird, die in der Nahe der ausgezeichneten Werte (90) liegen. Fur r=k, Wy) =(2/p) a, 
_besagt dies. da dic erregende Frequenz @ in die Nahe der doppelten EFigenfrequenzen und ganz- 
zahliger Teile davon fallen muB, damit Instabilitat entsteht. Diese Erscheinung hat sich bei 
allen bisher im Schrifttum behandelten Einzelproblemen der Stabilitat elastischer Schwingungen 
gezeigt und nunmehr allgemein bestatigt. Wir sprechen hier von Instabilitatsbereichen erster 
Art. Neu ist aber die Tatsache, da8 auch dann Instabilitaét auftreten kann, wenn die erregende 
Frequenz @ indie Nahe der Werte w)=(@,-++@,)/p mit r+k kommt, also ganzzahliger Teilbetrige 
der doppelten Schwebungsfrequenz zwischen der r-ten und der k-ten Kigenschwingung des 
elastischen Kérpers. Derartige Instabilitatsbereiche haben wir als. solche zweiter Art bezeichnet. 
Sie kommen hei denjenigen Einzelproblemen nicht vor, in denen synchrone Nachbarbewegungen 
aller Teile des Kérpers méglich sind mit einer Schwingungsform gleich einer beliebigen Eigen- 
schwingungsform. In diesen Fallen zerlegt sich das Gleichungssystem (11) in eine Folge von 
unabhangigen Mathieuschen Gleichungen, die bekanntlich nur Instabilitatsbereiche erster Art 
liefern (Beispiel: Der beiderseits gelenkig gelagerte, axial pulsierende belastete Stab). Im 
allgemeinen Fall dagegen sind die Gleichungen (11) mindestens zum Teil miteinander gekoppelt 
und fihren dann auf die genannten Instabilitaéten fiir die Schwebungsfrequenzen!. 


Zum SchluB dieses Abschnitts noch eine Bemerkung: Wer das verhaltnismafig komplizierte 
Verfahren verfolgt, das uns in Ziff. 7 bis 10 zur Abgrenzung der Instabilititsbereiche gefihrt hat, 
wird sich vielleicht an die elegante Methode erinnern, mit der man in der Theorie der Mathieu- 
_gleichung die Instabilitatsbereiche durch Berechnung der periodischen Lésungen festlegt?, und 
wird sich fragen, ob diese Methode nicht auch in unserem allgemeinen Fall zu einfacheren Ent- 
wicklungen fiihren wirde. Dies ware in der Tat so, wenn wir uns nur mit den Grenzen der In- 
stabilitatsbereiche erster Art beschaftigen wollten.. Auf-die Bereiche zweiter Art laBt sich die 
-Methode aber schon nicht mehr iibertragen. Auch die Berechnung der charakteristischen Ex- 
“ponenten leistet sie nicht. Wir haben es deshalb vorgezogen, hier das gegebene zwar kompli- 
ziertere, aber umfassende Verfahren allein zu entwickeln. 


- 12. Fortfiihrung des Beispiels. Die in Ziff. 4 begonnene Behandlung eines Beispiels soll nun 
soweit fortgesetzt werden, daB die Konvergenz und zahlenmafige Brauchbarkeit des im vorliegen- 
den ‘Abschnitt entwickelten Verfahrens wenigstens in einem Einzelfall deutlich wird?. Fir die~ 

“numerische Durchfihrung brauchen wir zunachst bestimmte Zahlenannahmen und wahlen dazu 
die Werte (vgl. Abb. 1) 1=100 cm, h=10 cm, b=1 cm, ferner die Materialkonstanten fir Stahl 

E=2,1-106 kg/em?, G =0,81-108 kg/em?, us =0,800- 10-4 kg sec?/em? . 


» Damit craibt sich* ~~ . ; : | 
ge) oa a 0.8038 tm4, - +: ty 2,901 cms | 
he Je =83,33 em4, yy =1,460-10% sec“, : 
pees oes I 33,180 eta) Way 1,314: 10% seen, : 
Be SN gs Mi, SOO1T0 ont ke. i 


‘ 


aba Micher den vAnsite (24) dreigliedrig (n =6) und setzen auBerdem 4M, =} Muipp fest. 


Wir neh 3 4) und | 2 
Dann folgen durch einfache Zahlenrechnung aus (28) die Eigenfrequenzen 
Be soi. tg = 1,263 10? seg, xe, 5.78810" see By | 
So 6, 61,920 -108 see>1, p= 3645-108 sect ; 
y= 1.294108 sec“, wg =9,639- 10? sect . 
ahettecnn vorhandenen Arbeiten on S. Woinowsky-Krieger, Ing.-Arch. 13 (1942), 
;. Hannover 1947, die sich tiber die Mathieusche Gleichung hinausgehend mit ge- 
beschiftigen, haben beide diesen Sachverhalt nicht klargestellt. . 
Krieger, Ing.-Arch. 13 (1942),S.90.0 Remon ari 
mich. Herr F. Weidenhammer wirksam unterstiitzt. 
re von. der Verdrehungsfestigkeit. | Berlin 1921.) “ 
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und aus (27) die Transformationsformeln 
w,= 0,9994v, + 0,03800 v2, 
Ww, = —0,03800 v,+ 0,9994 v,, 
w,= 90,9992 v, +.0,03982 v,, 
w, = —0,03982 v,+ 0,9992 v,, 
w= 0,9990 v, + 0,04406 vg, 
We =—0,04406 v;+ 0,9990 vg. 


Wir wollen ferner — und das ist der wesentlichste Punkt der Zahlenannahmen — der weiteren 
~Rechnung den Wert 


Wp = +H =2,055- 108 sec~1 
zugrunde legen, also den ersten Instabilitatsbereich zweiter Art untersuchen. Da dani @ stets | 


in der Nahe von @, bleiben wird, so kann man fir Ay. eine Naherung einfiihren: es ist 


Gon ee also 1 — (5) * l fair k=3. 


Ferner gilt auf Grund von (30) 


Also kommt 


=< mitvder Abkirzung » 


1 
B(a@) = 3 
ree) 
I iy 
Pa 7 
Um die Genauigkeit dieser Naherung zu iberprifen, nehmen wir 1—(2-) wl erst von k=5: ff 
an. Dann miissen wir 
: Sy ma ie ARE Se teenie i. nee 3x aoe 
a k on ers l 


setzen und erhalten eek kar ae als Faktor des rechts in (92) weggelassenen | : 
Gliedes sin 32x/l far w= =, den Wert 0,008. Die Kleinheit dieser Zahl dirfte die Vernach- 
lassigung des Gliedes sin 3 xl und aller héheren Glieder rechtfertigen. rae 
Aus (92) und (29) ergibt sich nun durch eine elementare Integration der mit € cos wt multi- 
plizierte Anteil von Ars Zu ; 


d ee Tt 26M rage = E 
; u ; Ae a ey here COs ot (w, 1 We a 9» Ws W, ale 25 Ws We) + or \ 
. 16feM tenes al Moral 7 1 1 9 ine 
| pe eet Cad ey Wy We — 35 5 105105 ++ rk 6 
a iehe et ie 7 Ws , Wy, — 79 9 Ws We — Bas alge WE Wo aie a a Fea oa ta) = 


In ‘diesen gd sind die oT cuit staan tforietoceln einzusetzen. Dabei. ellen wir noch, wi 
SEN bestimmten Zahlwerten zu gelangen, festsetzen, daB das Moment M,, dessen GréBe ja bisher 
_ ganz frei blieb, gleich dem Kippmoment Myipp ist, so daB also ¢ das Verhaltnis der Amplitude 
des tatsachlich auftretenden pulsierenden Moments zum statischen Kippmoment darstellt. 
Dann folgen aus waa die Koeffizienten Fem von (32). Sie gehen in die Koeffizienten | GY von 
(58) iiber, wenn man sie mit. @? dividiert und «=p setzt. Ferner gewinnt man phtepreevena: 
_ nach (55) und (57) die Taylorkoeffizienten Gy, Ge durch Differentiieren von B ( (w), und. ‘SeblieBlich — 

i folgen die Zahlen qy nach (56) durch Dividieres der Cp mit Wy. So ergibt sich die folgende Tabelle 
der Koeffizienten det Differentialgleichungen (58), soweit wir sie fir die Rechnung brauchen. | 


_ Dabei ist von vornherein Tier eases daB die Entwicklungen bis zur Potenz 8 getrieben werden 
oe sollen. : 
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Tabelle der Gy. 
Wegen c! = GO) braucht die Tabelle nur zur Halfte ausgefiillt zu werden. 


a i ee 
Pe 1 2 3 | 4 5 6 

1 0,5971- 10-2 

2 | —7,840-10-2 | —0,5971-10-2 

3 | —1,108-10-2 | 26,10-10-2 0,3253- 10-2 

4 —-2,861-10-2 | 1,150-10-2 | —4,075-10-2 | —0,3253-10-2 | 

5 | —0,4121-10-2 | 10,36-10-2 | —0,4809-10-7 | 8,623-10-2 | 0,9554-10-2 

6 0,3969-10-2 | 0,4725-10-2 | 3,091-10- | 0,5043-10-1 |10,81-10-2 | —0,9554-10-2 


Tabelle der Gi). 


ns 1 ; 62— —1,042 
~~~ (d0,2244-10-2 | 2947-10 gq, =0,06148 
2 2,947-10- | 0,2244-10-2 d= 0,9386 
3 0,2255-10- |—5,310-10-2 43=0,6297 
4 |—0,5821-10-2 |—0,2341-10-2 d4=2,816 
5 0,08385-10-2 |—2,107-10-2 Gar 
6 |—0,08076-10-1 |—0,09615-10-? _q,=4,690 


Mit diesen Werten kann man sogleich an die Aufstellung der Hauptgleichungen gehen. 
Beachtet man, dai} die ausgezeichneten Indizes jetzt k =1, r=2, p=1 sind und da (60) q,+q,=1 
lautet, so hat man nach (63) fir j =1 die erste Hauptgleichung K$), =0 in der Form (69) 


ah 1 1 
2 qo q@) = = Gf} cl) E ae cs] 


und die zweite Hauptgleichung KY) _,=0 in der Form (72) 
ee i ae 1 


—2 q, (G) — g) c) = 5 Ef, 


far j=2 die erste Hauptgleichung K®). =0 in der Form (70) 


‘ 6 
=e) = = 1 gq 
2 qn) = — G0 +5 D1 Oh ema, -1 + emda) + “gan 


und die zweite Hauptgleichung Ke) ;=0 in der Form (73) 
2 gy (G2) — 42) = — (G@— gq)? 4.2 gq (G2 —q®) o) + 


6 
Te 2 Pal cy, (oO), at: oho) rig Pate GY) ; 


AJarin ist 12) fire ka) 
| is Sa USE eae Breen e. 
man (at =n | malt GO), far 4 =0, —2 


und speziell c() _, =o, of, = 0. 

Da wir i j=3 gehen wollen, miissen die Hauptgleichungen fir j =1, 2 so ausfihrlich aus- 
gerechnet werden, dafB die g) und @‘) in ihrer Abhangigkeit von c und c( erscheinen. Durch 
Einsetzen der Tabellenwerte und Isolieren der q'/), q\/) erhalt man | 
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q®) = — 2,088-10-? &® , 
Tx] i 
qiescs | 2,088 e 4 81,90 =] -10-, 
q@ = — [0,1872 c( + 2,088 c@) + 4,382] -10-2, 
q? = | a + 00) (SP — 2,088) — 0,1872 0? — 1,749] -10-*, 
; 


0)? (0)2 


(entsprechend den Gleichungen (69), (78), (70), (79) in der am SchluB8 von Ziff. 9 angcRebeneny 
Gestalt). Daraus folgt Vigoh Anwendung der Vorschriften von Ziff. 10 


; om | — 43,908, a fe — 3,908, 
ef) = 2,426 (aus 22> = 0), = 5,163 (aus gf = qf). 
qi) = +,0,1632 , is ae 
gq?) == 0,1253 ; qo = —0,1231. 


Zur Berechnung der Glieder' j=3 miissen zunachst die Hauptgleichungen K),=0 und 
Ki) _ 7 = 0- nach(6353) serene: hingeschrieben werden, ferner die Ausdricke fiir die neu | 
darin auftretenden Koeffizienten Wee (f=1, 2,...6; 7=0, +1, —2). Setzt man darin fiir die | 
eV), g%), q@) die eben fir j =0, 1, 2 wigegebenou Zahlwerte ein, so kommt als Ergebnis ~ 


g®) = £0,03265,  J(q{) =— 0,09918 . 7S | 


Damit haben wir also die beiden maBgebenden Ausdriicke 
qo = 1¥ 0,1632 ¢ -— 0,1253 2% + 0,0327 63 + 0(e4) 
als Grew Laven des Instabilitatsbereichs und 


Sey 0,08162 ¢ + 0,09918 * + 0 (¢4) 
B(Os)an = = (24) x 1 — 0,1253 2 + 0(e%) 


ae Realteil des charakteristischen Exponenten auf der Mittellinie des Bereichs gefunden. Beide |] 
Ausdricke kann: man wohl unbedenklich bis zu etwa ¢=0,3 benutzen; man bekommt fiir ¢ =0,3- | 
ee 3 
qc = 1 0,04896 — 0,01128 + 0,00088 , 


0,02449 + 0,00268 
(On), = Sa, DLT 


0,02748 , i at 


_ woraus die praktisch befriedigende Konvergenz der | 
Reihen ersichtlich ist. Die Glieder mit ¢? bilden nur 
noch unwesentliche Korrekturen. In Abb. 5 ist der | 
so berechnete Instabilitatsbereich in der (q, )- ‘baw. 
der (w, ¢)-Ebene zeichnerisch dargestellt und durch. 
Schvatlue gekennzeichnet. Der Realteil des charak- 

_Aeristischen Exponenten (01), =0,02748 ist. kenn- 
zeichnend fiir die starkste Schwingungsaufschauke- 
lung. Die Zahl besagt, da®B sich die Amplitude der 
Zusatzverschiebung wahrend jeder V ollschwingung 
der Grundbewegung mit ¢?%- 0.02748 — 1,189. multi- 
‘pliziert. Nach 100 Schwingungen ist. die Amplitude 
schon auf das (1,189) —3,16-107-fache gewachsen. 
Da die Grundschwingung: eine Schwingungszahl von 
etwa 330 Hz hat, wird die Zerstérung des Stabes 
also rasch vor sich gehen, was bei der. erheblichen | 

Abb. 5. Der im Beispiel ES doheote Instabilitatsbereich Grobe des pulsierenden duBeren Moments ae | | 

zweiter Art, las Mae verstandlich ist. 


€. Grenzkurven 


e & 
= 


SchlieBlich ist noch das Konvergenzverhalten unserer Eosune ‘fiir verschiedéne n von intort 
_ esse. oe haben deshalb, da die gs, qi? von n unabhangig und die Glieder eq 3), 6 _ Jai 


i 


& 
“Al 
4 
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schon recht klein sind, den Koeffizienten gq noch fiir einen eingliedrigen und einen zweigliedrigen 
Ansatz (24) (n=2 und 4) berechnet. Das Ergebnis lautet 


—0,1071 fir n=2, 
q® = \ —0,1267 fir n=4, 
—0,1253 fir n=6 


§ = 4 te ® 
und zeigt, da8 unser Berechnungsverfahren fir n=6 und wohl auch schon fir n=4 einen hin- 


reichenden Grad von Zuverlassigkeit besitzen dirfte, wahrend der eingliedrige Ansatz (n =2) 
offenbar noch nicht ausreicht. 


13. Zusammenfassende Schlufbemerkungen. Hine hier als erste Mitteilung bezeichnete Arbeit 
des Verfassers entwickelte aus dem Hamiltonschen Prinzip der Elastizitatstheorie endlicher 
Verschiebungen durch Zusammenfassung der in den Zusatzverschiebungen quadratischen Teile 
eine Variationsgleichung, die das Stabilitatsproblem fiir die elastische Bewegung mathematisch 
formuliert. Die vorliegende Arbeit sucht das Problem zu lésen, indem sie einen Ritzansatz mit 
zeitabhangigen Koeffizienten in die Variationsgleichung einfihrt und die Eulerschen Differential- 
gleichungen fiir die Koeffizienten bildet. Da als Grundbewegung des elastischen Korpers eine 
harmonische erzwungene Schwingung angenommen wird, erscheint auf diese Weise ein lineares 
Differentialgleichungssystem mit periodischen Koeffizienten. Seine Integration gelingt in dem 
_erforderlichen Umfang durch eine doppelte Reihenentwicklung der Lésungen nach trigono- 

metrischen Funktionen und, dem Grundgedanken der Stérungsrechnung foleena. nach Potenzen 
eines kleinen Parameters ¢, der im wesentlichen die Amplitude der erregenden Kraft angibt. 
Aus dem Mechanismus der Reihenentwicklung heraus lassen sich in einer Koordinatenebene, 
auf deren Achsen die erregende Frequenz w und der Parameter ¢ aufgetragen sind, die Gebiete 
-abgrenzen, in denen die Grundbewegung stabil oder instabil ist, sowie die fiir die Instabilitat 
maBgebenden charakteristischen Exponenten berechnen. Nahere Einzelheiten sind in Ziff. 11 
| zusammengestellt. Kin Beispiel (Ziff. 4 und 12) zeigt, daB die in allgemeiner Form durchgefihrte 
_ Theorie in konkreten Einzelfallen numerisch brauchbar sein kann. Das Beispiel behandelt einen 
der hier neu gefundenen Instabilitatsbereiche zweiter Art, die mechanisch dadurch gekenn- 
_geichnet sind, daB die Aufschaukelung der Schwingung durch Zusammentreffen der Erreger- 
frequenz mit einem ganzzabligen Bruchteil der Schwebungsfrequenz zweier Eigenschwingungen 
des elastischen Korpers ermoglicht wird. Im uibrigen muBte die vorliegende Untersuchung eine 
_ Reihe von Fragen offen lassen: Das Problem der Konvergenz der lediglich formal eingefihrten 
_ Reihenentwicklungen ist noch ganz unerledigt und schlieBt sicherlich schwierige mathematische 
_ Fragen in sich. Ferner muSten die Sonderfalle, in ‘denen gewisse Gliednenner verschwinden, bei- 
seite bleiben und werden Spezialuntersuchungen notwendig machen. Weiter konnte die Damp- 
fung, die beziiglich der Ubereinstimmung zwischen Theorie und Wirklichkeit sicherlich eine sehr 
erhebliche Rolle spielt, lediglich mit einer kurzen Bemerkung qualitativer Art gestreift werden . 
(Ziff. 11). Und schlieBlich bleibt durch Behandlung méglichst vieler konkreter Einzelprobleme 
zu prifen, wie weit die praktische Leistungsfahigkeit der vorliegenden, Theorie bei der Durch- 
forschung dieses jungen Seitengebietes, der Elastizitatstheorie reicht, das wir mit dem Titel 
Stabilitat erzwungener Solwimeyaeen bezeichnen und das unzweifelhaft noch weit 1 in den An- 
fangen ; steckt. es 


pact 5 (Eingegangen am 15, Dezember 1948.) - 
Oy “Ansehrift des Verfassers: Professor Dr. Eberhard Mettler, (20 b) Clausthal- Zellerfeld, Ae aa 
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Uber ein nichtlineares Elastizitatsgesetz. 
@ Von H. Kanuderer. 


1, Einleitung. In der Elastizitatstheorie der kleinen Verformungen pflegt man von der An- | 


nahme auszugehen,; daf zwischen den Spannungen und den Verzerrungen ein linearer Zusammen- 
hang besteht. Es ist indessen seit langem bekannt, daB far eine groBe Anzahl von Stoffen diese 
Hookesche Annahme auch schon in dem Bereich der von erster Ordnung kleinen Verzerrungen 
nicht mehr gerechtfertigt ist, falls man, wie tblich, unter von erster Ordnung kleinen Verzer- 
rungen solche versteht, bei welchen die Quadrate der raumlichen Ableitungen der Verschiebungen 
3 sAeabhy den ersten: Potenzen dieser GréBen vernachlassigbar klein sind. Diese Abweichung 
vom Hookeschen Gesetz ist unter den in der Technik gebrauchlichen Stoffen besonders deutlich 
bei GuBeisen, zahlreichen Nichteisenmetallen, Steinen, Zement, Beton und den meisten organi- 
schen Werkstoffen; in geringem Ausma8 ist sie auch bei einigen Eisen- und Stahlsorten fest- 
stellbar. 

Schon mehrfach ist daher versucht onder durch eine geeignete Annahme eines nichtlinearen 
Zusammenhangs zwischen Spannungen und Verzerrungen dem elastischen Verhalten dieser 
Stoffe gerecht zu werden. Dies geschieht z. B. in dem Ansatz von Bach und Schiile1, die dureh 
eine letey ne der Form 

€= a0", 


wobei ~ und m experimentell zu bestimmende Materialkonstanten sind, den Zusammenhang | 


Zug- oder Druckversuch beschreiben?. 


Ein anderer, von Griibler® beniitzter Ansatz von der Form 


E 
ee | Oe Ad Bes 
mit den Materialkonstanten A und B, der von diesem auch fir allgemeinere Spannungszustande 


stand iibergefithrt wird. Der gleiche Nachteil haftet einem von Hencky® insbesondere fur groBe | 
. Verzerrungen vorgeschlagenen Gesetz an, was dieser auch selbst feststellt. 


SchlieBlich sind zwei Aufsaitze von Philippidis und Weber® zu erwabnen, die wihrond der 


Abfassung dieser Arbeit erschienen sind. Beide handeln von einem Elastizitatsgesetz, das in | 
seinem Aufbau grundsatzlich mit dem hier neu herzuleitenden ubereinstimmt. 


ne 
Piclippidis geht zur Ableitung dieses Gesetzes vom Hookeschen Gesetz fiir den dreidimensal 
sionalen | Ppenmunseppetand in peisocdier Form aus: aa 
Re: ue & 3 RK; O1—Oo ig acs ame p V (665) ae Gare ae ' 

(Se See Fei Gere etek Se a V (€;—€)? +. (€.—€p)? =F (e3—€6)? ae ; 


und die Hauptdehnungen, o, und ¢, deren Mittelwerte bedeuten. Er verallgemeinert dieses | 


pas Gesetz, indem er Bae Konstanten K und G durch beliebige, nur von Ey baw. von | 


BE pack Vel ee)* Hher—en)* + essed 1, | : | 


ue a ze B.C. Back und ve Baumann, Paes und Festigkeit, 9. Aufl. Berlin 1924. mage! | 
alist ine Zusammenstellung der bis 1897 versuchten Nah tze fii - t | 
_R. Mehmke, Z, Math. Phys. 42 (1897) 8. 327. pak daviaapei rin paae Shai! | 
8 M. Griibler, Z. VDI. 44 (1900) $. 1157. 3 | “Fl 
SEL: Schlechtweg, Z. angew. Math. Mech. 14 (1934) 5.5. ; a: Bee. || 
: ay ay Z. techn. Physik 9 (1928) S. 215. ‘ Oeste 

P ilippidis, Z. angew. Math. Mech. 25/27 (1947) S. 31. — C. Weber, Z. Mat 
8.189, sowie Z. angew. Math. Mech, 29 (98 s ey eee a er, Z. angew. Math, Mech. 28 (1948) | 


vans 


Bi 
p 
ee 
& 


_awischen dem Betrag o der Spannung und dem Betrag ¢ der Dehnung heim eindimensionalen jf 


~ wobei- K den Kompressions-, G den Schubmodul, o; und ¢; (i=1, 2,3) die Hauptspannungen | | 


erweitert wurde, wird, wie Schlechtweg* bemerkt hat, einigen wesentlichen Eigenschaften des |f 
Materials (z. B. Querdehnung) nicht gerecht. Er hat daher fiir spréde Stoffe den Griiblerschen ll 
Ansatz verbessert. Das Elastizitatsgesetz von Schlechtweg hat jedoch, wie wir spater nachweisen |} 
werden, den Nachteil, daB die sich bei ihm ergebende Formanderungsarbeit auch im Gebiet |] 
beliebig kleiner Verzerrungen vom Weg abhangt, auf dem der Kérper in den verzerrten Zu- | 


ard rye ares 
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- abhangigé’ Funktionen ersetzt, woraus er schlieBt, daB zwischen diesem fiir das elastische Gebiet 
aufgestellten Gesetz und der im Verfestigungsgebiet experimentell nachgewiesenen’ Beziehung 


V (=0,)2 )? + (02-0)? + (d3—0,)? = fLV ( (€1=£ )? +(€2—€p)? +(é3—&)? | ’ 
in der f eine beliebige Funktion bedeutet, kein Unterschied besteht. : 


Weber geht im wesentlichen von der Annahme aus, da8 sich (wie beim Hookeschen Gesetz) 
die Formanderungsarbeit in zwei Summanden zerlegen lat, von denen der eine nur von &, oe 


der andere nur von (€,—é,)? +(é— £9)” +(€3—€ )? abhangt. Das Pins tiga peereete erhalt hier- 
durch wieder die gleiche Form wie bei Philippidis. 


Beiden Herleitungen dieses Gesetzes ist somit gemeinsam, daf eine wallurlich, aber geeignet 
‘gewahlte Formulierung des Hookeschen Gesetzes herausgegriffen und in mathematisch nahe- 
liegender Weise verallgemeinert wird. Diesem mehr vom mathematisch formalen Standpunkt 
aus geleiteten Vorgehen soll hier ein Verfahren zur Ableitung des nichtlinearen Elastizitats- 

_ gesetzes gegentibergestellt werden, das die physikalische, experimentell gesicherte Seite des 
_ Problems in den Vordergrund riickt. Das Gesetz soll dabei zugleich auf eine mathematische Form 
_gebracht werden, die seine Handhabung sowohl bei der experimentellen Auswertung als auch 


_ bei der Anwendung auf die im folgenden behandelten elastomechanischen Probleme méglichst 
pomiach und ubersichtlich macht. 


5 


2. Aufgabenstellung und Uhersicht. Es soll also hier der Versuch gemacht’ werden, unter 
Pe chasnkune auf von erster Ordnung kleine Verzerrungen ein nichtlineares Elastizitatsgesetz 
_ herzuleiten, das auf beliebige, mit den Bedingungen der Statik vertragliche Spannungszustande 
anwendbar ist. Durch vier physikalische Forderungen an das elastische Verhalten, die besonders 
- die in der Einleitung erwahnten Mangel friherer Angize ausschlieBen, wird die Form des Ge- 
_setzes im_wesentlichen bis auf zwei willkirliche Funktionen des: Spannungs- oder des Ver- 
_zerrungstensors festgelegt, die das Material kennzeichnen und an Hand geeigneter Grundversuche — 
mu u bestimmen sind. : (eae 


Nach einem kurzen, wegen der Definition der beniitzten Begriffe notwendigen Uberblick 
ide die Invarianten und die Zerlegungsformeln des_ Spannungs- und des Verzerrungstensors 
Ziff. 3 und 4) werden die vier Forderungen an das elastische Verhalten, von denen sich drei | 
auf die Existenz und die Gestalt eines Ausdrucks fir die Formanderungsarbeit beziehen, formu- 
liert (Ziff..5). Hieraus wird das Elastizitatsgesetz (Ziff. 6) sowie der endgiltige Ausdruck fir die 
Formanderungsarbeit (Ziff. 7) hergeleitet. Das Elastizitatsgesetz wird zuerst auf einfache (groBten- 
ils gleichférmige) Spannungszustande angewandt, die eine exakte Lésung der elastischen 
Dias v: zulassen Boe 8). Sodann wird an drei Beispielen gezeigt, ™ wie man durch die 


ae. duveh das hice entwickelte Gabe: Beeson. ieBe, Fico schlieBt sich. eine ede 
erst oe mit den von Bach, Cans und Schlechtweg formulierten Gesetzen (2ift. 10). ‘Mit 


Biegung, Be sich hier nicht mehr iewleutne behandeln. labe (Ziff. 11), und a dies Torsion 

smatischer Stabe (Ziff. 12), soll ein erster Schritt getan werden auf dem Weg zur Entwicklung 
-exakten nichtlinearen Elastizitatstheorie kleiner Verformungen. i bach hevend fate a a 
W rden Hinweise ‘gegeben, in welchen Richtungen_ ae Theorie auszubauen ist. ce ie 


Meg ? 4 ‘i 


; Ties Spocnasueciiand.. Den cUpecemugiinven ee ein rechtshandiges kartesisches (x, ¥2)= 
em zugrunde gelegt werden. Wir fiihren die Normalspannungen Ox, Oy, Oz am Elementar- 
dx: nay dzinw licher Weise ein (als Zugspannungen mit positivem Zahlwert). Sind die : 
H uptspannungen, so seien sie auch mit 01, G2, 03 bezeichnet. Die Schub- — 
best Texs wobei. schon zum. Ausdruck gebracht ist, dab der Satz von der 
n. Schubspannungen gelten soll. Der Spannungerusrand JaBt sich dann 


/ 
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und voneinander unabhangige Ausdricke, die wir in der Form ® 
S, =0z+0y+0; = 0,+0,+03 5 
Sp = Ox0y+0y0z+020x—U gy —Vyz— Vex = 9492 +0203 +0 3% » 
| Ox Try Tx 
S3= | Txy Oy  Tyz | = 01 Op G3 
Tie Gyan Os 


schreiben wollen. Den Zustand T teilen wir in bekannter Weise auf! in den hydrostatischen 
Spannungszustand, der durch den Tensor 


Ty(%, 9.2) = 09° E (2) 
gekennzeichnet ist, wobei 
1 1 1 

aan (Ox+Oy+0s) = 3 (01+02+03) = 3 Ss, (3) 
den Mittelwert der drei Hauptspannungen und E den Einheitstensor bedeutet, und den durch 
Ox— O09 Txy Tax ; 
T" (x,y, 2) = Uny i Oyr Og Tyz (4) 
Tux Tyz 0z—O9 


bestimmten Spannungszustand. Da a, eine von der Wahl des Koordinatensystems unabhangige 
GréBe ist, gilt dies auch vom Tensor Ty (x, y, 2) und somit auch von T’(x, y, z) als zweitem — 
Summanden der Tensorgleichung : Be A 
eee | (5) 
Die analog den Invarianten S,, S,,.S3 aus dem Tensor T’ gebildeten Ausdriicke sind daher eben- 
falls koordinateninvariant. Man erhalt sie mit den Abkirzungen : 


vera 


On = Ox — Cos Oy = Oy~Gys 6, =6,—= 
in der Form 


eat ee | . (6) 
; : 

rn Sas = Oyt Oy o,+oa.o, Tay Tye une | 
= x a 2 of ‘)- yaya —T2, : eats 0 


oe Txy Tix 1 s E 

, , ; 

| S3=|tey GO, ty |= Ss—55S,Sy+# St. (Bye 
7 - 

Z ; : Tax Uys 0, | , ye * 


eden Hauptspannungen von T aegis wird 


S=- ; +(e (01 —9)* +(62—Go)” + (G3 —Gp) ake 


Die Invariante S; ist somit stets negativ; ihr Betrag ist das Anderthalbfache der mittleren 
quadratischen Abweichung der drei Hauptspannungen von ihrem Mittelwert oy. Sie gibt also 
ein skalares MaB dafiir, wie stark der Zustand T vom hydrostatischen Zustand T, abweicht. | 


Zwischen den Spannungen | und den Komponenten X, Y,Z einer auf die Volumeneinheit — 
bezogenen auBeren Kraft § bestehen | die pe eanton Gleichgewichtgbedingungen 


ee : ~~ sce : at ots =a oe X= ay 0, ; “3 | : o i 
i \y x . x Cea ie 

‘, \ ar - . 4 aes ee ‘ 

oe ey Bey yo, | Oi fare 
te Ota jin Ne ; = are t | 

ta at ay -—+ Z= { 
die sich, unter Verwendung des Begrilis der ipl Saat als Vektorgleichung in der F orm 
div T+R= 0 


. I = ¢ : 
1 Wir werden kiinftig aah durch einen Tensor T bestimmten Spannungs 
gustand k \ 32 oe 
bezeichnen und uns spiter beim Verzerrungstensor sinngemaB are “ae koe ie 
2 Siehe z. B. R. Gans, Vektoranalysis. Leipzig und Berlin 1929, S. 76. 
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2 echiiben laser In den Komponenten von Ty und ’ ausgedriickt, lauten die Gleichgewichts- 

_bedingungen 
Se £35 ee 44 a a 
- mr , 3 P) x EIN + f »¢ Fama’ 0 > 


“0 be 005 


0x 


— 


und zyklisch weiter, oder i in Vektorforn rieeriany | : 
div T’ + grad Fo t= ant 


Der Pistand. T’ ist, fir sich genommen, mit der:Volumenkraft @ zusammen demnach nur dann 


ein Gleichgewichtszustand fir den ganzen Kérper, falls grad ¢)=0, also Oo=konst ist. ° 
Xe 


3 a 4. Der Verzerrungszustand. Zur Beschreibung des Verzerrungszustandes gehen wir aus von 
a dem raumlichen Feld B(x yz z) der Verschieburgsvektoren, deren Komponenten im (x, y,z)-System 


h.. mit u, v, w bezeichnet seien. Da wir uns auf die Theorie der von erster Ordnung kleinen Ver- 
_ gerrungen beschranken wollen, fahren wir als VerzerrungsgréBen die drei Dehnungen 


tanh Cees en Wen Fie) eG eis NOR Br a es dw 2 
ct a 6 fs see ate | = — so Wy! | : 10 
me dx Cae oy” ay ie ay) 


D(x. ¥. 3) = 


f 


Stoic. Bagiigen ex, ey, e Heer aeene so. seien sie auch mit. es Eos &3 be- _ 
ich) wie zum. Spannungetensor Z ae zum Versérrungstensor D drei Invarianten, 


=e exbey ter = = -erheaches 1 eee 
= exert epee enes— 5 (Yay Vs +h): = aebeer tae 


z Pe 


u et Volumenelements bedeutet, cu 


+ Siri b 


Sey 
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die Invarianten von D’, namlich 


Jose te+e=0, ae (16) 
aS | Soa ee te tee—lwi ty tv) | 
ue =—Fle be ey aye 2) — 5 Way + Von + Pia) (17) 
3 = J,— di, i i 

, 1 1 
yee) Eg gy Pru 7 Y= 
, 1 
ss 3 Wxy éy = i er w j ies Jide 5a = Ripe, (18) 
he i pya é. 


i 


In den Hauptdehnungen ausgedrickt, wird 
ae Jo= — 3-5 [e—e0)? + (e260)? + (65 &0)"I - : | 


Die Invariante J; ist somit stets negativ; ihr Betrag gibt, ahnlich wie bei Sj, ein Ma8 dafir, 
wie stark der Zustand D vom gestaltstreuen Verzerrungszustand D) abweicht. 


Da sich die sechs VerzerrungsgréBen von D in den Ableitungen der drei Verschiebungs- 
komponenten u, v, w ausdricken lassen, miissen zwischen ihnen die sechs (selbst nicht von- 
einander yuabhaneigen) Vertraglichkeitsbedingungen} 


dee | Wey Way Mas 0 dye 8 ey Oey. 1oys 
Pe uae aye Toor Pete cal ee ee (124 


; und je zyklisch weiter, erfiillt sein. Aus dem Bestehen dieser Gleichungen folgt, daB die ent- 
eg sprechenden zum Tensor D’ gehorigen Vertraglichkeitsbedingungen nur dann befriedigt werden _ 
nk kénnen, wenn samtliche zweiten Ableitungen von é) verschwinden, also é eine lineare Funktion |} 
yon x, y und z oder eine Konstante ist. Dies bedeutet aber, daB die am Volumenelement vor- |} 
genommene Aufteilung des Verzerrungszustandes D in einen gestalts- und einen volumentreuen _ |f 
-Anteil nur dann auch am ganzen Kéorper moglich ist (unter Wahrung des Zusammenhangs 


der einzelnen Teile), wenn é) der genannten Bedingung entspricht. 


t 


5. Forderungen an das elastische Verhalten. Um das Elastizitatsgesetz, d.h. die Form des 
Zusammenhangs zwischen dem Spannungszustand T und dem Verzerrungszustand D zu er- 
mitteln, stellen wir jetzt an das Verhalten des Materials vier allgemeine, sich aus physikalischen’ 
: Hypothesen ergebende F orderungen. Die ersten drei davon werden sich dadurch in mathemati- 
cher Form ausdricken lassen, daB sich fir die skalare Funktion A, welche die am Volumenele- 

_ ment geleistete, auf die Raumeinheit bezogene mechanische Arbeit bedeuten mége, die bei der — 
Verformung aufzuwenden ist, eine ganz bestimmte Form der Abhangigkeit von den Verzerrungs- — 
gréBen ergibt. Aus dieser Form von A werden 1 wir dann in der nachsten Ziffer das Elastizitats- _ 

- gesetz herleiten. ° Le eerie lees rae Le ‘ : 


i 1. Fordérung. Die’ Verzerrungsarbeit A soll eine eindeutige Bunkeos: ee oe 
_ gustands des Elementes, also der sechs den Tensor D bestimmenden ‘Verzerrungsgrofen s sein. 
_ Sie soll ve verschwinden, wenn D gleich Null ist. - ete "NN 


Um nun aber durch’ diese Forderung das Elastizitatsgeset2 nicht auf Hieseeeal Store. ya 
a uate ZU miissen,. die vollkommen elastisch im Sinn der energetischen Definition dieses Be- re 
_griffes sind 2, 2, kann man etwa_ einrdumen, da fiir den Weg, auf dem das Volumenelement vom 
i vememangetreion Ausgangszustand in den Zustand D ibergefihrt wird, gewisse Vorschriften 
_ gemacht werden. Es wird sich hier hauptsdchlich darum handeln, daB man nur solche Ver- 
formungswege zulaBt, bei welchen die Betrage der Verzerrungskomponenten bis zum. Ubergang 

in den Endzustand nicht abnehmen. Ein Spezialfall hievon kann etwa dadurch definiert werde ie 
-daB der Verzerrungszustand zur Zeit t durch den Tensor #(t) + D bestimmt wird, wobei o see 
eine beliebige, zwischen dem Anfangszeitpunkt und dem Endzeitpunkt , Verzerrungevorganges - 
“monoton von Null auf Eins arrore moras Bunko: Lode Gapee eN 


1 


1 Siehe z. B. C. B. Biezeno und R. Crameols Technische Dyce: Berit 1939, s. 23, 


a - 2 Nach dem Vorschlag von G. Sachs, KongreBbuch Zirich | ges alee onaren, Verbandes fiir Material. 
lia eared 1932, S, aae FAO ta ; 
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2 2. Forderung. Der Stoff soll homogen und isotrop sein, also hinsichtlich seines elastischen 
Verhaltens keine ausgezeichneten Stellen oder Richtungen aufweisen. Es ist daher auch von der 


! Verzerrungsarbeit A zu fordern, da® sie allein von den koordinatenunabhangigen Invarianten 
Jy, Jz und J, des Verzerrungstensors abhangt?. 


3. Forderung. Die folgende Forderung wird nahegelegt durch die v. Misessche Plastizitats- 
hypothese?, die durch zahlreiche Versuche von Ro und Eichinger® bestatigt wurde. Nach dieser 
wird die FlieBgrenze eines Stoffes allein durch den volumentreuen Anteil T’ des Spannungs- 
tensors bestimmt. Es erscheint daher die Annahme berechtigt, da auch schon vor dem Eintritt 
in das plastische Gebiet die zusatzliche Verzerrung D,, die ein Volumenelement infolge eines 
Spannungszustandes T’ erleidet, unabhangig ist von einer etwaigen Verzerrung D,, die das Ele- a: 
ment vorher durch einen hydrostatischen Spannungszustand Ty erfahren hat. Die Verformung — ek 
infolge von T) kann aber wegen der Isotropieforderung nur in einer gestaltstreuen Verzerrung __ hanes 
bestehen. Denkt man sich nun den ‘Spannungszustand T in der Weise aufgebracht, daB zuerst | 
_  T, und sodann zusitzlich T’ am Element angreift, so wird die entstehende Verzerrung D sich 
_. zusammensetzen miussen aus einem Kugeltensor D,(T)), der sich in der Form cD, mit irgend- 
_ einer Konstanten ¢ durch den gestaltstreuen Anteil D, von D darstellen la®t, und dem Tensor 
D,(T’), der dann den Wert (1—c) D)+D’ annimmt. Da aber ‘ 


hese 
; (1—c) D) = —* D,(T,) 


_ eine Funktion von Ty, ist, darf dieses Glied in dem Ausdruck fiir den Tensor D,, der nur von 
T”’ abhangen soll, nicht auftreten, so daB c=1 und damit D,(T,)=D, und D,(T’)=D’, also 


D(T) = Do(To) + D'(T’) 


gesetzt werden muB. 33 : 

__ Diese Beziehung fiir den Verzerrungstensor muB jetzt.sinngema8 auch fir die Verzerrungs- 
arbeit gelten. Wir fordern also; daB sich auch A aufteilen 148t in einen nur von D, abhangigen 
- Anteil A, und einen allein durch D’ bestimmten Anteil A’. Beide Anteile miissen die zweite 
_ Forderung erfiillen, d. h. Ay darf nur von der (einzigen) Invarianten J, von D, abhangen, wahrend 
_ A’ eine Funktion von J; und J; sein darf. Man erhilt so die Zerlegung she poke 


a AU ee 


i) 


: oe Psy Ao(J;) +A’ (2-3 Jt Joya dy Jo+ 7 J) < | ts ny | Si a 


Ae ord erung. Durch Versuche ist erwiesen, daB das Verhalten der festen Kérper bei hin- i 
_reichend kleinen Verzerrungen beliebig wenig vom Hookeschen Gesetz abweicht*. Wir wollen | |— 
daher allgemein fordern, da das zu formulierende Elastizitatsgesetz fir jeden Verzerrungstensor — 


erster Ordnung nicht mehr vernachlassigt werden dirfen, => 


bei den in der Theorie der kleinen Verzerrungen iblichen Zahlwerten gegentiber * 


4a | 


erle tung der allgemeinen Form des Elastizitatsgesetzes fiir Kleine Verzerrungen. Wir _ 
1un, ausgehend von dem -Ansatz (20) fiir die Verzerrungsarbeit und unter Beachtung 
en Forderung, untersuchen, was fiir ein Zusammenhang zwischen den Komponenten 
ngs- und des Verzerrungstensors bestehen mul. — . ae 
erhalt man unter der Voraussetzung, daB der Verzerrungsvorgang isotherm — 
| rzerrungen die Komponenten des Spannungstensors als die partiellen 
‘erzerrungen ausgedriickten Verzerrungsarbeit A nach den entsprechen- 3 
s Ve errungstensors, wobei der Faktor + bei den Scherungskomponenten e be me 


ntare Mechanik. Leipzig und Berlin 1912, S. 571. 
. Géttingen, Math. phys. Klasse (1913), 5.582, RP 
che zur Klarung der Frage der Bruchgefahr. Ziirich 1926, sowie Dis-  — 
ngsanstalt in Ziirich, Nr. 28 (1928) und Nr. 34 (1929), 
4 (1906) S. 469, und F’. Kohlrausch und E, Griineisen, Sitzungsber, 
Klasse 1901, S. 1086. f Mesh . 


* f , rite hah ae j i 


oes 
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wegzulassen ist. Da wir auf Grund der dritten Forderung die Verzerrungsarbeit A am Volumen: 
element gema® den zwei voneinander unabhangigen Verzerrungsvorgangen Dp und D’ in die 
zwei Anteile A, und J’ aufteilen, so mu jener Zusammenhang auch fir die Komponenten der 
die Arbeitsanteile A,(J,) bzw. A’ (Jj, J3) leistenden Spannungstensoren Ty und T” je fir sich 
bestehen, d.h. es muB sein 


eT hah RCd Ay el Sea idee | (21) 
DORs G Reh Bod Suet oad ea ee 
fa (0 -At tise 0A’ ; ieee a i aes eee 5 Tees ’ ee hae (22) 
% de,’ , y 0 Ey’ 5 0 E,’ OW xy 0 Wyz 0 Pex 


(Der Faktor = ist in (21) deshalb einzufiigen, weil die drei Dehnungen fir die drei Koordinaten- 
richtungen samtliche mit dem gleichen Symbol ¢, bezeichnet werden und die Ableitung natirlich 
nur jeweils partiell nach einer dieser drei DehnungsgréBen zu bilden ist.) So werden wegen 


0 Ja’ Pane gh ahr ners 2 


1 2 
E VON Gg aan 
de, x? 0 x’ b amd A, Wyz 


und der hieraus durch zyklische Vertauschung entstehenden Gleichungen die Spannungen 


ee Gh Pena ea he (sg ON :) 
0, = ITs. 0 Ex! see Fp 0 Ex’ rR: Ex “fe oJ, y oz A Pyz 3 
aan 0A’ Bilan 0,.A’ Osis ae 0A’ / 3 A’ ( Le Ia. it 3) 
oy Files ey’ + 0 Je’ 0 ey’ 0 Jo’ éy 0 Js’ Ez Ex, 4 Wrx} > 
ene 0 A’ Ce Pe 0 A’ alg fxs “% 0A’ , 0A’ ( Pia T, 1 2 ) 

EOE PETER Oe Te Fe ae ee 


Da sowohl nach (6) “ey 
Og 4-0, + 6, = 0 
als auch nach (16) 
ete te =0 
ist, folgt wegen (17) aus der Summation der drei Gleichungen, daB 


0A’ / / 4 / y Lf 1 aA 3 
ae. a €, Ex te oaey tap 4 Ve + Yet vy) | = aay J =0 


wird. Dies bedeutet aber, daB A’ von J; unabhangig sein muB, da ja J; nicht identisch verschwin-_ 
den kann. Dann wird aber wegen Z . 


PA ones & 
/ : 4 \ \ 0 Wry. 2 Vay 
die Schubspannung : 
BAe BIg die 


IN aA Va oe EES 


und zyklisch weiter. Somit lassen sich die Beziehungen zwischen Spannungen und Dehnungen 
mit Hilfe der beiden, nunmehr die Form A,(J,) und A’(J;) annehmenden Ausdriicke fir Vo- | 
Jumen- und Gestaltanderungsarbeit je Volumeneinheit in den beiden folgenden Gleichungen | 
zusammenfassen Lge ¢ 

dA. eer yee 
ot We Oe Cie ad cee pcs | | 
Die beiden Funktionen dA,/dJ, und dA'/d J, bestimmen daher vollstandig das elastische | 
Verhalten des Kérpers. Ihr Verlauf hat allein der vierten Forderung zu geniigen, darf aber sonst — 
willkirlich sein. Es ist zweckmaBig, an Stelle von J, und J; als neue Veranderliche die in (14) : 
eingefihrte mittlere Dehnung ¢, und eine neue reelle Verzerrungsgrébe Wo deren Quadrat 


(Oi 


os 


2 


Basha te + ety tet 3) 4 


‘sei, zu beniitzen. Dann 1aBt sich dA,/dJ, auch als Funktion von ¢, und dA'ldJi als Funktion © | 
von p> darstellen. ; - 

. Nun lautet bekanntlich das Hookesche Gesetz, nach 0» bzw. T’ aufgelést, mit dem Kom- |} 
pressionsmodul K und dem Schubmodul G pais 


0,=3Ke, 3 Ley QDs (24), (25) att 


- é ; 
1 Siehe z. B. H. v, Helmholtz, Vorlesungen iiber theoretische Physik, Bd. II. Leipzig 1902, S. 77. 


DP te the MIS Ta SUS 
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Setzen wir also dA,/dJ, und dA’/dJ, mit den zunichst beliebigen Funktionen x(g,) und Y (y%) 
in der Form 


dA, dA’ 


Tye 3 Key: Heo), GH = —26- yh (26) 
an, so muf} offenbar, um der vierten Forderung zu gentigen, 

lim %(€) =lis4 lim y(y%) = 1 

E90 po 


sein. Wenn fiir x(é,) und y (yo) Potenzreihenentwicklungen méglich sind, werden diese somit 
- die Form, haben missen 


: %(&) = 1+ %,& + Hy & xe + +: » | (27) 
: 7 (po) =1 + yao + vawo + Yoyo +--+ > J 
und das Elastizitatsgesetz erhalt dann die Gestalt 
: Gp =3Kx(G)eq, TT’ = 2Cy(yi) D'. (28), (29) 


Wir bezeichnen die Funktion x(e¢)), mit welcher der Kompressionsmodul K in (24) zu multi- 
plizieren ist, um (28) zu erhalten, als ,,Dehnungsfunktion‘’ und entsprechend die Funktion 
y (yo) als ,,Scherungsfunktion“. Der Verlauf dieser beiden Funktionen bzw. deren Entwicklungs- 
koeffizienten x; und y; sind, zusammen mit.K und G selbst, durch geeignete Grundversuche 
experimentell zu bestimmen. Wie sich eine solche Bestimmung durchfihren laBt, wird in Ziff. 9 
gezeigt werden. 

Addiert man die mit dem Einheitstensor E multiplizierte Gleichung (28) fiir, O zu der Glei- 
chung (29) fiir T’, so folgt unter Beachtung von (2), (5) und (13) 


: T= T,+T' =3 Kx(&) Dy +2 Gy(y2) D’ . (30) 
- Fir die einzelnen Komponenten von T ergibt sich dann 


Ox = 3 Kx(Eq) & + 2 Gy(yi) (Ex—&>) » 


dy =3 Kx(é)& +2 Gy(y%) (Ey—£p) » ; (31) 
dz = 3 Kx(é) & +2 Cy(y) (e — £0) ; Ng 
Tay = CY (Yo) Pry Ts =GY (Yo) yz» Tee = OY (Wo) Pax - (32) 


Um die Gleichungen nach den Verzerrungen aufzulésen, miissen wir die Dehnungsfunktion 


_und die Scherungsfunktion durch Funktionen von SpannungsgréBen ersetzen. Hiezu bilden wir ~ 


in der Gleichung (29) beiderseits die zweite Tensorinvariante; dann erhalten wir 
8, = [2 Cy(yo)]? Ja 


oder, wenn ‘wir eine neue wens Spananieserobe T) durch 


3 B= Sa aly (ORF og + 204) thy + He te J (33) 
einféhren, wegen (23) 
ee ; = Cy (yi) Yo} eee (34) 
| Mit den dimensionslosen Se ey 
oe peter Pars Sopa ae. Z 1. (35) 
: ay ae a) G 
Jauten 5 somit die selielen gee (28) und (34) if ‘ 
SS ‘ Sq = (Eq) £0 > to = (Yo) (36) 


Wir Hilwen, ee mit Hite. zweier Funktionen von sy ee i, namlich der ,, Kompressions- 
“tes k(s,) und der RCL UP IARRUOR 8 (to) die Umkehrungen dieses funktionalen Zusammen- 


hangs in oer Form | | 
2 . eo = z= soh(s So) > Yo = b8(to) ; (37) 


2 die Potenareihen fiir k(so) und g(t?) die Gestalt. 


a 


Pe ahioes (So) = 1 + hy 89 + ikass + has + + as ce (38) 
a SNe y (to) = 1+ Bato + Bato + Beto + 


i Daw alee Cleichung auch im Ta Tieh nb berbicll hasiekt. haben Rok, Eiichinger er a, O.) sowie 


vehi Nee inate 3 eo? S. atl cercementell Daeheewiesen, 


ein. ‘Sind die Fuakaean #(E 9) und y(wo) in Potenzreihen nach (27) entwickelbar, so erhalten — 


aN 


ree 
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wobei nach den Regeln tber die Umkehrung von Potenzreihen 


ky =—%,,; hy = 2 ui —%e; ky = —(5 Hi —95 #1 Hq + %3) » 

keg = 14 x4 — 21 23 ty + 6 vty 2g +3 0g — qs (39) 

ks =—7(6 we—l2 xi wg t4xingt 4x, 1B — Hy Hy — Ugg + 7 %s) ; | 
=e, Bs=3¥a—Yar-- ee (40) | 


zu setzen ist. Die Gleichungen (39) und (40) behalten auch noch ihre Giltigkeit, wenn wir die ot 
Koeffizienten x; mit k; und y; mit gi vertauschen. Aus dem Elastizitatsgesetz (28), (29) wird so 


r 1 1 
: oa Mju se Dee 41), (42) 
) ; oder, werin man die mit dem Einheitstenor E Bagi cgape Gleichung zur zweiten addiert, 
ae x 
D=D,+D= 3K k(S9) DW en liar ce a g(t) T’. (43) 


Fur die einzelnen Komponenten von D folgt dann 
be = aa bls) 9 + se Bl) (2-0) » 
(oy = Ze Rlse) oy + eq Blt) (Gyo) » | (44) I 
ex 2 oe k.(89) 0 + se a(t) (Gz—Op)5 i nt é | 


: ‘ 1 By 

Yer = GUA Tay Y= Teles Ye — slit. (48) 

Die 15 Gleichungen (9) bis (11) mit (31), (32) oder mit (44), (45) stellen zusammen mit den 
Definitionsgleichungen (3), (14), (23), (33), (35) fiir 09; &. Wor To» So» to das System der 15 partiellen . 

Differentialgleichungen fir die je sechs Spannungs- und Verzerrungskomponenten und die drei 
__-Verschiebungskomponenten dar. Sie sind, und dies wird sich als die hauptsachliche Schwierig- il 
_ keit bei ihrer Integration erweisen, nicht mehr linear i in den Spannungen und NereC Or ANE Be wie — | 
dies beim Hookeschen Gesetz der Fall war. a | 


‘Wir kénnen die Spannungen aus dem Gleichungssystem eliminieren, indem wir mit den 
Plehangen (31) und (32) in die Gleichgewichtsbedingungen (9) eingehen. Es ergibt sich ise 


3K22[; x (6) + SH a) «]+26[y wa (A oe Yat 4 ie er 


dy r) ; Aesg 1 
eg sop ( we fat 3 je war ty Se vx) | +X=0. 
und ey kliceh’ weiter. Wiirden wir hier noch fiir die Cie © die partiellen Ableitungen 
der. Komponenten u, v, w des Verschiebungsvektors 8 einfuhren, so hatten wir das System von 
e TS: Gleichungen auf ein System von drei partiellen Differentialgleichungen fur u, vu, w zuriick- 
-gefuhrt, wie dies in der linearen Elastizitatstheorie iblich ist. Dieses Vorgehen empfiehlt sich in — 
ngerem. Fall jedoch nicht, ‘das besonders infolge. der komplizierten Gestalt, welche die Inva- 
riante yj annimmt, wenn man sie in den Ableitungen der Verschiebungskomponenten darstellt — 
“diese Cleichungen zu umstandlich und fir ee praktische Anwendung ungectairet, werden. " . 
e 2 

ca ‘Die Formanderungsarbeit. Una die jevehiensnee Darctclludeen der F orminderungsarbeit 
_A je Volumeneinheit herzuleiten, gehen wir aus von den Gleichungen (26) fir die Ableitungen 
der Volumenanderungsarbeit A, nach J, und der Gestaltsinderungsarbeit A’ nach ty Sie | 
_lauten, wenn wir in ihnen zemas (14) und (23) J, durch eo eal Se durch Yo ersetzen, Meet 


a : \ . Ni dA, { Vaart st salle e ‘% ‘ RM ns | 
ee RO, ee 
 woraus dureh Integration folgt Sree isn pe Ra eee yew 


q 
# 


| | 
- 


Ea 


aa | 
a | 


(40) 


lated) 


i 


/ 


Gay ic aa ae ies i cy ous Pee niien al 
ee Ao(éo) 9K fs Ex(é )de, ; A Wi= Ze J v(y?) dip? ‘ 
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ied wenn wir x(£) und »(y?) nach (27) als Potenzreihen einsetzen und die Integrationen aus- 
thren, 


1 aeaG' i 

Ag(eo) = 9K (pets met amet s+) =a Keill42 (set 4 me+--)], 
, 3 1 

Ay) =e +s mvt st) —26y 


1 
[1 +2 (7 yaya +> yas + ++) . 
Die gesamte Formanderungsarbeit' je Cee ae ist dann 
A (&q, Yo) = Ao(Eo) + A’ (yi) - 


Die Arbeit A lat sich aber auch in den Spannun 8 
gsgroBen o und T, bzw. sy und tf darstellen. 
Wegen (36) und (37) wird ja : * ; ; \ Be 


Ais,) = ox sd[sh(s)] = 9K [ [Ms) +s 2 | as, \ 


es 


Ai) = $6 [dg y= gel | g(t?) (t?) ele <8 at 
oder, wenn wir nach (38) Poialeechen fir k(s) und g(t?) einfihren und integrieren, 


9 E ) 
Adlso) = Ks [142 (Fh ot phsi+shs+--)|, 
, 3 7 
A (ioe GZ Sh [1 +2 (Feoh+eesth+tg.8 + +) : 
Hiermit wird die gesamteé Formanderungsarbeit je Volumeneinheit 
Cy _ A(Sq) 9) = Ao(S) + A’(t) - 
; Im Fall eines nichtlinearen Elastizitatsgesetzes liefern die partiellen Ableitungen von A 


nach den Spannungen o2,..., Txy,... nicht mehr, wie dies beim Hookeschen Gesetz der Fall 
war, die Verzerrungen é,,..., Wxy,,.-3 wir erhalten vielmehr nach elementarer Zwischenrechnung 


04  dAy d5 dA’ Ot” 
oy ds 0a dit 0G.,. 


a 


1 I 
BK [1 + 2h, s, + 3 a8 Sora er: 96 [1 + 3g.t +++] (oxo), 


und zyklisch weiter, sowie 
0A dA’ 0t,? 1 
eS eee ee etek Bay Tc) Tey 
- und zyklisch weiter. In den entsprechenden Gleichungen (44) und (45) fir die Verzerrungen 
éx und Waxy stehen an Stelle der eckigen Klammern die Funktionen k(sq) und g(t). 


8. Einfache exakte Lésungen der Grundgleichungen. Wir wollen uns nun den einfachsten | 


_ exakten Lésungen der elastischen Grundgleichungen zuwenden. Es liegt auf der Hand, als erste 
_ Gruppe die gleichférmigen Spannungszustande zu untersuchen. Ist der Spannungstensor un- 
-abhangig von x, y, z, so gilt dies offenbar auch vom Verzerrungstensor und umgekehrt, so daB 


die Gleichungen (9) fir R=0 und (19) erfillt sind. 
Ein leicht tiibersehbarer Fall ist der hydrostatische Spannungszustand. Wir setzen hier 


Ox=Oy=0z;=0 9; Tay = Tye = Tax = 9. 


‘Dann wird é 
1 : EG = 0, z= e== 0, ois 


2 y 

und somit — 

oe ae 

4 : = aks 0. 

ee ag Verzerrungen werden dann nach (44) und (45) 

q © 8p b= by= 82 = BIC M50) 5 Pay Pye Yow = 08 (48) 


* 


q yiir erhalten also, wie dies schon in Ziff. 5 anlaBlich der Formulierung der dritten Forderung 


= 


- 


f 


7 


a 


; auch aus (48), daB nach (31) ~ | 
Ng : Go Aa =oy = Os = SK x(e0) &0° Txy = Tyz= Vax = 0 


‘sein muB,. 


_ zum Ausdruck gebracht worden ist, eine gestaltstreue Verzerrung. Umgekehrt folgt natiirlich 
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Als nachsten Fall behandeln wir den eindimensionalen Zug-Druckversuch (in x-_ 


Richtung), indem wir 


setzen, so dah 


/ 2 / / At ; \ 
= ao Gg a. Os Oy fGen ee (49) | 
und nach (33) iS , : 
2 
ae ies o, (50) 
om o oe 
80 OR cies Oe: 
wird. Dann erhalt man fiir die Dehnung in Zugrichtung 
ik I o l 202\; ‘ 5] 4 
pate = slax' lox) to 8(sez)| > 2 ist 
_und fiir die Querdehnungen : ean i 
; 1 1 a i A 20 52 4 
== 336 * (ox) ag &(ger)| & ; be | 


wahrend die Scherungen yxy, py: und wx verschwinden. Das Querdehnungsverhaltnis m wird d 
hier eine Funktion von o, namlich 


Dono 1 20? 
as sK* lox) +78 (se) 
yore 1 oe ee ae ne 
aK" (ox) a6 & (sea ; 
Das ‘Gi peasvack hierzu bildet der Fall des in eine starre Vertiefung whe : 
Zylinders, dessen Achse zur x-Achse parallel sei. Hier ist anzusetzen 


Ex = Ee 5 by = €z= Way = Wye = Yux = 9, 


(53) 


also 
1 ; 2 y ‘ jak 
- i 3 Ne eRe Sell ae Tie en Sa 
und_nach (23) f : 212, t 
: Yo = 


‘i Far die Spannung in Richtung der Pane. erhalt man dann nach (31) 


02 =|Kx($ s)+% cy(Le)]e Te 
- und fir. die Spannungen, die auf den Zylindermantel wirken, 
: By Se Ban5 
oy=0.=|Ku (5) = 3 Gy (4 et)]e. 
Der maximale Betrag der ipa) ore ti SO_ 
} ’ | oh 78 
|t]max = Rig tle dy = Gy (> e*) : le| é 
Auch die Uberlagerung_ einer hydrostatischen. Spannung o G mit einer (etwa i in der 


x-Richtung wirkenden) einachsigen Zug- oder Druckbeanspruchung o lab ster einfach _ 
behandeln. Es ist. hier anzusetzen: pa 


} Ox=o9+0 o's fe Oy=O 2 0:=0 3 Cny = Uys Tan 0 . coat 
So wird (.- 4) 


Pat 


nia He send Acbeato sipyhannienichonpbie weeny Daily ronan 


hiss 3 (o+3.) j ; a 
und, wie beim einachsigen Zug-Deuck- Versuch, ; ee | ‘ 
su BND ck ite Bie 
RD Gorn Ge meng : 
so daB die Elafchangen fiir die Dehnungen die Form aie Ab ee 
. ied [clu ey eo RLF EC Cake hacen ieee ns ay te 
em slat oe )etse+ee(ser)o], 8h) 
Oe ee US aa o+36 20? nies Atha’ Wee 
sie firtg (ee i ( aK )(o+35 3) — ze 8 (oe) o| i) 


Pedenaey wahrend die Scherungen wiederum verschwinden. no wea jae ie 


y ( = ‘ / 
" te Pat “ = 


Die reine Schubbeanspruchung mit 
= We i Ox =Oy= O62 =Tyz=Tix =0 ’ Txy =T 
ergibt : i, 


; 2 
4 6o>=9 ’ To x. /3 T 5 ; 
also a 
j ao ot : oe 2 Tt 
Dann werden die Verzerrungen 
- 1 aimee 
Ex= by = Ex= YPys= Yrx = 0, Wxy =e8lsa) T; 
wir erhalten also eine volumentreue Verformung. : 


> Drehen wir die x- und die y-Achse um 45° um die z-Achse, so nehmen mit |t| = |o| bekanntlich 
die Spannungen die Werte an 


| 


5; é Cz—0, Oy=—O, Os Tay bys — Te ’ : a 
so daB : ee ae ‘ 
: Cp Ua te zo 
Oa oe epee (re eee i 
Bate Stes Nl ats : 1/2 «o ss 
I ara ; Be ; = —Sy=9 x = V5 e : ss 


wird. Die Verzerrungen werden jetzt . zee ; s 


ae im ane G2 Noa 
¥ Es res er OS IG 5 (5 or) &:= Puy = Y= au 
eS Der Potachote ungleichférmige Verzerrungszustand, der sich in geschlossener Form be- 
handeln laBt, ist die Torsion eines Zylinders mit kreis- oder kreisringférmigem Quer- 
schnitt vom AuBenhalbmesser_ a und vom Innenhalbmesser }20. Wir versuchen den in der 
linearen Theorie iblichen Ansatz fiir die Nereeripnseue der die Vertraglichkeitsbedingungen (19) 
erfullt, namlich — ee ee the . idee pie 


| be ty=8= Yay =O, : Yoon, } es oy, is ; (56). 


Ds 


tr aor mit (neat ni Rie = (87) 


’ : 2 : y ky? ee 


i 


Poot 7 ( 4 Ae es ? f f } Al 
der a ceenaesobane <aAuyeas Reda Sc Shc aM us 


s braucht jedoch nicht Gabemiage! wie dies bei der Me 
r Rad anzunehmen. Da a Jers nach (33) 
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a Beniitzen wir gemaB (27) die Reihenentwicklung 
ae ip) Wey Vee devaen tates : 
: so erhalten wir fir M die Potenzreihe 


out M(w) = Ge (Jpo + Jp2@? V2 + Jpao* ya + Ippo Ve + °°)» 


wobei 
Ipo = = (a*—b4) 


das polare Flachentragheitsmoment des Zylinderquerschnittes' bedeutet und die iibrigen Koeffi- — | 
zienten folgende Werte haben: 


1 8 
Jon = 75 (8-0), Ty (8B), Ing 7yg (=P). 


Wir wollen auch die im Zylinder von der Lange Eins aufgespeicherte Formanderungsarbeit A i 
berechnen. Da sie eine reine Meme a petal se erhalten: wir sie nach (47) in der Form | 


A= [ai Wo) * toed ee fra Jaye de® 


r2=b? w,2=0 ( : 
oder, da das reer: er M(o) Sores berechnet ist, hier auch durch das Integral 


as f me )dw. 


OE A eas 


In beiden Fallen wird, wenn wir r wieder die Be ete eckieoee fiir y(y>) bzw. fiir M(w) benittzen 
und die Integrationen ausfihren, 


ae) = 7 6u | Joo +> Jp2 vet T Jy4004 i es sk a rot] ; 


_ Haben wir insbesondere einen Zylinder, dessen W ahactarke h klein ist im Vemiaieuian zu 
seinem mittleren Halbmesser R, so konnen wir die Vernderlichkeit von oy) unberticksichtigt, | 
pce lassen und das Moment ansetzen zu 


| M~ 2a Rthe(R) 22 nGy (2 08h) 9 m 
oder, wenn wir hier Jpo= 2a Reh nehmen, 
M= -CIpoy (F wR?) ow. 


“Die Tinkcbvune. dieser Bezichung lautet dann CHE 
M 2- R? M? 

G dre” a(F GC Ipo? Ine) 
Bei Messungen wird man i. a. den Zusammenhang zwischen M fad 7) direkt bestimmen, 
Bei Vollstaben 146t- sich aber nach einer bekannten Formel!, deren Richtigkeit man auch al 
Hand der Darstellung (62) bis. (64) fir M und (60) far t(r) leicht unmittelbar nachprifen n 
und die bei jedem. beliebigen Elastizitatsgesetz anwendbar ist, die Speman sc. Ta am Rand des | 
6 Zylinderquerschnitts als Funktion von w ee Die Formel lautet —— a a 24 


ee sett alae gakae Wales = nar (® 

9. “Fixpertmontelle Bestimmung ce, ivitinaccsee ‘Um fir einen bestimmten 
Dua. angeben zu kénnen, ist der Wert. der Konstanten K und G sowie d I 
der F unktionen (¢9) und vy). oder k(s,) und 8(t0) experimentell in Grundversuchen zu 
Als solche eignen sich in erster Linie diejenigen, bei deren rechnerischer Auswertung 
die eine der ‘beiden zu -bestimmenden | Funktionen auftritt, also Versuche, bei denen die ¢ 
_zelnen Elemente des Versuchskérpers allein ihr Volumen oder allein ihre Gestalt ander Man 
. kommt, wenn man auf Kontrollen verzichtet, mit zwei Versuchen aus, bei seam eweils der 
maa einer mafgebenden SpannungsgréBe mi einer pec Maia 
“ _ermitteln ist. " eder der ane Deruehe liefert denn den Verlat -ein 


“< 


pares 


» 


fae 


Sn oe : est a 
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-- Von den in der vorigen Ziffer behandelten Spannungszustanden kamen demnach hier zunachst 
_ einerseits ein hydrostatischer Zug-Druck-Versuch zur Bestimmung von K und x(é) oder k(s,) 
und andererseits ein reiner Schubversuch zur Bestimmung von G und (v2) oder g(t?) in Betracht. 
Die Versuchstechnik fiir Kompressionsmessungen bei hydrostatischem Druck ist sehr gut 
_ ausgebildet, und es liegen zahlreiche Ergebnisse vor!. Die Funktionen %(€)) und k(s)) werden 
-sich somit fir negative Werte von ¢,). und s, durch solche Versuche leicht ermitteln lassen. Fir 
positive €, und sy ist dies jedoch auf einem entsprechenden Weg kaum moglich, da Versuche mit fi 
_allseitig gleichformigem Zug bis heute technisch nicht befriedigend durchfahrbar sind. be 
Ebenfalls schwierig ist es, an plattenférmigen Kérpern reine Schubversuche anzustellen, sei ee 
es, daf} man den rechteckigen Versuchskérper an seinen Randern auf Schub oder auf Zug und Beet 
gleich grofSen Druck entsprechend den beiden in der vorigen Ziffer behandelten Spannungs- ‘ 
_ zustanden beansprucht. Hier bietet jedoch der Torsionsversuch mit dinnwandigen kreiszylindri- 
schen Rohren, der ja ebenfalls eine (wenigstens annahernd) gleichférmige reine Schubbeanspru- 
chung erzeugt und sich verhaltnismafig leicht ausfithren 148t, einen guten Ersatz. Natirlich 
_ kann man auch aus einem Torsionsversuch mit einem massiven Rundstab G und y(wo) oder 
_ &(%) bestimmen; die Auswertung eines solchen Versuchs ist jedoch etwas mihsamer als beim 
_ dimnwandigen Rohr. ae : 
Ist auf Grund des hydrostatischen Druckversuchs K und x(€) fiir e950 oder k(s,) fir s) SO © ? 
und dazu mit Hilfe des Torsionsversuchs G und y(wo) oder g(tj) bekannt, so benétigt man noch epee. 
einen weiteren Versuch, um die Funktionen x(¢)) fir e97>0 oder k(s,) fiir s)>0 zu bestimmen. Bae 3 
Hierzu eignet sich am besten ein einachsiger Zugversuch, bei dem es geniigt, die Dehnung ¢éx. : 
Is Funktion der Zugspannung zu ermitteln. Nach (51) ergibt sich dabei die Funktion k(s,) mit 
Hilfe der schon bekannten Funktion g(tj) und hieraus schlieBlich auch x(é,). tel 
Es sei hier ausdriicklich darauf hingewiesen, daB der tbliche eindimensionale Zug-Druck- 
Versuch mit Messung von e;(o) allein ohne Bestimmung der Querdehnung iiber das elastische 
Verhalten des Kérpers unter einer anderen als der Versuchsbeanspruchung tiberhaupt noch keine 
ussage erméglicht. Wird jedoch die Querdehnung ey als Funktion der Spannung mitgemessen, 
so 1a48t sich sowohl K und k(s,) als auch G und g(t) auf Grund der aus (51) und (52) herleitbaren 
5. Beziehungen = is . youn : : . 


oe a ee ce 
4 ©) ae  ee)- ee) le en, 


ry 


mz 


wie ma: etaor t Eas AR eh yA a an LOSS, ae ay 1 ; 
und der Schubfunktion gewinnen kann. =  * 


W. Hort, ‘Handbucl der physikalischen ‘und technischen Me- 
‘alst wait erie Hie ay iS 4 N iS tueal } Nae: Va Anes 


ee Wes Gye | 


Ba al 


- In Zahlen wird somit hier 
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a) Mit massiven Rundstaben aus Aluminiumbronze wurde ein Zug- und ein Druckversuch 
je mit Messung der Langsdehnung é, allein sowie ein Torsionsversuch ausgefiihrt. Hierbei er- 
hielt man folgende Ergebnisse’. 


Zug- und Druckversuch: : 
to} 0 +1000 +2000 +3000 kg/cm? 
éi = 0 +0,96 + 6,25 + 36,3 - 10-3 
Torsionsversuch (Werte mit Hilfe von (65) ermittelt): 
t(a)=%_ = 0 >. 5007. -1000 1500 kg/cm? 
a0 =e = 0 “1,19 5,85 30,6 - 107-8. 

Wir wollen zuerst die Versuche analytisch ayswerten. Hiezu stellen wir ¢, und Wa als Polynome 
in o und T. dar, welche die in der Tabelle angegebenen Werte sane, Mit Hilfe des wblichen | 
Differenzenschemas? erhalten wir | 

Ex = Xo + h,0" + a3 0° 4- aon + %,0° ’ 


Ya = By Ta + Bg te + Bste 


I 


mit den Koeffizienten 
% = 0,78 10-Femtkg, a=0, ag=0, 049 - 10- 12 om®)/kg® 5° Og==0;, o,=0, 104-10" *8 em?/kg®; 
By=2,10 - 10-° em?/kg, fy=0,25 - 10-38 em9)kg®, B,—3,50- 10-18 em'®/kg®. 


Da nun nach (61) fir r=a 
Yo = =6(3%) ta = Ta + Be Ge5 Te + Bart ie tioe 
sein mu8, so folgt durch ean aces Potenzen von Ta 
C=+., & = CBs, gs = > C5B,, b= Beans Oe oe (68) 


oder in Zahlen | ; 
G = 0,477 - 10° kg/cm?, 82 == 0,040:-105c gp OF195 A022". 


Die Funktion g(tj) wird so : 
g(t) = 1-+0,040 - 106 2 +.0,195 - 10124. 


Aus dem Zug-Druck-Versuch kann man jetzt die Funktion k(sy) und K ermitteln. Wir beniitzen. 
dazu die Gleichung (51), die, nach ks) outer lautet | 


o 3 &x 1 20% 
: (ox) =3K[-F-eslse)]- 
Fuhrt man hier fir k, g und e, die Reihen bzw. Polynomdarstellungen ein, so ergibt sich 
j ere 68 oO o \4 : 
thse thse) tals) th) to 
_ 3K {3 (a4 1,003.0 + Og 0? + 045 0%) - -—— cals a poe ier we (sc oa) I} : 


a ee 
ad tae 


92 


Durch Vergleich der Koeffizienten gleicher Potenzen von o folgt hieraus, wenn man noch die 
Beziehungen (68) bericksichtigt, i: 
1 1 cs 

‘ eee =3 "ache (69) 

Ed G 5 


eens Meer as ie C7. 
ky =(9K)?a, ly =(9K)® )3 (%3— $f). a= (9K)ta45 geass a HP.) ks=hg= >: =0, (70) | 


d : 


; K=1,35- 10% kg/om?, ,=0, ky—38,0- 108, kj=0, Iy=120,4- 1008, 
so daB k(sp) die Form | : \ feces P re i 

: (so) =1+38,0 - 10%s§-+120,4 - 101254 
annimmt.: : : 


1 


ied 
ie ee) Se rWitht von M. Rof und A. Richinger, Diskussionsbericht der Eidg. Muterjalprifongsanstelth 


* Siehe z.B. G. Schulz, haiiers eras zur praktschen Mathematik. Berlin. und Leipzig 1937, Ss. 67. + 


= 
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Versucht man mit Hilfe der Umkehrung der Formeln (39) und (40) die Koeffizienten x; und 
yi der Potenzreihendarstellung der Funktionen (¢) und y(p>) zu gewinnen, so zeigt es sich, daB 
diese Reihen nur in sehr engen Bereichen von eq und yo gut konvergieren, weshalb sie fir prak- 
tische Zwecke kaum brauchbar sein dirften. 


Es wird daher niitzlich sein, wenn wir dem soeben gezeigten analytischen auch noch ein graphi- 
sches Auswertungsverfahren zur Seite stellen, das jetzt beschrieben werden soll. 


Wir gehen aus von einem Diagramm des Torsionsversuches, in dem Ta in Abhangigkeit von 
Ya eingetragen wird (Abb. 1). Die Skalen fir wa und Ta ersetzen wir zufolge der Beziehungen 


> ne 
Yo = =z Ya=0,816 pa, = \/t ta= 0;816-46 + 


die sich nach (58) und (33) fiir die Werte von yy und 7, an der Oberflache des Stabes ergeben, 
durch y, und t,. Die.Ursprungstangente mége die Gerade y,)=c, wobei c eine willkirlich in ge- 
eigneter GroBe wahlbare positive Konstante ist, im Punkte A schneiden. Dann ist die t,-Ordinate 
von A gleich cG. (Im Beispiel ist c=10-° gewahlt; dann wird cG=477 kg/em?.) Wir zeichnen 
nun auf der Ordinatenachse eine y- und auf der Abszissenachse eine g-Skala ein, wobei die Ein- 
heit fir y durch die Ordinatenstrecke und die Einheit fir g durch die Abszissenstrecke des — 

_ Punktes A darzustellen ist. AuBerdem tragen wir auf der Ordinatenachse eine t,-Skala auf, . 
deren Einheit wegen t)=1,/G gleich der mit c dividierten Ordinatenstrecke von A ist. (Im Bei- 
spiel entspricht somit diese Ordinatenstrecke selbst 10-3 Einheiten von t).) 


Um nun zu einem beliebigen Punkt P des Torsionsdiagrammes den Punkt Q der Kurve 


y(y>o) mit der Abszisse y) und der Ordinate y und den Punkt Q der Kurve g(t?) mit der Ordinate 

t) und der Abszisse g zu erhalten, zeichnen wir den Ursprungstrahl OP, der die Ordinatenparallele 
_ durch Ain R und die Abszissenparallele durch A in R schneiden mége. Der Punkt Q hat dann 
_. die Abszisse von P und die Ordinate von R, wahrend der Punkt Q die Abszisse von R und 
' die Ordinate von P hat. Der Beweis fiir die Richtigkeit dieser Konstruktion ergibt sich leicht 
B ynit’ Hilfe ‘des Strahlensatzes, wenn man bericksichtigt, daB die Abszisse t) von P nach (36) 
gleich oy (y?) und die Ordinate y, von P nach (37) gleich t)g(t}) sein mu. Man konstruiert so 
_ die Kurven y(y;) und g(t;) punktweise, wobei zu beachten ist, dafs die erstere die Gerade y=1 
_ fiir yy=0 und die letztere-die Gerade g=1 fir t)=0 berihrt, 


A, Diese Methode zur Konstruktion der Kurven fiir y und g aus dem (yo, T)-Diagramm 1]aBt sich 
a sinngema$ ubertragen auf die Konstruktion der Kurven fiir x und k aus dem (é, 0,)-Diagramm; 
 hierbei ist lediglich wy durch é , T durch oo, t) durch sy und G durch 3K zu ersetzen. 


Um nun aber im vorliegenden’ Beispiel das (€9; Go)-Diagramm zu erhalten, benétigen. wir 
~ auBer dem Diagramm des Torsionsversuchs noch das Zug-Druck-Diagramm (Abb. 2). Das erstere 

kénnen wir dabei in den Mafistabsverhaltnissen unverandert tbernehmen, wenn wir zur Ein- 

-tragung des Zug-Druck-Diagramms eine o-Skala beniitzen, deren Einheit (1 kg/cm?) die gleiche 

> Lange hat, durch welche der Wert t)= a Ye: kg/cm?=0,4714 kg/cm? auf der T>-Achse dargestellt 

_ wird, und ebenso eine ¢,-Skala, deren Einheit der Lange entspricht, durch welche / 2 Einheiten 

der wo-Skala dargestellt werden. Das Torsionsdiagramm bedeutet in diesem (Ex, 0)-System 

dann den Ausdruck | ae 5 

PO se sees: 4 2 PEO Fee | (45) ¢ 

raat : f temaee Byer 9G , 
der, von éx(0) abgezogen, gemaB (51) die Funktion 


el 
as \ 
3 (HID EN GAGS Was SUN Lot AO 
| OK I (ox) = 2X0) — se 8 (sar) ° 

liefert. Das Kurvenbild dieser Funktién laBt sich aber jetzt direkt als Darstellung von Do(Ep) 
verwenden, wenn wir die ¢,-Achse als ¢,-Achse mit der gleichen Skala bentitzen und die o- Skala 
gemaB oo = i o durch eine Skala fur o, ersetzen, auf der 1 kg/cm? durch diejenige Strecke dar- 
gestellt wird, welche 3 kg/cm? der o-Skala entspricht. Man hat also lediglich nach Einfihrung 
der Skalen fir o, ‘gs ex und €) noch die Zug-Druck-Kurve o(€:) einzuzeichnen und von deren 
A bszissenstrecken diejenigen des Torsionsdiagramms, die zur gleichen Ordinatenhéhe gehéren, 
abzuziehen. Die so entstehende’ neue Kurve stellt dann in den entsprechenden Skalen direkt 
(€5,09)-Diagramm dar. 
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‘ Abb. 1. Aluminiumbronze. Torsionsversuch, 
uh ‘Bestimmung von Scherungs- und Schubfunktion. 


eae eds BO ag ant ot eo tS 


we 


Kk Bg aie K(s 


Abb. 2. Aluminiumbronze. Torsions- und Zug-Druck-Versuch. 
ane Bestimmung der Funktion o,(&). 


X 


iy — 


\ 


ie Abb. 3 faites nun hieraus die: Diagramme’ 
far 2(€) und Ks) konstruiert, wobei das Konstruk- 
tionsverfahren, wie schon ober n-erwahnt, ganz analog 
A YEAR:: dem der Abb. 1 verlauft.. Die Abb. 1, 2.und 3 sind 
nur der - Ubersichtlichkeit wegen getrennt gezeichnet; ort 
bei der praktischen Durehfihrung_ des graphischen : cael 
Verfahrens kann man gut mit: einer ee ak Zeich- Rise ein | Rome, 
rang auskommen, oie a SS ak einer eae om 
Betrage der. Dehnungen i m vorliegenden 
viel bei Zug und bei gleich groBem Druck tiber- Ey a aoe ; 
si mmen, wurde die Darstellung der Funktionen xund k far. negative Werte von 
Beide Kurven sind hier symmetrisch zur  Ordinaten- baw. Al zissena se 


b) An. Stiben. aus Siemens-Martin-Stahl’ haben Ros und Eichinger! ie den a iiche 
( 

 Zug- und Druckversuchen die Loe ee Ex ce die Querdebaung ey gemessen, i 

eisonte ere ergaben — Mie s B Oe 


yen f Oy feta : i 
om = Rol ae ae pete. J ; wee : 
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Die noch nicht befriedigende Ubereinstimmung der Betrage von e,—ey fir gleich groBe Zug- 
; und Druckspannugen laBt sich aus dem Diagramm von Ro§ und Eichinger dadurch erklaren, 
_daB die Dehnungswerte nicht auf den spannungsfreien, sondern auf den mit ungefahr 160 kg/cm? 
_ (beim Zugversuch) bzw. mit ungefahr 100 kg/cm? (beim Druckversuch) vorgespannten Versuchs- 
_ korper bezogen sind. Beriicksichtigt man dies und verteilt sodann die noch verbleibenden kleinen 
_Differenzen entsprechender Betrage von &x—éy gleichmaBig auf die MeBergebnisse, so erhalt 
man die folgenden verbesserten Versuchswerte 
; og = +2000 +1000 0 —1000 —2000 kg/cm? Penge 
pes +1,000 +0,468 0 —0,445 —0,977 - 10-8 : arias 
fy = —0,264 —0,121 0 40,144 +0,287 « 10-8 

Ex—éy = +1,264 +0,589 0 —0,589 —1,264 . ee 


_ Man stellt We im Blick auf die Gleichungen (66) und (67) die Ausdricke €x+2ey und €,—éy 


I 


als Polynome in der Form : 


/ 


i : éyt2ey = a,0+a,0?+a,0%+a,04, &x—&y = bo +b, 03 2 (71) 
dar. Hierbei werden die Zahlwerte der Koeffizienten 
ae a,;= +0,183 - 10-6 em2/kg, ay=+0,0425 - 10-® cm4/kg?, 
eee a;= +0,009 - 10-12. cm8/kg?, a,=— 0,0085 - 10-15 cm8/kg4; 
se. ur aes bee +0,5747- 10-§ cem2/kg; 6, = +0,0143 - 10-12 em®/kg’. 
Andererseits ist nach (66) und (67) nee et cs 
[| ttm seh mak AGE) + aha) tne) 
5 Bee No 3 


sy 


72-10 


- 70 


Fey LS Py Cae eee] 
- Abb, 4. Siemens-Martin-Stahl. Kompressionsfunktion. 
3 ve Ley pale: ; , 


pegs 


; Jergleich derselben Potenzen von o in (71) und (72) die folgenden Gleichungen zur 
yung des Kompressions- und des Schubmoduls sowie der Koeffizienten der Reihen- 


hi 


ng von. Kompressions- und Schubfunktion_ liefert: Feld é a 


> a3 1 aS . f ae 
i y= 218 Ta, Kt=27-7%, kp=ke=-+- =0; eK ane 
Se een a. ae et Ah pantie) ; 
par 848s : 


: nae : eA to ¥ 
te ae Yar ieee OVS elie tag RaiRlia, eaves Mf 

Abb, 5. Siemens-Martin-Stahl. Schubfunktion. 
2+ 108s2—204,6 - 10%5 
ee ti . z a nt a 672 : ve 
g(t) =14+0,085 + 108 


d (je eils innerhalb des Versuchsbereichs) in den Ab bods eee 


: pat 13,2 


5 ay 
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: c) Vom Kupfer ist bekannt, daB es — wie zahlreiche andere Metalle — beim hydrostatischen — 
Druckversuch ein lineares Gesetz befolgt, so da man die Kompressionsfunktion h(so) fir negative 
Werte von s, konstant gleich Eins setzen kann. Der Wert des Kompressionsmoduls wird hierbei | 
K=1,33 + 108 kg/em®. Zur Bestimmung der Schubfunktion g(t?) und der Kompressionsfunktion — 
k(s,) fir positive Werte von sy steht noch ein Torsionsversuch an einem massiven Rundstab mit 


folgenden Ergebnissen* 


tT = 0.100 200 kg/cm? 
Ya = 0 0,26 0,80: 10-3 


sowie ein Zugversuch? zur Verfiigung, bei dent eine Spannung ¢ =490 kg/cm? eine Langsdehnung, 


€x=1,0- 10-% zur Folge hatte. 
Die Scherung Yo | laBt sich in Shin cele: von der Schubspannung ta durch ein Polynom | 


“Ya = By tat Bsto 
gti) darstellen, dessen Koeffizienten sich TAY 
py=2,13.-1078 cm?/kg , 
B3=46,7 - 10-17 cm®/kg? 


ergeben. Nach (68) wird dann der Schub- | 
modul 2 : 


C= a7 = 0,470 « 10° kg/cm? 
1 


und der einzige nicht verschwindende- Koeffi- | 
zient der Schubfunktion A] 


; 3 . : if | 
ee a ae se = 2.6%f, = 7,26 - 108, 
a ate Ure Oath ee also diese selbst (Abb. 6) — 
g(t) =1+7,26- 10%. 


‘ i ) , \ re 
. Um nun auch den Zugversuch durch ein Polynom von der Form 


ata Abb. 6. Kupfer. Schubfunktion. 


i = 


& =%40+ 4,05 — ; ; 


-darzustellen, bestimmen wir zundchst «,, das sich aus der Gleichung (69) zu 


NS ak oe =o 


Hies mit deri Wert %,=0, 793 - 10-6 cm?/kg errechnet. Der Koeffizient a, folgt dann aus de 
Rania da das Polynom fir e, bei dem oben angegebenen Zugspannungswert den dure! 
Versuch ermittelten Langsdehnungswert liefern muB, zu «3;=5,20- 10-12 em*/kg®. gic werde 
Sela (70) die Koeffizienten der Kompressionsfunktion k(so) 


k= 0, ky=11,979(5,20—5,19) - 108, eats eS oes BY ae 


eMiie wiehe, daB k, innerhalb der Versuchsgenauigkeit ohne Bedenken gleich Null diverse 
_ werden kann; man darf demnach | auch fiir positive Werte von fo die Bompreesiou 


(so) konstant gleich Eins setzen. ive 
10. Elastizitatsgesetze von Bach-Schiile, Griibler und Sbhleshrwep.. Wir eoeaiien unters 

inwieweit die mathematischen Ausdricke, die Bach-Schiile, Griibler und Schlechtweg 3 aufg. 
haben, um das experimentell ermittelte elastische Verhalten von, festen Kérpern zu beschr 
‘als Spezialfalle unseres allgemeinen Elastizitatsgesetzes aufzufassen sind. = 2S 
_ Bach und Schiile bentitaten,, wie enon in ler Biel part Crwanat, far ihre Zug-Dri ck- 
suche ue ee ts cae Se Soap aS 


ts aaa aA 5) aubey. in dem. erivialew Fall. m: rr nicht geniigen ‘und somit ; auch 
ne unseres: be Need a Gesetaes (51) sein kann: 


: lw Siehe FuBnote 1 os Seite re Py i Sei) * Nivea aah cist te ey 
: cans Nédai, a. a. O, S. 56, Abb. 53. _ sheng AAR tate etc as Rte a eat aes 
} ® Siehe Hapaote qe. nad 4 auf Seite 450. Pitti . ‘a Pema We came ter fae 2 


et 
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Der von Griibler zunachst fiir den Zug-Druck-Versuch beniitzte und ebenfalls in der Kinleitung 


erwahnte ,,hyperbolische‘‘ Ansatz lautet in unserer Schreibweise, nach der Langsdehnung seule 
gelost, 


oder, in eine geometrische Reihe entwickelt, 

és = Ao(1+ Bo+B?o?+--:), 
wobei A und B Materialkonstanten sind. Er stellt also tatsachlich einen Spezialfall unserer all- 
gemeinen Formel (51) dar, aus der sich jedoch natiirlich nicht die hier zugrunde liegende Kom- 


pressions- und Schubfunktion in eindeutiger Weise ermitteln lassen; denn dazu wire ja, wie 
wir gesehen haben, auBer dem Zug-Druck-Versuch noch mindestens ein weiterer Versuch notig. 


Die Verallgemeinerung des Griiblerschen Elastizitatsgesetzes in der Form 


A Tix 
Wo ore Ty. 
wobei Dy und Tx, entsprechende Komponenten des Verzerrungstensors D (12) und des Span- 
nungstensors T' (1) sind, la8t sich jedoch mit der von uns benitzten allgemeinen Form (30), 


die man sich ebenfalls in Komponenten geschrieben denken mag, nicht mehr in Hinklang bringen. 


Schlechtweg verbesserte diesen Ansatz, indem er statt der Tensorkomponenten Invarianten 


-einfihrte, und gab ihm schlieBlich — in unserer Schreibweise — mit den Materialkonstanten 
B, C und L die Form 
t 1 7 
f= 8, D'=5 : T’, (73) 


ee 9BKGs—V3Ch +913 LK VCs Vy 


- die, wie er nachwies, insbesondere dem Verhalten des GuBeisens gerecht wird. Der Vergleich 
zeigt, da zwar Gleichung (41) mit k(s,)=1 die entsprechende Gleichung des Schlechtwegschen 
Gesetzes liefert — sie stellt ja noch die lineare Hookesche Beziehung dar —, daB aber die zweite 
Gleichung von Schlechtweg sich nicht als Spezialfall von (42) darstellen laBt, da der bei T’ stehende 
Faktor auSer t, auch noch sy enthalt. Tatsachlich ist auch bei GuBeisen die Gestaltsanderung 
von einer tiberlagerten hydrostatischen Spannung nicht unabhangig, weshalb auch das von 
uns beniitzte Elastizitatsgesetz sich fiir GuBeisen nicht verwenden la8t. Man erkennt dies z. B. 
an den Ergebnissen der Zug-Druck-Versuche von Ros und Eichinger, fir die das in der vorher- 
gehenden Ziffer genannte Kriterium, wonach ¢,—ey eine ungerade Funktion von o sein mub, 
nicht erfullt. ist. 


Jetzt soll noch gezeigt Gerdes daB schon der (auch von Schlechtweg beniitzte) vereinfachte 
Ansatz, den man erhalt, wenn man die Materialkonstante L in (73) gleich Null setzt, es nicht 
2 erméglicht, einen vom Verformungsweg unabhangigen Ausdruck fir die Formanderungsarbeit A 


Materialkonstante B verschwindet, so daB der Faktor bei T’ in (73) eine Funktion von ty allein 
wird und der Ansatz fir T’ somit zu einem Sonderfall unseres Elastizitatsgesetzes wird). 


ist at B. das Bestehen der folgenden Beziehung 


EVEL ie . - (74) 
ot : 0 &y dex” 
denn beide Seiten der r Gleichong stellen ja, wie in Ziff. 7 dargelegt wurde, den Ausdruck 
; : eA 
déx 0 Ey 
dar. Wir weisen neki) ‘daB diese Beziehung beim Ansatz von Schlechtweg in der Form 
1 Pees 
eae armen wanes (75), (76) 


2 Re eae eee 


i Reha YE =. 2(C— 9BKGs,— V3 Ct) D’. 3 : . (77) 


etzt ‘sind in NS Faktor von D’ die SpannungsgrdBen Sy und ty durch Verzerrungsgré8en zu er- 
en. ae Sq ist dies bereits i in a) geschehen. Um auch fir ¢, einen Ausdruck in den Verzerrungs- 


31* 


je Volumeneinheit aufzustellen (abgesehen natiirlich von dem Sonderfall, da auch noch die 


Eine notwendige Bedingung fir die Existenz einer Funktion A im oben genannten Sinne 


a 
Away 
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gréen zu erhalten, bilden wir in (76) beiderseits die zweite Tensorinvariante und erhalten nac’ 


kurzer, aus (23), (33) und (35) folgender Plipie hace oa 


‘G = Poo 
mit der Abkirzung 
= 1—9BKe 
: Let Vi~ Yo 
So nimmt (77) die Form an 
P26 OP's 


und es wird 


Ox = 3Keq+2Gqo(Ez—&)s Oy = 3K eg +26 Hy (Ey—£>) - 
Die Ableitungen in Gleichung (74) werden mit 


B fy) Lilo A 18 Vo 8 ay 38 
park aie ca a ’ = oy — & 
Dex © O88 Mee 2 2 Yo (ex—€9) by op (€y— £9) 
zu 
0 Ox 3 0 2 2 OP 
0 €y aa Kk +6|— (€x—€p) (€y Eo) aa iG o| ap ce G(éx a) far ’ 
0 Oy hee 0 Qo 2 y IQ 
hee K +6|—-(e,—e) (€x—£9) Faperea ee ee G(éy— Eo) tect 


Gleichung (74) ist fiir beliebige Werte von ¢, und ¢y somit nur dann erfillt, wenn @ @o/d 9 ver- 
schwindet, d.h. also fir B=0. Fay 


11. Biegung prismatischer Stabe. Die Biegung prismatischer Stabe fiir nichtlineare Elasti- | 
zitatsgesetze ist schon von Bach}, Herbert? und Nddai? untersucht worden. Samtlichen Unter- 
suchungen werden die Annahmen der technischen Biegelehre zugrundegelegt, die wenigstens, 
-hinsichtlich des Ebenbleibens der Stabquerschnitte auch bei Stoffen mit nichtlinearem Elastizi- 
tatsgesetz mit geniigender Genauigkeit zuzutreffen scheinen?+. 

Man wird nun bei den Stoffen, die das hier entwickelte nichtlineare Elastizitatsgesetz be- 
folgen, nach einem Anhaltspunkt suchen, der allein aus der Kenntnis der Kompressions- und 
der Schubfunktion heraus zu entscheiden gestattet, ob man bei einem bestimmten Material damit 
rechnen darf, da die technische Biegelehre in der Form, wie sie von den oben genannten Ver- 
’ fassern fir das nichtlineare elastische Verhalten erweitert worden ist, Ergebnisse von ahnlicher 
Genauigkeit liefern wird wie die gewéhnliche technische Biegelehre bei linearem Elastizitats- 
gesetz. 

Hierbei mége uns folgender Gedankengang leiten: Bekanntlich stehen bei linearem Elastizi-. 
tatsgesetz die Ergebnisse der technischen Biegelehre i im Fall der sog. reinen Biegung (Beanspru- ~ 
chung der Stabendquerschnitte durch freie Momente) mit den Grundgleichungen der Elastizitats- | 
theorie im Einklang. Es ist nun anzunehmen, daB auch bei nicht linearem Elastizitatsgesetz | 
die erweiterte technische Biegelehre dann das gleiche leisten wird wie die gewéhnliche Biegelehre 
bei linearem Gesetz, wenn auch sie den Fall der reinen Biegung exakt wiedergibt. Wir werden 
also untersuchen miissen, unter welcher Bedingung fiir Kompressions- und Schubfunktion ein 
Spannungszustand als Folge einer reinen Biegebeanspruchung moglich ist, wie ihn die technische | 
Biegelehre fiir den Fall der reinen Biegung bei nichtlinearem Elastizitatsgesetz annimmt. : 

Hiezu-gehen wir von einem noch allgemeineren Spannungszustand aus, der den Fall der 
_ technischen Biegelehre sicher umfaBt; wir setzen namlich voraus, daB bei fehlenden Volumkraften 4] 
der Spannungstensor T nur die einzige Komponente o,=o enthalte, die jedoch zunachst noch © q 
von allen drei Ortskoordinaten abhangen darf, und nehmen an ~ 


Ox=6 (x; > z), Oy =O2= Tay = Tys= Tix = 0, g= 0. d ' : , (78) ay 
Man sieht sofort, daB die Gleichgewichtsbedingungen (9) nur dann erfillt werden kénnen, wenn. 


do,/0x verschwindet, also o nur von y und z abhangt. Dann worees mit den eee in (49) und | | 
(50) fir konstantes o bestimmten aon ae F- 


- 2 “ay 
Ree = 4 Oly 2), ty = 2 oly, 2)’ 


1 Siehe FuBnote 1 auf Seite 450. ; 

2 H. Herbert, Forsch, Arb. Ing. Wes., Heft 89 (1910). 
8A, Nédai, a.a.O., S.114. . 

4 C. Bach und R. Baumann, a. a. O., S. 267. 


- ih > an 7 


XVII. Band 1949. -Kauderer: Uber ein nichtlineares Elastizitatsgesetz. _ 471 


— 


-_ entsprechend (51) und (52) die Verzerrungen . 
ete fod o 1 20? : 
(= lag k(ax) tq e(ser)|7- eng I) 


9 
erly) = ely.) = 3 [ae k (Ge) — ge a (GS)] 0. (80) 
Pry = Pyz = Yux = 0. 


Diese missen die Vertraglichkeitsbedingungen (19) erfiillen, die sich mit Ricksicht darauf, 
dafi die Dehnungen nicht von x abhangen, auf folgende Gleichungen reduzieren: 
Oty oe Siac ea 0? ex ; 0? 0? 
ay? =0, dy 02 =u) =O) As,=A 6.=0 mit A= ata: 


Da samtliche zweiten Ableitungen von ¢, verschwinden, kann ¢, nur eine lineare Funktion von y 
und z sein: 


&x = Ay +Bz+C (A,B,C = konst). (81) 


«Es gilt also notwendig auch Aeé,=0. Somit miissen die beiden Gleichungen. 


Aes = Fe Ao +5 [(22)'+(22)'] =0, Shea 


Oy dz 
oo Aey = Gedo + 52 | (32) +(32)'] =0 


gleichzeitig erfillt werden. Dies ist aber nur méglich, wenn entweder 
; do\2 OG\ 2H 
Ge) +Ge) =% | 
also g=konst und somit der Spannungszustand gleichformig ist (welcher Fall uns hier nicht mehr 
interessiert) oder 


_ also éy eine lineare Funktion von é, ist, falls man o als Parameter fiir die beiden Veranderlichen 
CO und éy auffaBt. Da aber fir o=0 nach (79) und (80) sowohl ¢, als auch ey verschwinden, muB 
notwendig Pie 
emote 1 
éy =*konst SOx ee a ek 
At ee eine : y ™po ' 
werden, wenn man entsprechend der Querdehnungszahl m in der linearen Theorie auch hier eine 
_Konstante my einfihrt. Ein Spannungszustand der Form (78) ist also mit nicht konstantem o 
héchstens dann méglich, wenn die Querdehnungsfunktion m(¢) in (53), die man ja durch einen — 
-einfachen Zug-Druck-Versuch auch direkt ermitteln kann, eine Konstante ist. 
_. Wie man aus (53) erkennt, wird aber m(o) nur dann eine Konstante, wenn fir alle o die 
nktionen g und k sich héchstens um einen konstanten Faktor unterscheiden, der dann aber 


ach (38) nur die Zahl Eins sein kann, also wenn 


Be Bax) oe Gel 
Die Dehnungen werden in diesem Sonderfall : 

ee e=a(se+c)*(sx)% | eek ae 
eel 


nd die ; Gleichung fir ex lautet, nach o aufgeldst, 
oe Boe ee ORG: oh Cee 3s 
Re rae, o—sez6 * (aKa) © } cian 
-den somit annehmen diirfen, daB die technische Biegelehre bei nichtlinearem Elasti- 
tz um so bessere Ergebnisse liefern wird, je geringer die Schwankung der in (53) dar- 
en Funktion m(o) innerhalb des in Frage kommenden Bereichs von a ist. Im allgemeinen 


Fa 


ein o gerade Funktion m(o) dieser Forderung wenigstens fir kleine Werte von o besser 

erecht werden als eine ungerade, die man dann erhalt, wenn k in o lineare Glieder oder, was 

s gleich bedeutet, &x in o quadratische Glieder aufweist, falls man sich die rechte Seite von 
nach Potenzen von o entwickelt denkt. Peay ei 


’ 


os ww % 
Nae Sates z 
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Far Aluminiumbronze nach Beispiel a) von Ziff. 9 erhalt man etwa den in Abb. 7 dargestellten 
Verlauf von m(o). Hier kann man fiir |o] < 1000 kg/em®* noch gut mit einem konstanten Mittel-. 
wert 2,95 fiir m rechnen. 

Nun wollen wir zeigen, wie man die exakten Differentialgleichungen fiir die reine Biegung 

des prismatischen Stabes. bei nichtlinearem Elasti- 

m0 C7 zitatsgesetz aufzustellen hat, und sodann deren In- 

tegration fiir eine besonders einfache Form dieses 

Gesetzes naherungsweise durchfiihren. Wir gehen — 

von der Annahme aus (deren Richtigkeit sich erst 

am Schlu8 der Untersuchung erweisen wird), dab 

26 o simtliche SpannungsgréBen von. der Abszisse x (in 

7500 1000 $00 0 500 1000 «7500 Kg Richtung der Stabachse) unabhangig sind und die 

ye Schubspannungen Txy, Tx (und somit auch die 

ste Ss ue a eee Scherungen Wey, Wsx) verschwinden. Wir setzen also 
die Spannungen an 


PA Ox(¥s2)+ Gy(¥s2)s O2(Ys3)s Ty2(¥s3), Tay Toe =O. - (83) 
Y Die Gleichgewichtsbedingungen (9) werden bei fehlender Volumkraft samtliche befriedigt, wenn _ ‘| 


wir dy, 6; und Tye aus einer Spannungsfunktion F(y,z) in der Form | 

: eF ar SOs a a8 Be eae 
SO EN One ae 63 — qy , Cyne Oy oz (aay | 
herleiten. Wir erhalten dann nach (3), (33) und (35) | 

a ~ io 

: 1 : sy 

: $9 = BR eT (89); | 
Be Ne Tag cee ge ke ye ak \2 Cg sta ae eee 2). + 
Case ite | (3) : 2 sya) +(3) p hoes A eee (86) | 


pio’ 


_ Greift im Fall der reinen Biegung an den Stabendflachen Q ein freies Moment M an, dessen Vektor | is 
die Richtung der y-Achse haben mége. so verlangt man bekanntlich, daB 


JJ 0% dyd z=0, JJ ony dydz=0 SS Ont dy dz= ee WR (87), (88). (69) 


wird. AuBerdem muB die Mantelflache oo Stabes eee sein; d. h. es muf dort gelten, a 
falls « der Winkel zwischen der (nach auSen weigenden) Flachennormalen und der Richtung der 
We -Achse ist 


dy cosa. ty: sina = 0, Tyz cosa. + a.sing = 0. i Ss gee 


Drticke man diese Forderungen in der Spannungsfunktion F(y.z z) aus, so erhalt man nach einer 
it _ Umformung, die vollig analog zu derjenigen verlauft, die fur die Spannungsfunktion des ebenen 
_ Spannungszustandes fiir den lastfreien Teil des Randes durchgefihrt zu werden PiRese das 


_ Ergebnis, da in den Randpunkten 
Tee et . 
anes 5‘ : (s— == konst, 


RP ; 
) es te = konst 
fn) yY Rand ~ 


02 ary 


sein muB, wobei die Konstanten gleich Null gewahlt werden dirfen. ‘Wenn wir nun mit (8 ) 
Pe (86) i in die Gleichungen (44) und Vag fir die ae eingehen sO stellen sich diese 


; die Vecigienliche x ¢ enthalten. Mit. Bsckuchs fiierauf sowie auf Pee Metschwindon, der: Sohavonben ¥ | 


Yxy und l wx reduzieren sich die von den. Verzerrungen zu erfiillenden Vertraglichkeitsbedingungen 4 
aaemeas : RNG ine <i eae 


Ey 


oh ‘Die ersten drei Gleichungen verlangen Wieder) daBe &x eine lineare Funktion von bg und a sein | 
i mul, also wie in (81). ° 
ee Ay +Bz+C (4, oe - konst). 


; In ae poise ausgedriickt, liefert dies eine Berichung. nwischen A, B, Gs o mee 
zweiten Ableitungen von LF, welche die Form hehe ent x 


aF BRL 98 
fi (2 ad ar 


: ; ay’ ' i= OLY az : tr) = Ays BAC ek % 
a i ‘Siehe Zs ‘Be L. Féppl, Drang und Twang, 3. Bd. ‘Miinchen_ 1947, 1B. a) ve ; ; ; | Ler: i 


a ¢ f Meta Tn Wes oe eS 
t oo . 83 2 cal fie Vath ay ? * i 
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‘annehmen wird. Die vierte: Vertraglichkeitsbedingung (91) ergibt, wenn man sie in den Span- 
‘nungen ausdriickt, eine partielle Differentialgleichung, welche die Ableitungen von F'(y, z) bis 
zur vierten und die von ox bis zur zweiten Ordnung enthalt, von der Form 


: Ia. 80, ao, 2B 3F 46 
if. (3 i 0" 0, K 0 i fr) é Ok }=0. (94) 


ee aa ey yee ye aaa Tye” 


Grundsatzlich kénnte man nun ox und seine Ableitungen aus (93) und (94) eliminieren und erhielte 
so eine partielle Differentialgleichung vierter Ordnung fir F(y,z z), die y, z und die zweiten, dritten__ ney. 
und vierten Ableitungen von Fy, z) enthalt; ihre Lésung, in der noch die Konstanten A, Bund C ie 
stehen, ist dann so zu bestimmen, daB sie die Randbedingungen (90) befriedigt. Hierdurch sind 
i. a. die zweiten Ableitungen von F'(y,z) als Funktionen von y und z eindeutig festgelegt. Geht 
man hiermit in (93) ein, so erhalt man ox als Funktion von y und z sowie den drei Konstanten 
A, B und C. Diese reichen nunmehr gerade noch aus, um die Bedingungen (87) bis (89) zu be- 
friedigen, und somit erweist sich unsere Annahme (83), von der wir ausgegangen sind, als ge- 
rechtfertigt. Hierbei kénnen sich noch sémtliche drei Konstanten, A, B und C verandern, wenn 
‘man den Wert von M variiert, und dies bedeutet praktisch, daB Lage und Richtung der neutralen 
Achse Ay+Bz+C= 0 vom Wert von M abhangt. ies 


Wir erkennen also aus dieser allgemeinen Untersuchung, daB im Fall der reinen Biegung 
bei hichtlinearem: Elastizitatsgesetz mit folgenden, in allen Stabquerschnitten gleichen Span- 
ungen zu rechnen ist: Normalspannungen o, (auf Ebenen senkrecht zur Balkenachse) sowie 
Normalspannungen Gy, 6; und Schubspannungen 1,, (auf Ebenen parallel zur Balkenachse), 
die samtliche nicht mehr linear verlaufen werden. Papa. treten dann folgende, ebenfalls. von der 


mehr see ZU Ex ed aad Haber Se: Yyz Bur Folge haben. 
i . { 


Die es des eee age (93), eet und ‘die ae der Lésung an die 


‘18 - == | Ag (A= konst) 


t. Die Konstante ) wird also hie identisch mit 8». Dabei sei noch zusatzlich : angenom- 
} nur wenig vom Wert Eins abweicht, d.h. da Glieder von der GréBenordnung — 
-solchen von der Pits At a cera werden dirfen. 


Bk wind” afer F( Gee 0 wid wenn wir voraussetzen, daB die yehchae ae 
tod piebaneay ata) ist, Ox= Fe z, wobei . 


Lad 


L Botts pecwten: der y-Achse | bedeoutet, Wir op 
‘rion Ay: 2) und Fly z) an ye eats ye if 


\ ms + 


Foy ay ae saute ces ae (85) 09, 
pied , “ £ , wf pS 3 ‘ye fod. Mi A : ' . 5 wae et 
a Tye = z Th ay ae oe Be ri a i ett oe Ua joe 


 Oaddzate der. mit» z 'hobebieten Korrektionsglieder_ 
f ssigt Werdew) Fabre ‘man die cee AUnE. unter diesem 
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Nach den Grundgleichungen (44) und (45) werden dann die Verzerrungen 


= nae AO) tog [2peta(2x—-AS+ oer B=") |. (98) | 


by = gic [Bet Az+4f)] + oe [Beta (8-2E 4 — Af - 8 pst), (99) 


y f 
es = ge [Bet M+ 4] tog [- 644 (35; sq B%2*)]» (100). 
at OE 
| Vat e: ae Pax = Yay = 0- ase 
Da der Wert von &z bei Hookeschem Elastizitatsgesetz sich zu 
byl 1 
en = 3 (¢ + aK) 
ergibt, empfiehlt sich hier fir ez der wegen (92) in y und z lineare Ansatz 
Ex = Exy tA (ay+bz-+c) 
mit den Konstanten a, b und c. Mit ihm wird aus (98) 
Posed 1 1 Dane te ; 
(+3) *— (ae ~ aK) 4I+ gs BOs = Flay tbe te). See 


Ebenso wird mit (99) bis (101) aus der Vertraglichkeitsbedingung (91) nach einiger Umformung | 
die Gleichung si 


~ (gh — Alar +(e + pi) ats~-deptens. 0) 


Um die Funktion y aus (102) und (103) zu eliminieren, wenden wir auf (102) die A- Since 
an und erhalten : 


1 I 1 1 4 , LAY 

Ge fae) t Ge ae ee (104) | 
Siirashier! man jetzt Ay aus (103) und (104), so ergibt sich fir f(y,z) die Differentialgleichung jf 
8 6° | q 
Hate 3G(3K+46) ~ rue Ue) if 

deren Lésung die aus (90) folgenden Bedingungen . 
of of | | 

Cees. (So) oe (106) ; 


erfiillen muB. 


Hat man die Lésung (105) gefunden, so ist oy, ¢, und tyz nach (97) berechenbar. Die Funktion i 
x(y.z), die man zur Bestimmung von o;, bendtigt, es aus (102) zu | 3 
3KG 1 
108) = see (se a +x) Af — 525 82*+3( ay-+bz-+o)] . aon 4 


Zur Bestimmung der drei Konstanten a, b, c stehen jetzt noch die drei Forderungen (87) bis i 
' (89) fiir o, zur Verfigung, die sich fir x(ys2) in der Form a 


SS ulyz) dy dz=0, Sf xly.z)ydydz=0, Sf y(y,z)z dy dz=0- (108), (109), (i1oy'§} 
eX) 0 ooet Ae: 
' ausdricken. ‘* 


Als Anwendung ‘bohaiidela’ wir die Biegung eines zylindrischen Stabes mit- -elliptischem | 
Sea oe Rand durch die Gleichung Hee Rope x | 


y2+022—1? = 0 a 


mit beliebigem r und beliebigem +0 bestimmt ist. “Dic Losang der Gleichung eee mit den on | 
Randbedingungen (106) 1aBt ee in diesem Fall leicht erraten; sie ist 5 


" flys2) = Cz(y? +8222 18)? 


mit einer Saoeh zu bestimmenden Konstanten C. Die ersten Ableitungen 


7a =4Cyz(y 24022212) , } if = Cly* +5 940471 Gort) 
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__verschwinden offensichtlich auf dem Ellipsenrand. Die Spannungen werden dann nach (97) 
Oy =ACz- 497 (3y24+5 9222—3r2), o, = ACz- A4(3y?+ 022212), 
Ty: = —ACy + 4(y?+3 22272) , 
Af=4Cz [3(1+ 0?) y2+92(1+5 02) 22—(143 7) 72] (111) 
AAf= sh Le aa a 


so da man, um (105) zu befriedigen, 


und es wird 


sowie 


B oS 
= C= Fate Esa CEKTIG,. iy) 
- zu setzen hat. Nach (107) wird dann 
2 ae Ae 216° K (ay-+ba+o) 6G @K— 26) Cz [3(1+0%) v2.62 (145.0%) 22 (1+30%)r9] 2 K p28 
: 3G62(3K+G6) 
_hierbei sind noch die Konstanten a, b und c zu bestimmen. 
; Aus Symmetriegriinden muf nun beim elliptischen Querschnitt o, eine gerade Funktion von 
_y sein, woraus sich nach (95) dasselbe auch fir x(y,2) ergibt. Da gemals (111) Af in y gerade 
ist, kann aber, wie Gleichung (107) zeigt, y(y,z) nur dann in y gerade sein, wenn a verschwindet, 
“was zur Folge hat, daB dann Gleichung (109) erfillt wird. Nunmehr ist. y(y,z) nach (113) als — 
umme aus einer in z ungeraden Funktion und einer Konstanten, welche c als Faktor enthalt, 
darstellbar. Bei der in (108) verlangten Integration kann somit wegen der Symmetrieeigen- 
schaften des Ellipsenquerschnitts der Integralwert nur dann. verschwinden, wenn auch ‘die 
Konstante ¢ zu Null wird. Um die noch fehlende Konstante 6 zu bestimmen, steht die Gleichung 
(110) zur Verfiigung, die sich nach (113) mit a=c=0 auch in der Form 
jaded nites +)P—32 2 - 
tye Se f ere*K6246643K- -2G) C2[3( 1488) y2-4-02(14599)24— (1+3 8?) r?|— 
pene ee . — 2K 6323\ zdydz = 0 


Edicethen laBt. Die Auswertung dieses Integrals fir den Ellipsenquerschnitt gibt, wenn man 
den Wert. von C nach (112) einsetzt, fiir b die Gleichung 


(113) 


(333 r2 
| = STEre 
Geht man mit , diesem ‘Wert in die Gleichung (113) fur x ein, so erhalt man 1 nach (95) fiir die 
- Bspeennune:: Og = Ba tA Cz(Cyy?—C,22+C,r?), 


obei die’ Koeffizienten Cy, C,, Cs die Werte — 


(14-82) (3 K—26) 
oe ee 
3 oe p MUFSO) KES 4419415 09 KE4 245 09,0 
ue Ties on é G(3K+6) f 
Coe C 904204509 HO E+ROLTI KEL 4909 #0 
2 Uy oa See PCGKLG) | 


: IM poosertional dem Mice ec eeu ist, wird der Faktor C nach (112) proportional 
n Potenz des Biegemoments, ' und dasselbe gilt somit von den mit dem Faktor AC 
fte n Anderungen samtlicher Spannungen sowie vom Faktor b, der ja die Anderung der 
mung der r (sich i in einen Kreisbogen verformenden) Mittelfaser y=z=0 bedeutet, die wegen 
wie vor die ‘Eigenschaft besitat, das ihre Punkte spannungsfrei bleiben. 

| pezialfall einer kreisrunden Welle v vom Halbmesser te wird mit C= 1 und bis = pre 


a 


semis mg 
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’ Eine genauere Untersuchung dieser Spannungswerte zeigt, da® der groBte Betrag, den eine 
Normalspannung einer Schnittflache parallel zur x-Achse annehmen kann, an den Stellen ies o 
z= +rvon der Spannung dy erreicht wird. Es wird dort 

1 M3 1 ; 
OS ecg pe Sea eee) A 


Fihrt man (als Vergleichswert) die Spannung 


oy(9, 


M 
One eres yo a 


ein, die sich in den au®ersten Fasern des Querschnitts bei linearem Elastizitatsgesetz unter der 
Wirkung des Momentes M ergeben wiirde, so wird an den Stellen y=0, z= +r das Verhaltnis 


nur® 


obo. 9 GBK+L4G) 


Bei positivem A(g,) haben somit oy und 04, gleiche, bei negativem A dagegen verschiedene 
V orzeichen. 
An den gleichen Stellen nimmt der Spannungswert 


Oy 1 oo y 


00% = Ox— Zs 


M 
Pie 
um den sich die Langsspannung gegeniiber ihrem Wert bei linearem Elastizitatsgesetz andert, 
sein Extremum an. Man erhalt dort auf Grund von (114) das Verhaltnis 


do, 1 18K®+21KG6+2C? Soy mae! 
Gore 8 26(K+G) G(3K+46) ae 
Bei positivem / wird also der Betrag von oy, in den auBersten Fasern verkleinert, bei negativem A 


vergroBert. 
So ergibt sich z. B. fir G4, = 100 kg/cm? mit den Werten 


x A= gen 1,c0- 10°)" 605470 10" ke] em Base 106 kg/em? 
fiir uae nach Ziff. 9, suse c) As 
Sh te eee 9Y _ 9.992 - 10-2, SOK 8 yg Ege 
esi F ea Sbo 
‘Die zusatzlichen Spannungen halten sich demnach hier in sehr niedrigen Grenzen; es ist jedoch , 
zu beachten, da®B sie mit der dritten Potenz des. Biegemomentes anwachsen. 


Nees be po Neca aes ei pak 


Da auch der Beton ein nichtlineares Elastizitatsgesetz befolgt, ist es denkbar, da in Eisen- 
betontragern bei starker Biegebeanspruchung eine zusatzliche Zugspannung oy senkrecht zur 
Eiseneinlage auftritt, welche die relativ niedrige Bruchspannung des senkrecht zum ‘Trager 
nicht durch Eiseneinlagen verstarkten Betons erreicht und so Langsrisse auf der Oberflache, 
des Tragers hervorruft. : : BC wee 4 es 


12. Torsion yon Stiben mit beliebigem Querschnitt. Als weiteres Problem és Elastizitats 
theorie mége noch die Torsion eines zur z-Achse parallelen prismatischen Stabes von beliebigem 
panetscanity Q unter einem Verdrillungsmoment M behandelt werden. 


Wir wollen untersuchen, ob sich mit dem in der linearen Elastizitatstheorie ablichen Ansata 
‘far die Spannungen und die Volumenkraft ~ 


: Gioia, cg Vcianctayls eRe a 
Sher gesigneter Wahl yon Ty: und tx die Gleichgewichts- und die Vertraglichheitsbedingunge 7a {| 
auch. hier erfillen lassen. Die ersteren, die sich nach (9) auf : 


put, sen, ae Tx 0 Tyz ‘ f us bat foe Ye xu 
ax PAA es Oy : 4c SATS eat ty oe tik 


5 ‘reduzieren, werden genau dann erfullt, wenn man die Spannungen gemab 


eho yeaa @ 
: lin Se oa es hs 


von einer ‘Schubspannungsfunktion © (x,y) herleitet. Auf der Mantelflache des Seabees: roialh deme 
bekanntlich @=const. sein, und man wird i. a., da es ja bei der Bestimmung von Ty: und Tx. n 


auf die Ableitungen von ®(x, y) ankommt, auf einfach vusammenhangenden puuese 
-pandern ®=0 annehmen. : 


ay ang 
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In O(x,y) ausgedriickt, wird 


= $= ser((32) +35) 


so da} nach (44) und (45) die Verzerrungen 


1 0@ 


Sie Oe 1 ID 
Ex= Ey = Es= Way =0, Pys= — & 8 (t) ax” Pax = G 8(t0) 


ay 


(115), (116), (117) 


werden. Die Vertraglichkeitsbedingungen (19) reduzieren sich dann auf 


0 (9Yse 8 Pye) _ 9 (9x Oy 
eae Pe Slay are 
oder 
0 Pex 0 Pyz i 
| heel, Se (c=konst ), 
was wegen (116) und (117) der Gleichung _ 
Ry en ON is 8 Cie Oe 1 
es [a SS | +55 [2 l= (118) 
gleichwertig ist. Fihrt man die Differentiationen aus, so ergibt sich mit den Abkirzungen 


0? 
dy? 


folgende quasilineare partielle Differentialgleichung zweiter Ordnung fir die Spannungsfunktion 


D(x,y): 4 D2 xO ID aD Yo ® 2@ 
2). 1/42), 0 af dP 2 Sh dla eal = 
&(to) A®D+g' (tt) 3@ (S>) ax ax dy dxdy ova ry | ao: 


Zur Berechnung der Verschiebungen versuchen wir, wie in der linearen Theorie, den Ansatz 


(= gas), 425+ 
& oo de 8 0)? cae 


(119) 


‘2 U=WY2, V=—WX2, w=w(x,y), 


wobei w den Torsionswinkel je Langeneinheit der Stabachse bedeutet. Dann sind offenbar die 
Gleichungen (115) mit (10) und (11) erfillt. Aus (116) und (117) folgt 


aw Ee pg Oe — 8w ekg le @ 

a GE, Pe gy 2 OC Bag 

und somit : 
dw 1 (eo et dw 1 0 
eon a an OO Ot 


© als System von Differentialgleichungen zur Bestimmung der Verwolbung w(x, y)- Leitet man die 

_ erste. der beiden Gleichungen nach x, die zweite nach y ab und setzt die rechten Seiten einander 

' gleich, so erhalt man 

1 () a® 0 aD 

sae os [seo $2] +4 [ee 4$]} --20 

oder wegen (118) ea f 
Scho 

ro OED Tlie 

also ein Ergebnis, das ganz dem bei linearem Elastizitatsgesetz entspricht. 

Das Verdrillungsmoment M ergibt sich aus dem Integral 


-M=S/S (tyx—txy) dx dy 
Q 


unter der Voraussetzung, daB auf dem Rand O=0 ist, in gleicher Weise wie bei der linearen 
Theorie zu ges . 

M=2/f (x,y) dxdy. (120) 
j . t Q i 
Das Torsionsproblem besteht somit im wesentlichen darin, Lésungen der quasilinearen. par- 


E tiellen Differentialgleichung zweiter Ordnung fir (x,y) (119) zu finden, die auf dem Rand des 
- Querschnitts verschwinden. Da diese Aufgabe jedoch infolge der komplizierten Form dieser 


in der vorigen Ziffer behandelten Sonderfall entwickelt werden, daB die Funktion g(¢j) in der 


Se aa itda © , (121) 


_ Differentialgleichung schwierig ist, soll hier wiederum eine Naherungsmethode fir den schon | 
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darstellbar ist, wobei At} < 1 angenommen werden dae Die Forderaae k(s)) = 1, die beim | 
Biegeproblem der vorigen Ziffer noch zusatzlich zu stellen war, braucht hier jedoch nicht er- 
hoben zu werden. : . 4 
Wir nehmen an, das Torsionsproblem sei fir das lineare Elastizitatsgesetz (mit A=0) fur 
den gleichen Verdrillungswinkel @ schon gelést und habe die auf dem Rand varschindeie 
Spannungsfunktion (x,y) ergeben, welche die Gleichung : 

A® =c=—2Cw (122) | 
befriedigt. Die zu dieser Funktion gehérigen Werte von Ty:, Tx und ty seien ebenfalls ‘durch 
Uberstreichen gekennzeichnet.-Wir setzen dann die Lésung des nichtlinearen Problems in der 
Form 


D(x,y) =D (x,y) +Ap(%,y) 
mittels einer zu bestimmenden Funktion p(x,y) an, die auf dem Quersesrar penne verschwinden 
mu. Offenbar wird dann 


is a® a ae P : ; 

te op he Ie (123) 
DBE de ee 4g BG Ot 

Seige Fairs aie Tox ta ae (124) | 

gu 2) (fee agp \? a® ap\?] = 4 (3D dp , © dy « 

i= gels tae h te + ase) | oA aes Pee ay as eee 

und ‘ eo: 

‘ : AD =AD + IAAyp=c+Adg. . Sake | 


i SEE aus (118) = 
Nidief eas af oP ap 2 » (9® aP\)} 
am +48) ( Pot ara dy a4 Ati) ( paeae a hy 

oder, wenn man von den mit A behafteten Gliedern nur die linearen bericksichtigt, 

: = a sa®-, | aw a @ -, 

46+4| 55 (G8 + $2) +3, Gy tay)l = 

Nees (122) folgt hieraus fir 0) die lineare partielle wisi cleaNeo zweiter Oesrtung 


7. 

Tes ‘a® an? a ane 
A 0 
PS Si Ie" Py ay 


7 


==), * . (125) 4 | 


Als Anwendung bexerdan wir die. Verdrillung eines Stabes von elliptischem Ovcrochtene & 
mit oder Randkurve Mall 


x? by 2_ p=) “(r, 0 —konst). 
Die de St. Venantsche Losung fir das lineare Elastizitatsgesetz ist bolanmehen 


D(x) == rage (t+ Oty? ys 


woraus ca Teh a ? yiant fie: a 
ae aD Cx 2 aD 1 COE a | 
8 es Ee cae ; y 148?” : . 
ee one E 9 ty? 4c? penne re Aer ‘Oy = 4 
i sie 2 é 2 PAS ; 0 et. 
a On Ge (148)? (x? 404 y we je 7 3@ +e) ; ay = 3EF aro aie 
und somit nach us) fiir vl, i die Digs ebiie aces terete cig ay NO er lee Sie arial 
o 4 A rey Ag +3@ ae : a [(3+.0%) x24. 94(143 #4) 1] = at 0 Rice, 
_folat Die res dee ) die auf dem Rand verschwindet, mul, wie man sich leicht berlegt, ott 
die ‘Gestalt Gib . aah | | 


23 1 
g(a iy) = ~ 3e@ +a (x? Lory 272) Jat byt yr?) 
omit den drei Konstanten ap und as annehmen. Es wird namlich dann 


Chea A . : 
ee : 36 nye lore yet ol Bae noone yt aL-a p+ ae syle 4 


ee 
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S und man berechnet hieraus durch drei lineare Gleichungen fiir x, 8 und y, daB die Differential- 
: gleichung fiir (xy) befriedigt wird, falls man 


pete Dee SUE ee Egan OAD OE SOF ‘ 
a ee a Boo erte | : 
_ 1 341662410044 16.98+3 98 | (27) | 
12 (140%) (+6 0-- 0%) 


* -setazt. Die Spannungen werden nach ee bis (124) 


2G 
Tys = Te we{1 pee Ei ay eee. aR [2a a2 + (9 a+B)y (o= zy) rl}. (128) \ 3 
2682 16 
_ bax Deatin 1+ wy) Aw? 369 (1+ 6) [(d2a-+8) yu? 42 B O24 2 ( (B—0? y) )r]}. 


Fir das Torsionsmoment M erhalt man durch Integration aber die Ellipsenflache Q gemaB 
ae (120) nach einiger Zwischenrechnung a 


= 


MS 25/0 x, ¥) tedy —21 ff Gtay) de ay +S oven dats 


aGwr* h a 2 4G+20°+354) , pa | (130) 
ee 2 00TE) 1 oO” 90 +08 al 


Fur den Fall des Kreisquerschnitts vom Halbmesser r(@=1) wird aes eis so daB an ~ 
der Stelle SBT Y= Oe die; SRE aNE. Tyz nach (128) den Wert 


Gor (1-10? - r?) 


annimmt,-der sich aus dem exakten Wert (59) offenbar dann ergibt, wenn man die Potenzreihen- 


entwicklung der Funktion Ly (F 7) 27?) nach dem zweiten Glied abbricht und deren Koeffizienten y,, n 
der ja nach (40) mit dem negativ genommenen Koeffizienten g, der Reihe fur g aibereinstimmt, : a 
durch (— A). ersetzt. Im gleichen Sinn erkennt man auch die Ubereinstimmung des Moment- aes Sea > 
werts oe von ee mit Be oueen von gage (62). 


chnitts far das Hookesche Gasets hue bekanntlich die Eigenschaft, da sie auf Ellipsen, 
die mur Randellipse des” Querschnittes ahnlich sind und 4hnlich liegen, konstante Werte an- 


immt, s SO. da i)  Blciohrertig die ee ee fir SO as darstellt, die von 


ae ‘Dann anid T=10 cm, j= a ond: wir erhalten nach (2%) 


ok TS ar og BES ee 171 IE 
Be Coe B= aE pan eqea 


oS 


den Konstanten fiir Siomens:Martn-Seaht A Bt 
gs C= 0, 870 - - 108 kg/em?, Ai Bois = 0, 085 - 106, . 
Dail > » von Lift 9 erhalten. haben, wird somit bei einer Verdrillung 


as wee 10-t emt ee SS aa . 


6112 
~~ 15375 was 
= 1920 0,135) kgjem? ©. 


ade : me . a 


2 


ate 
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Die maximale Schubspannung verringert sich demnach gegentiber dem Wert, den sie bei gleichein: 
Verdrillungswinkel und gleichem Schubmodul G bei Hookeschem Elastizitatsgesetz annehmen 
wirde, um 13,5°% und das Torsionsmoment ebenso um 7,1 %. 


13. Schlufbemerkung. Im Zusammenhang mit dem vorstehenden Versuch, die Grund- 
gleichungen einer Theorie des elastischen Verhaltens fester Stoffe bei kleinen Verzerrungen in 
nichtlinearer Form aufzubauen, erdffnen sich im wesentlichen zwei groBe Gebiete, die eingehen- 

{ der zu durchforschen waren. So ist es einerseits nétig, die gebrauchlichen Werkstoffe daraufhin 
zu priifen, wieweit sie das hier entwickelte Elastizitatsgesetz befolgen. Hierzu waren vor allem. 
hydrostatische Spannungsversuche, Torsionsversuche sowie Zug- und Druckversuche an Staben 
mit gleichzeitiger Messung von Langs- und Querdehnung durchzufiihren. Die Verzerrungsgr6Ben 
sind dabei in erster Linie fir die Spannungen unterhalb der jeweils technisch zulassigen Héchst- 
werte zu ermitteln, da ja das Elastizitatsgesetz gerade fir die Anwendung in diesem Span- 
nungsbereich bestimmt ist. 
Bernd Andererseits miBte die Theorie der nicht elementaren Belastungsfalle, die ja hier nur durch _ 
die zwei Beispiele der reinen Biegung und der Torsion prismatischer Stabe vertreten ist, noch weiter _ 
entwickelt werden. Weitere geschlossene exakte Lésungen fiir nicht gleichférmige (aber praktisch 
vorkommende) Spannungszustande werden sich wohl kaum noch finden lassen; hingegen kann 
man die Behandlung von Spannungszustanden, die nur von einer Ortskoordinate abhangen’ 
(z. B. solchen an dickwandigen Rohren unter Innen- und AuSendruck oder an umlaufenden 
Scheiben mit Rotationssymmetrie) i. a. auf eine gewohnliche Differentialgleichung zurtick- _ 
fihren, die fir beliebige Kompressions- und Schubfunktionen eine rechnerische Losung nach 
Differenzenverfahren oder eine geeignete zeichnerische Lésung zulaBt. Allerdings wird die—4} 
Anpassung der Lésung an die Randbedingungen wegen der komplizierten nichtlinearen Form 4 
der Differentialgleichungen ziemlich groBe Schwierigkeiten bereiten. nee 
-Zweidimensionale Probleme werden sich jedoch wohl nur noch far bestimmte einfache For: 
men der Kompressions- und der Schubfunktion lésen lassen. Dies gilt z. B. vom Problem des 
ebenen Spannungszustandes, bei dem sich die partielle Differentialgleichung fur die Spannungs- 
-funktion leicht angeben laBt, wobei aber eine Weiterbehandlung dieser Gleichung nur bei kleiner 
Abweichung vom Hookeschen Gesetz (wie beim Problem der reinen Biegung) méglich zu sein 
ae : 
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